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Foreword 



This volume in the Academic Press Electromagnetism series 
presents a self-contained introduction to magnetic resonance imaging 
(MRI). This technology has been evolving at a remarkable pace and it 
has become a very valuable and even indispensable diagnostic tool in 
various areas of medicine. This remarkable expansion of MRI technol- 
ogy has been achieved as a result of coordinated progress in such areas 
as physics of nuclear magnetic resonance, invention of novel image 
acquisition and image reconstruction techniques and development of 
sophisticated MRI instrumentation. Thus, it is becoming increasingly 
important to give the exposition of MRI technology with the emphasis 
on its interdisciplinary nature. This is exactly what the author of this 
book has accomplished. 

This is a short and concise book that nevertheless covers an extra- 
ordinary amount of technical information. It is intended for readers 
without extensive experience in the area of MRI technology. It reviews 
the basic underlying physical principles of nuclear magnetic reson- 
ance by and large within the framework of classical electromagnetism. 
The book also covers the main MRI imaging techniques with the 
emphasis on such topics as image contrast, signal-to-noise ratio, 
image artifacts, rapid imaging, and flow imaging. The final chapter 
of the book contains a brief discussion of MRI instrumentation. 

I believe that this book will be very attractive as a clear and concise 
introduction to the MRI technology. As such, it will be a valuable refer- 
ence for beginners and practitioners in the field. Physicists, radiologists, 
electrical and biomedical engineers will find this book very informative. 

Isaak Mayergoyz, Series Editor 
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Introduction 



Since its first implementation by Lauterbur [1], Magnetic Reso- 
nance Imaging (MRI) has become an important noninvasive imaging 
modality. MRI has found a number of applications in the fields of 
biology, engineering, and material science. Because it provides unique 
contrast between soft tissues (which is generally superior to that of 
CT) and high spatial resolution, MRI has revolutionized diagnostic 
imaging in medical science. An important advantage of diagnostic MRI 
as compared to CT is that the former does not use ionizing radiation. 

MRI is based on the phenomenon of Nuclear Magnetic Resonance 
(NMR), independently discovered by Bloch et al. [2] and Purcell et al. 
[3]. Although the potential of NMR spectroscopy was recognized 
almost immediately and its development started soon after the discov- 
ery of NMR, it took more than two decades to implement NMR-based 
imaging. Even after the pioneering work by Lauterbur and the develop- 
ment of basic imaging techniques by Kumar et al. [4] and Mansfield [5], 
several more years were required to design and develop imaging hard- 
ware at the level necessary to produce high-quality diagnostic images 
of the human body. Despite its relatively slow beginning, MRI has 
become an indispensable diagnostic tool since the early 1980s. 

The foundation of the NMR phenomenon is the interaction 
between an external magnetic field and nuclei which have a nonzero 
magnetic moment. According to the classical theory of electromagnet- 
ism, the motion of individual nuclear moments in a static magnetic 
field B is a precession about B at an angular frequency io Q , known 
as the Larmor frequency, which is proportional to the strength of the 

1 



magnetic field. Another well-known fact is that the energy of inter- 
action with B depends on the direction of nuclear moments in such 
a way that the minimum energy corresponds to the state in which 
the moments are parallel to B . As a result, in thermal equilibrium 
the majority of nuclear magnetic moments are aligned along the 
external field. The alignment of magnetic moments gives rise to non- 
zero magnetization in macroscopic samples of solids, liquids, or gases 
containing a large number of nuclei (e.g., ~10 23 cm" 3 of hydrogen 
nuclei in water). 

The NMR phenomenon is observed when a macroscopic sample in 
a static magnetic field is irradiated by an oscillating magnetic field of 
frequency u> that equals the frequency of precession, uj . The NMR phe- 
nomenon can be best explained by using the arguments put forward by 
Bloch [6]. Suppose that a macroscopic sample is placed between the 
poles of a magnet that produces a static magnetic field, B . Under 
the influence of B the sample becomes magnetized. In thermal equili- 
brium the nuclear magnetization in the sample can be expressed as 

M = X Bo, (1) 

where \ is the nuclear susceptibility. Suppose that the sample experi- 
ences an oscillating magnetic field B 1 produced by an alternating cur- 
rent in a coil encompassing the sample, and assume that the oscillating 
field is perpendicular to B . It can be shown that at resonance, 
achieved when the frequency of B^ equals w , even a weak oscillating 
field applied as a pulse can rotate the magnetization in the sample 
and place it in the transverse (i.e., perpendicular to B Q ) plane. 1 After 
the excitation pulse ends, the transverse magnetization in the sample 
precesses about B . During the precession the transverse magnetiza- 
tion decays because of nuclear interactions and nonuniformity of B 
(see Chapter 1). From Faraday's law of induction it follows that the 
time-varying magnetization induces voltage in the coil. The induced 
voltage, which identifies the presence of the transverse magnetization 
in the sample, can be Fourier transformed to obtain a NMR spectrum. 
For many applications of NMR spectroscopy it is critically important 
that the observed signal includes contributions from nuclei in different 
chemical environments (e.g., hydrogen nuclei in water and fat). Such 
nuclei normally have slightly different frequencies of precession 
depending on the chemical composition of a sample. NMR spectra 



1 In contrast, off resonance excitation normally produces negligibly small components of 
the magnetization in the transverse plane (see Chapter 1). 



can be used to identify chemically different populations of nuclei and 
to determine their relative amounts in a sample. 

We are indebted to F. Bloch for a set of widely used phenomeno- 
logical equations that describe the dynamics of nuclear magnetization 
in media [6]. In particular, Bloch has predicted that thermal fluctua- 
tions would cause exponential relaxation of the longitudinal (i.e., 
parallel to JB ) magnetization, M ton , to the equilibrium state in equa- 
tion (1): 

M lon - X B oc expH/TO, (2) 

where 7\ is a characteristic relaxation time. 

On the other hand, Bloch has suggested that magnetic inter- 
actions between neighboring nuclei cause exponential decay of the 
transverse (i.e., perpendicular to B ) magnetization, M tr , with a char- 
acteristic time constant T 2 : 

M fr «expH/T 2 ). (3) 

Although the Bloch model has important limitations (e.g., it fails to 
describe NMR in solids), it has proven to work well for liquids and bio- 
logical systems in general. Throughout this book the Bloch equations 
(discussed in more detail in Chapter 1) are used as the main tool for 
analysis of the NMR phenomenon and associated effects important 
for MRI. 

The pioneering work of Lauterbur [1] has demonstrated that the 
NMR signal acquired in the presence of external magnetic field gradi- 
ents can be used to obtain MR images of the transverse magnetization 
in the object. The basic effect that makes MRI possible is that in the 
presence of magnetic field gradients the frequency of precession, uj q , 
becomes spatially dependent. The observed NMR signal is a sum of 
many signals produced by nuclei at different locations in the sample. 
Each component of the signal acquired in the presence of magnetic 
field gradients is characterized by its unique frequency and phase. 
How can we produce an image from the signal that is given by a sum 
of different frequency components? To answer this question we need 
to review the basic concepts of spatial encoding and image reconstruc- 
tion in MRI. 

It can be shown (see Chapter 3) that the fundamental relationship 
between the NMR signal, S, and the transverse magnetization in the 
imaged object can be expressed as 

S(ft)oc / M tr e jkr dV. (4) 



In this equation k is a spatial frequency denned by external magnetic 
field gradients used for imaging, and the integral is taken over the 
object's volume. This equation is extremely important for MRI because 
it shows that the signal is given by Fourier transform of transverse 
magnetization in the object. At this stage we need to recall the well- 
known property of Fourier transform: if f(k) is the Fourier transform 
of f(x), thenf(x) is the inverse Fourier transform of /(&). After examin- 
ing equation (4), it becomes apparent that an image of the transverse 
magnetization in the object can be produced by computing inverse 
Fourier transform of the signal, S. 

In general, spatially encoded signals in MRI are obtained by using 
a number of repetitive excitations of the nuclear magnetization in the 
object. After each excitation the NMR signal is sampled a number of 
times (typically 256) during a short acquisition interval (limited by 
the decay of the transverse magnetization) in the presence of external 
magnetic field gradients. Most of the complexity in the reconstruction 
of MR images exists due to the fact that in practice the continuous 
NMR signal can only be represented by a series of its discreet samples. 
The problem of reconstructing a continuous function from its discreet 
samples is well known in signal processing. In particular, the famous 
sampling theorem [7] defines conditions under which a continuous 
function can be completely recovered from its samples. Principles of 
signal sampling in MRI as well as the relationship between imaging 
parameters (including strength and duration of imaging gradients, 
number and interval between excitations, etc.) and such important 
factors as scan time, image resolution, image contrast, and signal-to- 
noise ratio are discussed in Chapters 3 through 5. 

Hydrogen nuclei ( 1 H) appear to be the best target for in vivo MRI. 
They have the following advantages: a) among all nuclei present in 
tissues hydrogen nuclei produce the greatest NMR signal; b) X H MRI 
in vivo achieves excellent contrast between different tissues. The 
remarkable feature of *H MRI is that it provides a variety of ways to 
manipulate image contrast. Since the inception of MRI it has been 
apparent that the major source of contrast in MR images is the depen- 
dence of 7\ and T 2 relaxation times, and proton density on the bio- 
chemical composition of different tissues. The studies of Damadian 
[8] and other investigators [9] have demonstrated that in certain 
instances malignant tissues have longer 7\ and T 2 relaxation times 
than normal tissues. Because of the exponential dependence on the 
tissue relaxation times, the NMR signal is extremely sensitive to 
small variations in 7\ and T 2 . Moreover, by adjusting acquisition par- 
ameters (see Chapter 4) the dependence of the signal on the relaxation 



times can easily be varied in order to improve image contrast. For 
many applications of MRI it is important that contrast between 
different tissues can be further enhanced by using certain materials 
that can significantly shorten 7\ and T 2 in tissue. It is noteworthy 
that in some cases administered materials are preferentially absorbed 
by malignant tissues making it possible to improve identification of 
malignancies in MR images. Different mechanisms of contrast as 
well as options for manipulating contrast in MRI are discussed in 
Chapter 4. 

When acquiring MR images the major concern is the relative 
amount of noise present in the signal. The signal-to-noise ratio (SNR) 
is one of the most important image characteristics because it defines 
the observer's ability to distinguish between different structures in 
an image. Noise in NMR can be considered as a component of the 
signal produced by randomly fluctuating currents in a receiver coil 
and imaged object. It is interesting that the relative importance of 
these two sources of noise depends upon the strength of the applied 
magnetic field: at low field strengths the noise from the coil dominates; 
at high field strengths the noise from the object becomes more 
important (see Chapter 5). In principle, the signal-to-noise ratio in 
MRI depends on two different kinds of imaging parameters: intrinsic 
parameters (e.g., 7\, T 2 , proton density), which are typically beyond 
our control; and user-controlled parameters (e.g., spatial resolution, 
acquisition time, and signal averaging), which can vary depending 
on the goals of a particular MRI study. For example, in many instances 
high spatial resolution might be necessary to identify small structures 
in images. On the other hand, since the SNR in MR images decreases as 
spatial resolution increases, the quality of images acquired with high 
spatial resolution might be very poor. As a result MR imaging 
frequently entails a compromise between the conflicting requirements 
of high spatial resolution and high SNR. Image SNR can normally be 
increased by using signal averaging. The main limitation of signal 
averaging is that it increases scan time in MRI. Moreover, signal 
averaging is relatively ineffective because the SNR increases slowly 
with the number of times the signal is averaged. Therefore, it is 
important to know if and how SNR can be increased by varying 
other imaging parameters (e.g., matrix size and acquisition time). 
The SNR dependence on basic imaging parameters and sampling 
schemes used for acquisition of MR data is discussed in Chapter 5. 

A well-known problem in MRI is the existence of various artifacts 
that can significantly degrade MR images and can mistakenly be iden- 
tified as real structures in the object. Artifacts in MRI are most often 



caused by magnetic field inhomogeneities and different physiologic 
motions of patients during scanning. 

Despite the fact that when designing and manufacturing MR scan- 
ners special care is taken to ensure a high degree of homogeneity of the 
main magnetic field B , in practice magnetic field inhomogeneities are 
always present due to spatial variations in the magnetic susceptibility 
of the object. Magnetic field inhomogeneities cause two different kinds 
of artifacts in MR images: geometric distortion and signal loss. Geo- 
metric distortion results from the presence of intrinsic magnetic field 
gradients that alter the dependence of Larmor frequencies of nuclear 
spins on their spatial locations. Signal loss is caused by dephasing of 
nuclear spins in a nonuniform magnetic field. Image distortions and 
signal loss can vary from mild to severe depending on the degree of 
magnetic field nonuniformity, imaging parameters, and technique 
used to acquire MR images. A detailed discussion of these MRI arti- 
facts is presented in Chapter 6. 

One of the fundamental assumptions of MRI is that the intrinsic 
dynamics of transverse magnetization in the object does not change 
during image acquisition. In reality this assumption can often be vio- 
lated as a result of different motions in the imaged system of nuclei. 
Examples of such motions, important for in vivo MRI, include blood 
flow, cardiac and respiratory motions. Motion-induced artifacts can 
vary in appearance. One example is a series of equally spaced artificial 
structures produced as a result of periodic motion (e.g., blood flow 
pulsations). Another example is image blurring caused by the motion 
of the object as a whole. Chapter 6 contains a discussion of different 
motion artifacts and various techniques used to reduce these artifacts 
in MR images. 

Analysis of spatial encoding in MRI as well as practical experience 
clearly indicate that motion artifacts can be reduced by shortening 
scan time. Historically, the need to decrease motion artifacts was 
one of the main factors responsible for the development of rapid 
MRI. It is interesting that one of the fastest acquisition schemes in 
MRI, known as echo-planar imaging, was proposed by Mansfield soon 
after the inception of MRI [5]. However, it took more than a decade 
to develop special hardware and to improve the data acquisition 
scheme in echo-planar imaging before it became widely used. The 
main advantage of echo-planar imaging over the conventional 
approach is that the former, in principle, allows acquisition of a 
whole image after a single excitation of the transverse magnetization 
in the object, making it possible to obtained echo-planar images in a 
fraction of a second. Such short scan times are achieved at the expense 



of prolonged signal acquisition. Unfortunately, increased decay of the 
transverse magnetization resulting from long acquisition time often 
causes significant artifacts in echo-planar imaging. An alternative 
approach for rapid MRI (based on a more conventional acquisition 
scheme) achieves short scan times by reducing the interval between 
successive r.f. excitations. The latter approach has become commonly 
used because it does not require special hardware and imposes less 
stringent requirements on magnetic field uniformity than echo- 
planar imaging. This and other techniques for rapid MRI are consid- 
ered in Chapter 7. 

Because this book is intended primarily for those without exten- 
sive experience in MRI, its main goal is for the reader to develop an 
understanding of MRI and its applications. To achieve this goal it 
seems appropriate to begin with a review of the basic physical princi- 
ples of NMR, including relaxation of magnetization in media, spin 
echoes, and attenuation of the NMR signal due to diffusion. In Chapter 
1, discussion of these subjects is based on the classical theory of elec- 
tromagnetism with the addition of certain concepts from quantum 
mechanics useful for understanding the NMR phenomenon. Excitation 
of nuclear magnetization in a sample is the subject of Chapter 2. Chap- 
ter 3 considers principles of spatial encoding and the conventional 
acquisition scheme used in MRI. The material from the first three chap- 
ters is used in Chapters 4 through 8, which discuss particular aspects 
and applications of MRI including image contrast, signal-to-noise 
ratio, image artifacts, techniques for rapid MRI, and flow imaging. 
Finally, Chapter 9 contains a discussion of the main components of 
MRI instrumentation including magnets, gradient coils and radio- 
frequency coils. 
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CHAPTER 1 



Basic Principles of Nuclear 
Magnetic Resonance 



1.1. THE PHENOMENON OF NMR 

The phenomenon of Nuclear Magnetic Resonance (NMR) was 
independently discovered by two groups of physicists headed by 
F. Bloch and E.M. Purcell [1-4]. The basic physical effect at work in 
NMR is the interaction between nuclei with a nonzero magnetic 
moment and an external magnetic field. The NMR phenomenon is 
observed when a system of nuclei in a static magnetic field B experi- 
ences a perturbation by an oscillating magnetic field. The frequency 
uj of the oscillating field must satisfy the following condition: 

hu =\Ei-Ey\. (1.1.1) 

In this equation E t and E? are two Zeeman energies of the magnetic 
interaction between a nucleus and B . To determine the resonant 
frequency in (1.1.1), recall that the interaction is described by the 
Hamiltonian 

H = -jiB , (1.1.2) 

where \i is the magnetic moment of a nucleus. The moment \i can be 



\i = ihi, (1.1.3) 

where M is the angular momentum of the nucleus and 7 is a constant 
called the gyromagnetic ratio. The energy of interaction in (1.1.2) 
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Figure 1.1. Zeeman energies for the case / = 1/2. 



depends on the direction of the magnetic moment relative to the 
applied field. For example, in the case when nuclei have spin 1/2 (see 
Figure 1.1), the allowed energies are given by 



(1.1.4) 



J -lhB /2, if n is parallel to B 

\ ihBo/2, if n is anti-parallel to B . 

In a general case the Zeeman energies can be 
E m = —lwyhBQ, where m = —I, — I + 1, . . . , I. From the perturbation 
theory it follows that an alternating magnetic field can cause transi- 
tions between states m and m, only if m - m = ±1 [5]. In other 
words, only transitions between adjacent Zeeman levels separated by 
energy 6E = jhB are permitted. This result makes it clear that the 
resonant frequency in Eq. (1.1.1) is given by 

u = 6E/h = 7-B - (1-1-5) 

It should be noted that for typical magnetic field strengths used in 
NMR imaging and spectroscopy (~1 Tesla), the resonant frequencies 
belong to the r.f. region of the electromagnetic spectrum (see Table 
1.1). 

In NMR experiments an applied r.f. field induces transitions 
between adjacent Zeeman energy levels. As an example, we consider 
hydrogen nuclei with spin / = 1/2. In thermal equilibrium the 
probability of finding a nucleus in a particular state with energy E m 
is given by 

^eM-EJkT)^ (116) 

where Z = Y^m ex P( - ^m/^) i s known as the partition function, k is the 
Boltzmann constant and T is the temperature. Because Ei/ 2 < E_y 2 , 
the equilibrium population of hydrogen nuclei (/I1/2) with the lower 
energy E x / 2 exceeds the population of nuclei (^-1/2) with the higher 
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Table 1.1. NMR constants for several nuclei used in biological 

applications of NMR 

Element Spin Gyromagnetic Natural Resonance frequency, 

ratio (s _1 T" 1 ) abundance v = u/2ir (MHz) 
(%) at 1 Tesla 



X H 


1/2 


2.675 x 10 8 


99.98 


42.573 


13 C 


1/2 


6.726 x 10 7 


1.11 


10.705 


19 F 


1/2 


2.517 x 10 8 


100.0 


40.052 


2a Na 


3/2 


7.077 x 10 7 


100.0 


11.263 


31p 


1/2 


1.083 x 10 8 


100.0 


17.237 



energy E_i/ 2 . Consequently, there are more induced transitions from 
the lower energy level to the higher energy level than the reverse. 
Although the difference in population of the energy levels is very 
small (e.g., (n 1/2 - n_ 1/2 )/n 1/2 « jhB /kT ~ 10" 5 at 1 Tesla), the large 
number of nuclei participating in the induced transitions between 
these levels makes it possible to observe absorption of r.f. energy in 
macroscopic samples of solids, liquids or gases. Since the probability 
of induced transitions reaches its maximum when the r.f. frequency 
reaches jB [5], we can expect maximum energy absorption at the 
resonant frequency. 

Purcell, Torrey and Pound [4] studied absorption of r.f. energy in 
paraffin containing hydrogen nuclei. In their classic experiment the 
frequency w and amplitude of an oscillating field 3i were fixed while 
the strength of a static magnetic field B , applied perpendicular to 
Bi, was gradually varied. A sharp peak in the absorption spectrum 
was indeed observed when the strength of B reached u>/j. 

We shall see later in this chapter that the NMR signal from a 
macroscopic sample is proportional to the equilibrium magnetization 
in the sample M = ra/(|i), where n 7 and (ji) are the concentration and 
average magnetic moment of nuclei, respectively. Because the energy 
of magnetic interaction in (1.1.2) reaches its minimum when n is 
parallel to B , the probability that a nuclear magnetic moment is 
parallel to B is greater than its being antiparallel to B . Consequently 
in thermal equilibrium more magnetic moments in a sample are aligned 
with the external field than against the field. On the other hand, 
because the energy of interaction does not depend on the component 
of ft in the plane perpendicular to B (referred to as the transverse 
plane), the average transverse component of n is zero. It is therefore 
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clear that the average nuclear moment (ji) is parallel to the applied 
field. Because M is proportional to (|i), the equilibrium magnetization 
satisfies the following conditions: 

M z = M Q > 0, M x = and M y = 0, (1.1.7) 

where the z-axis is chosen in the direction of B . The equation describ- 
ing magnetization in a system of nuclei with a spin /and gyromagnetic 
ratio 7/ can be written as 

(1.1.8) 



3kT 

assuming that the energy difference between adjacent Zeeman levels is 
much smaller than thermal energy of nuclei, that is, hfjB /kT «d 1 [5]. 
Note that this condition is normally satisfied in NMR spectroscopy and 
imaging experiments. 

As we will see in Section 1.4, the observed NMR signal is created 
by the transverse component of nuclear magnetization, M tr (M tr J. B ), 
in a sample. It is therefore important to verify whether the transverse 
magnetization is indeed excited by a r.f. field applied in the presence of 
a static field B . We can prove excitation of M tr in the case when an r.f. 
field B x is perpendicular to B by using the equation for the absorbed 
r.f. power P: 

P = - f M^idV. (1.1.9) 

Since it has been established that P is nonzero during r.f. irradia- 
tion [4], we can deduce that the transverse component of magnetiza- 
tion must somehow be excited by the irradiating field. We shall 
consider the dynamics of spins during excitation in the following 
section. 

1.2. MOTION OF MAGNETIC MOMENTS 

Although only quantum mechanics can completely describe the 
NMR phenomenon, some of its features can be explained by the classi- 
cal theory of electromagnetism. In the classical theory the motion of a 
magnetic moment \i in an external field B is described by the following 
equation: 

dfi/dt = 7 fi x B (1.2.1) 



1.2. Motion of (i in an External Magnetic Field 



k(B -co/y)+iB rf 




Figure 1.2. (a) Precession of the magnetic moment n about B in 
coordinate system; (b) precession of (i about the 



rigure i.z. ^a; rrecession 01 un 

the laboratory coordinate system; yu) precession or (i aooux; ine 

„«?„„+; +:„ «„u wd , jj^ _)_ iB r{ in the rotating coordinate 



effective magnetic field k(5 ( 
system. 



In the case of static B , a solution of the last equation can be written 



(i x (t) = (i x (0) cosu t + fi y (0) sinu> t, (1.2.2a) 

H y (t) = n y (0) cosco t - n x {0) smu> t, (1.2.2b) 

^(0=^(0), (1.2.2c) 

where co = jB and the 2-axis is chosen along B . According to 
(1.2.2a-c) ]i precesses about B on the surface of a cone at a Larmor 
frequency 7B0 (see Figure 1.2(a)). 

In the presence of an r.f. field B x of frequency to, the motion of |i in 
the laboratory coordinate system is described by the equation 



dn/dt = nx>y(B + B 1 ) 



(1.2.3) 



To obtain a solution to Eq. (1.2.3) we consider a coordinate system 
rotating around the z-axis with an angular frequency — to. The relation- 
ship between the components of an arbitrary vector a in the laboratory 
coordinate system (a x ,a y ,a z ) and the rotating reference frame 

(a x > , ay , a z i ) is given by 



a x i -- 



a x cos u:t — a y sin uit, ay = a x sin uit + a y cos cot, a^ - 



(1.2.4) 
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We now wish to consider a case when the field B : is linearly polar- 
ized and perpendicular to B . That is, 

B lx = B l cosujt, B ly = 0, B lz = 0. (1.2.5) 

The linearly polarized magnetic field in Eq. (1.2.5) can be considered as 
the sum of two circularly polarized components rotating around the z- 
axis at the same frequency but in opposite directions. We shall focus on 
the effect of the component rotating in the direction of the Larmor pre- 
cession in the laboratory coordinate system, because the counterrotat- 
ing component only slightly perturbs the motion of nuclear magnetic 
moments and can therefore be neglected. By using Eqs. (1.2.3)-(1.2.5) 
we obtain the following equation of motion in the rotating reference 
frame: 

dii/dt = fi x 7[k(B - w/7) + ifi,.f.], (1-2.6) 

where i and k are the unit vectors in the x and z' directions in the 
rotating frame of reference, jB rf = B 1 /2. According to Eq. (1.2.6) the 
dynamics of n is defined by the effective static magnetic field 
B eff = k(B - uj/-y) + i-B r .f.. Therefore, in the rotating reference frame, 
(i precesses about B eff at an angular frequency 7[(Bo — ^ll) 2 + ^r f ] 
(Figure 1.2(b)). 

Equation (1.2.6) greatly simplifies the analysis of the motion of \i 
by eliminating the time dependence of the r.f. field. Notice that at 
resonance {uj = jB ) the effective magnetic field in the rotating 
reference frame equals iB r{ . To consider the effect on resonance r.f. 
irradiation we assume that initially (i.e., before the beginning of the 
irradiation) ji is parallel to the static field B . During the irradiation 
the magnetic moment will precess about iB r{ in the y'-z' plane in the 
rotating reference frame at an angular frequency -)B ri . For example, 
a 90 degree rotation of the magnetic moment will be completed at 
time t = tt/2^B i f . when the moment is along the j'-axis. If the r.f. irra- 
diation ends immediately after the completion of the 90 degree rota- 
tion, then \i will precess about B in the transverse plane in the 
laboratory coordinate system. Such excitation is known as a 90 
degree pulse. An excitation pulse, which rotates \i by 180 degrees in 
the rotating reference frame, is referred to as a 180 degree pulse. 
Notice that at resonance even a relatively weak r.f. field can cause 
rotation by an arbitrary angle in the y-z' plane. 

Because nuclear magnetization M is proportional to the average 
magnetic moment (\i), Eq. (1.2.6) can be used to describe the dynamics 
of M in a sample under the assumption that internuclear interactions 
during r.f. irradiation can be neglected. Based on the observation that 



1.3. The Block Equations 



B eff makes an angle 6 with B such that tan 6 = B r{ /(B Q - w/7), we find 
that in the case when B T f <C \B — uj/j\ nuclear magnetization initially 
parallel to B will remain aligned along B with very small transverse 
components of M produced during excitation. It then follows that the 
resonance condition lo = 7_B is necessary for effective excitation of 
the transverse magnetization by a relatively weak r.f. field B r{ <C B 
normally used in NMR imaging and spectroscopy. 

1.3. THE BLOCH EQUATIONS 

The classical model of motion of a free magnetic moment (see 
previous section), although useful for understanding of the NMR 
phenomenon, cannot explain many of its important features defined 
by the interactions between nuclei. To overcome some of the limita- 
tions of the classical model, F. Bloch introduced in 1946 the phenomen- 
ological equations describing the dynamics of nuclear magnetization 
which have become an extremely useful tool for theoretical analysis 
in NMR imaging and spectroscopy. 

The foundation of the Bloch model [3] can be understood by con- 
sidering relaxation of nuclear magnetization M in a sample after an 
excitation r.f. pulse. We already know that in thermal equilibrium, 
M lies in the direction of an external static field B . Bloch suggested 
that the equilibrium state is established because of two different 
processes governing the dynamics of M: thermal perturbations and 
internuclear interactions. In the Bloch model [3] it is the thermal 
perturbations which cause relaxation of the longitudinal magnetiza- 
tion M z (assuming that the z-axis is taken along the direction of B ) 
to its equilibrium state with the minimum energy of spins 
E = -B M, while interactions between nuclei in the sample cause 
decay of the transverse magnetization M tr = {M x ,M y } without affect- 
ing E. 

According to reference [3] the relaxation of M z is described by the 
equation 

dMJdt = (M - M z )/T,, (1.3.1) 

where T l is a constant known as the spin-lattice relaxation time. A solu- 
tion of this equation can be written as 

M z = M + [M 2 (0) - M ] expH/Ti), (1.3.2) 

where M z (0) is the initial value of M z . Relaxation of M z after a 90 
3 pulse is shown in Figure 1.3(a). 
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Figure 1.3. (a) Relaxation of M z after a 90 degree pulse; (b) decay of M tr 
after a 90 degree pulse. 



In a system of identical noninteracting nuclei placed in a static 
uniform magnetic field, nuclear magnetic moments would precess at 
the same Larmor frequencies. In reality different nuclear moments 
precess at slightly different Larmor frequencies as a result of mag- 
netic interactions between neighboring nuclei. The presence of micro- 
scopic magnetic fields created by the nuclei themselves or by 
electrons causes dephasing of nuclear spins in a sample. This results 
in an exponential decay of M tr that occurs with a characteristic time 
constant T 2 known as the spin-spin relaxation time (Figure 1.3b). The 
dynamics of M tr in the laboratory coordinate system is described by 
the following equations: 
dM x /dt = ~M x /T 2 + ^M y B , dM y /dt = -M y /T 2 -~/M x B . (1.3.3) 

Equations (1.3.1) and (1.3.3) are known as the Bloch equations. By 
examining (1.3.3) we obtain the following solutions for M x and M y : 



M x = exp(-Z/T 2 )[M,(0) cos(u t) + M y {0) sin(w 0L 
M y = exv{-t/T 2 )[-M x {0)sm(u t) + M y {0)cos{oj t)\. 



(1.3.4) 



where ui = -yB , M x (0) and M y (0) are the initial values of the com- 
ponents of M tr . To express Eqs. (1.3.4) in a more compact form let 

M xy = M x +jM y , (1.3.5) 

where j = -y/^l. Using M xy we can rewrite Eqs. (1.3.4) as 

M xy = M xy (0) expHu; « - t/T 2 ), (1.3.6) 
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where M xy (0) is the initial value of M xy . This equation describes 
both T 2 decay and precession of the transverse magnetization about 
B . It should be noted that the actual value of T 2 (as well as T x ) depends 
on the properties of nuclei and surrounding media, and the strength of 
B . It can be shown that, in general, T 2 < T\. The quantum-mechanical 
theory of nuclear relaxation is quite complicated and beyond the 
scope of this book (see [5,6] for an in-depth discussion of relaxation 
mechanisms and references to the relevant literature). 

The Bloch equations can be modified in order to take into account 
the presence of an oscillating magnetic field Bj. Everywhere in this 
book we will assume that Bj is much smaller than B . Under this 
assumption the Bloch equations in the laboratory system can be writ- 
ten as follows: 

^ = 7 Mx(k 5o + Bl )-M tr /T 2 -kfc^. (1.3.7) 

1.4. BASIC NMR EXPERIMENT 

To avoid r.f. interference, detection of the NMR signal is typically 
performed after excitation of the transverse magnetization in a sample. 
According to the principle of reciprocity [7], the induced emf in the 
receiver r.f. coil is given by 

In this equation B l /I c is the magnetic field produced by a unit current 
in the coil at the location of M. The integral in Eq. (1.4.1) is taken over 
the sample's volume. For simplicity we consider a long cylindrical coil 
of length L with N turns encompassing a sample of volume V. In this 
case B 1 /I c = fi (N/L), and the r.f. field is parallel to the axis of the 
coil. We further assume that the coil's axis is perpendicular to B . 

Based on the results obtained in the previous sections we can now 
describe a simple NMR experiment in which a system of nuclei is 
initially excited by a short 90 degree pulse with duration t p <g:T 2 ,T 1 . 
Immediately after the pulse the nuclear magnetization lies in the 
transverse plane. Using Eqs. (1.3.4) and assuming that the Larmor 
frequency of nuclei w > 1/T 2 , we obtain 

emf = HoNL' 1 Vu; M e~ t/T2 cos(ui t + ip), (1.4.2) 

where the phase ip is dependent upon the excitation r.f. pulse. This 
equation shows that the induced signal, known as the free induction 
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decay (FID), is proportional to u> M . Since lo M oc n^ 1(1 + l)B , the 
sensitivity of signal detection in NMR experiments can be increased by 
using high field strengths and exciting nuclei that have large 7 and are 
most abundant in the sample. Because J H nuclei satisfy the latter 
conditions in biological systems, they appear to provide the highest 
sensitivity in NMR studies in vivo. The subject of overall imaging 
sensitivity in terms of the signal-to-noise ratio is further discussed in 
Chapter 5. 

NMR Spectra 

The observed NMR signal can be conveniently analyzed via 
Fourier transformation of the signal. In practice the so-called phase- 
sensitive detection technique is used to shift signal down in frequency 
by ujq prior to analog-to-digital conversion (see Appendix). Assuming 
that signal acquisition starts at time t = 0, the NMR spectrum in the 
case of a single resonance can be written as 

S{w)= f K e~ t/T2 +;V ^ dt, (1.4.3) 



where if is a real constant. Under the condition e J,p = 1, the real 
S re aiM and S imag (o;) components of the spectrum are given by: 

*- (w) = TTifc>' S '- (W) = TT#- (1A4) 

The shape of S rea ](w) in (1.4.4) is known as the Lorentzian lineshape. The 
linewidth of the spectrum denned as the width of S rea i at half height is 
2/T 2 . In general, spectral linewidths as well as lineshapes also depend 
on other factors such as magnetic field heterogeneity and chemical 
exchange. 

One of the main reasons responsible for the observed variety of 
NMR spectra is the chemical shift effect (CSE). CSE refers to the 
difference in Larmor frequencies of identical nuclei in different 
chemical environments. The basic phenomenon giving rise to CSE is 
the presence of microscopic currents induced by an external field B . 
The induced currents in atoms and molecules shield nuclei from B by 
creating an additional magnetic field that is proportional to B . As a 
result, the effective field acting on the nuclei can be expressed as 

B = Bo(l - a), (1.4.5) 
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Figure 1.4. H NMR spectrum obtained from a sample containing fat 
and water. 



where a is a constant (sometimes called the shielding or screening 
constant). It is a common practice to specify chemical shifts by using 
the dimensionless unit of parts per million (ppm). The chemical shift 
in ppm is defined as 

6 = Ui ~ Ur x 10 6 , (1.4.6) 

where w,- is the resonance frequency of the given nuclei and uj r is an 
arbitrary chosen reference frequency. By using Eq. (1.4.5) and taking 
into account that normally a < 1, we obtain 

6={a T - <Ti) x 10 6 , (1.4.7) 

where a t and a r are the shielding constants of the given nuclei and 
of the reference nuclei, respectively. For example, due to stronger 
shielding 1 H nuclei in the fat tissue have a slightly lower Larmor 
frequency than that of *H nuclei in water. The observed fat/water 
chemical shift is approximately 3.5 ppm (Figure 1.4). 

The chemical shift effect is important because it "fingerprints" 
chemically different populations of nuclei in a sample. That is, the pre- 
sence of chemically shifted peaks in NMR spectra makes it possible to 
determine the chemical composition of a sample as well as to measure 
the relative amounts of nuclei present in different chemical environ- 
ments. 

1.5. T 2 * DECAY 

Based on the Bloch equations, we could expect decay of the NMR 
signal to occur with a characteristic time T 2 . However, in practice a 
more rapid FID, known as the T\ decay, is normally observed. To 
explain this faster decay we first notice that the local magnetic field 
in media is always heterogeneous because of the nonuniformity of an 
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Immediately after the pulse 

t = to t = t,>t t = t 2 >t 1 

Figure 1.5. Dephasing of nuclear spins (indicated by arrows) in the 
rotating coordinate system after a 90 degree pulse. 

applied field or spatially varying magnetic susceptibility in media. As a 
result, the Larmor frequency, u , and phase of precession, <p = w t, 
become spatially dependent. Second, we notice that in all NMR experi- 
ments the measured signal is actually a sum of signals produced by 
numerous nuclei in the excited volume. Since the spread of phases in 
the volume increases with time (Figure 1.5), the destructive interfer- 
ence of the signals produced by different nuclei causes faster decay 
of the net signal as compared with the T 2 decay. 

To take into account the effect of dephasing due to macroscopic 
magnetic field inhomogeneities as well as T 2 decay caused by micro- 
scopic interactions between nuclei, the resultant decay of the transverse 
magnetization is frequently approximated by a phenomenological 
equation (compare with Eq. (1.3.6)): 

M xy = M xy (0) exp(-./u; i - t/TZ). (1.5.1) 

The rate of T 2 decay can be written as 

\/T 2 = 1/T 2 + 1/T 2 . (1.5.2) 

In this equation T 2 is a time constant of the decay that occurs due to 
the presence of magnetic field inhomogeneity, 8B. It should be noted 
that the exponential form of T 2 decay can only be considered as an 
arbitrary assumption which makes it possible to conveniently describe 
the effect of magnetic field nonuniformity in certain cases. 1 



1 An example of a nonexponential decay caused by magnetic field nonuniformity i 
presented in Chapter 6. 
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Magnetic field inhomogeneity, like chemical shift, is often expressed 
in terms of parts per million. The value of SB in ppm is defined as 

6B = ^x 10 6 , (1.5.3) 

where B is an arbitrary chosen reference field. A rather crude estimate 
for T 2 is given by (■y8B)~ 1 . In the case when T 2 is short compared to T 2 , 
the resultant signal decay is defined primarily by dephasing of spins. 
For example, T 2 « 25 msec for hydrogen nuclei at 1.5 Tesla and with 
local field inhomogeneities of 0.1 ppm. This T 2 value is significantly 
smaller than typical T 2 in biological tissues (40-150 msec) or T 2 for 
the 1 H nuclei in water (about 3 sec). 

1.6. SPIN ECHOES 

Fast dephasing of nuclear spins due to magnetic field hetero- 
geneity causes a significant loss of signal and thereby presents a 
serious problem for NMR imaging and spectroscopy. With regard to 
NMR spectroscopy, another detrimental effect of dephasing is line 
broadening which impedes identification of closely spaced spectral 
lines. In practice magnetic field homogeneity is improved through 
the use of an auxiliary magnetic field that can be adjusted to reduce 
spatial variations in the resultant magnetic field in a sample (this 
approach known as shimming is further discussed in Chapter 9). 
However, it turns out that even with shimming residual magnetic field 
inhomogeneities remain. It is therefore extremely important that 
dephasing of spins can be significantly reduced by using the remarkable 
phenomenon discovered by E. Hahn and known as spin echo [8]. 

To describe this phenomenon, we consider a simple model system, 
composed of two noninteracting magnetic moments \i x and \i 2 , in a 
static nonuniform magnetic field, B = kB. We can express B as a 
sum of a uniform field B and spatially varying field B' 

B = k(5 + B')- (1.6.1) 

In the absence of r.f. irradiation, the dynamics of the magnetic 
moments in the rotating reference frame {uj = ^B ) is governed by the 
effective magnetic field B - u/j = B' (see Section 1.2). Suppose that 
the magnetic moments, initially aligned in the z direction, experience 
a 90 degree pulse such that immediately after the pulse \i x and \i 2 lie 
along the /-axis (Figure 1.6(a)) in the rotating reference frame. 
Suppose also that a 180 degree pulse is subsequently applied along 
the x'-axis at time r after the excitation. Figure 1.6(b) depicts ji x and 
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(c) x r (d) 

Figure 1.6. Refocusing of spins in the rotating coordinate system. 



fi 2 immediately before the 180 degree pulse assuming that fij and fi 2 
have different Larmor frequencies wi < 7B0 and w 2 > 7-Bo. respectively. 
The 180 degree pulse inverts the /-components of jij and ji 2 (Figure 
1.6(c)). By examining the diagram in Figure 1.6(c) we can see that at 
time 2r (referred to as the echo time) the magnetic moments would be 
aligned again and would point in the negative /-direction in the rotat- 
ing reference frame (Figure 1.6(d)). Based on the fact that this result is 
valid for arbitrary moments \i x and ]i 2 , we can predict the phenomenon 
of refocusing of nuclear spins in a macroscopic sample that contains a 
great number of nuclei. The refocusing of spins manifests itself as a 
"spin-echo" signal (Figure 1.7). 

To take into account the effect of internuclear interactions, we 
can consider the Larmor frequencies of spins as the sums of two com- 
ponents u. < m acro an d <^micro> where oj macro is defined by a macroscopic 
(static) magnetic field in media and u> micro is defined by microscopic 
(time-dependent) magnetic fields produced by individual nuclei or by 
atomic electrons. Refocusing of spins following a 180 degree pulse 
requires that the Larmor frequencies of spins remain constant before 
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Figure 1.7. FID and spin-echo signals produced by a sequence of 90°- 
180° r.f. pulses. 

and after the pulse. In reality these frequencies constantly change due 
to the presence of microscopic magnetic fields. The 180 degree pulse 
does not change the irreversible phase dispersion caused by temporal 
variations in uj micro , and, therefore, does not alter T 2 decay of magneti- 
zation. In contrast, the phase dispersion caused by spatial variations in 
^macro is made zero at echo time It. Consequently, spin echoes can be 
used for accurate measurement of T 2 relaxation time. For example, 
in the pulse sequence shown in Figure 1.7 the magnitude of the 
observed echo signal, S, can be written as 

S(2t) = S(0) exp(-2r/T 2 ), (1.6.2) 

where S(0) is the signal magnitude immediately after excitation. This 
equation can be used to calculate T 2 from the measured spin-echo 
signals acquired with different echo times 2r. 

It is a common practice to use a series of spin echoes created by a 
train of 180 degree pulses applied after a single excitation pulse. This 
increases the accuracy of T 2 measurements without lengthening the 
total measurement time. One of the most frequently used pulse 
sequences for T 2 measurements, known as the Carr-Purcell-Meiboom- 
Gill pulse sequence (CPMG), can be written using the notation 
{90*- - r - 180y - 2r - 180y -It- 180 y . . .}. Notice that in the 
CPMG sequence the phase of the refocusing 180 degree pulses is shifted 
by 7r/2 with respect to the phase of the excitation pulse [9]. It can be 



Chapter 1 Basic Principles of Nuclear Magnetic Resonance 



shown that the CPMG sequence does not lead to accumulation of 
errors due to imperfections of 180 degree pulses in contrast to the 
Carr-Purcell sequence {90°- - r - 180°- - 2t - 180°. - 2t - 180° - . . .} 
[10], in which all r.f. pulses have the same phase. 

It should be noted that Hahn has studied spin echoes produced 
by repetitive 90 degree pulses [8]. Hahn has found that in a sequence 
of two identical r.f. pulses separated by a time interval r, a spin echo 
also occurs at time 2t as in the case described earlier. Three identical 
r.f. pulses generate a more complicated pattern of spin echoes. Hahn 
has discovered that in the case when a third pulse is applied at time 
T after the first pulse, additional echoes occur at T + r, 2T — 2r, 
2T -t, 2T. Moreover, Hahn has also found that the amplitudes of 
the echoes generated by the third pulse, known as the stimulated 
echoes, in general depend on 7\ and/ or T 2 . 

1.7. SIGNAL ATTENUATION DUE TO DIFFUSION 

In his classic work devoted to spin echoes Hahn has shown that 
diffusion in a nonuniform magnetic field contributes to decay of the 
transverse magnetization in media [8]. To analyze the effect of diffu- 
sion we start with the equation that describes the dynamics of magne- 
tization due to self-diffusion of particles in isotropic homogeneous 
media [6]: 

dM/dt = DV 2 M, (1.7.1.) 

where D is the diffusion coefficient. Torrey [11] has demonstrated that 
in order to incorporate the effects of 7\ and T 2 relaxations and preces- 
sion in an external magnetic field, the corresponding terms from the 
Bloch equations can be added to the right part of Eq. (1.7.1). Below 
we will obtain an exact solution of the Bloch equations in the case 
when diffusion of spins occurs in the presence of linear magnetic 
field gradients. 

In the presence of a non-uniform magnetic field B given by a sum of 
constant and linear terms: 

B = k(B + Gr), (1.7.2) 

the dynamics of magnetization in the rotating reference frame is 
described by the modified Bloch equations: 

dMJdt = -M x /T 2 + TGrAfy + DV 2 M X: (1.7.3a) 

dMy/dt = -M y /T 2 - 7 GrM ;c + DW 2 M y . (1.7.3b) 



where f(t) is an arbitrary function of time. By substituting the last 
equation into Eq. (1.7.4) and taking into account that /(0) = 1 we 
obtain 



(1.7.6) 
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Using M xy = M x + jM y again we obtain from the preceding two 
equations 

dM xy /dt = -M xy /T 2 -jfGrMxy + DV 2 M xy . (1.7.4) 

We seek a solution of this equation in the following form: 

M xy = M xy (0)f(t) exp I -t/T 2 - jyr j GdA, (1.7.5) 

n arbitrary function of time. By sub 
Eq. (1.7.4) and taking into account 

fit) = exp -D7 2 ( j G dt" dt' 

In a particular case when diffusion occurs in the presence of a time- 
independent gradient G = kG 2 we have f(t) — exp(-I>7 2 G^ 3 /3). In 
this case M xy is given by 

M xy = M xy (0) exp(-t/T 2 -hG z zt - D 7 2 G 2 i 3 /3). (1.7.7) 

Equation (1.7.7) defines the attenuation of the transverse magneti- 
zation as a function of diffusion coefficient and gradient strength. In 
principle, this equation can be used to calculate the diffusion co- 
efficient from a series of FIDs observed at different gradient strengths. 
However, in practice it can be difficult to differentiate between the 
effect of diffusion and signal attenuation due to dephasing of spins in 
the presence of an applied gradient. A better approach for NMR 
diffusion measurements was developed by Stejskal and Tanner [12]. 
The diffusion related attenuation of the NMR signal in the Stejskal- 
Tanner approach is achieved by applying two strong gradients 
symmetrically with respect to a 180 degree pulse (Figure 1.8). Because 
of the presence of a 180 degree pulse, the transverse magnetization of 
static spins is refocused at time TE after excitation. Therefore, the 
amplitude of the spin-echo signal from the static material is modulated 
by T 2 decay only. However, diffusion of spins in the direction of the 
applied gradients causes irreversible phase dispersion, which leads 
to additional signal attenuation. Under the assumption that the 
applied gradients are much greater than any intrinsic magnetic 
field gradients present in the sample, the effective attenuation of the 
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Figure 1.8. Stejskal-Tanner pulse sequence for diffusion 
ments. 

spin-echo signal is given by the well-known Stejskal-Tanner formula: 
S(TE) = S(0) exv[-TE/T 2 - Dj 2 G 2 6 2 {A - 6/3)], (1.7.8) 

where S(0) is the signal magnitude immediately after the excitation 
pulse; G is the amplitude of the applied gradients; 6 and A are the 
duration and time interval between the pulsed gradients, respectively 
(Figure 1.8). Measurements of the diffusion coefficient in the Stejkal- 
Tanner method can be carried out by varying the area under the 
gradient pulse, G6, or by varying A. 
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CHAPTER 2 



Excitation of the Transverse 
Magnetization 



In NMR imaging and spectroscopy it is a common practice to 
excite repeatedly a macroscopic system of nuclei by a series of r.f. 
pulses. For example, in the conventional approach used for spatial 
encoding in MRI (see Chapter 3) a number of excitations (e.g., 128 or 
256) are needed to collect all spatially encoded signals required for sub- 
sequent image reconstruction. One of the main assumptions of this 
approach is the existence of a steady-state, under which repetitive 
excitations would produce the same transverse magnetization in a 
sample. Using the Bloch equations we demonstrate that such a steady 
state can indeed be established as a result of the evolution of the 
nuclear magnetization subjected to a series of r.f. pulses. The last sec- 
tion of this chapter describes an example of spatially selective 
excitation that creates nonzero transverse magnetization only in a 
chosen slice of material. 

2.1. DYNAMICS OF REPEATEDLY EXCITED 
MAGNETIZATION 

In this section we consider the dynamics of nuclear magnetization 
in an external magnetic field B in the presence of a train of identical 
r.f. pulses (Figure 2.1(a)). To simplify further derivations we assume 
that the transverse magnetization in a sample is negligibly small 
immediately before the beginning of each successive excitation. This 
assumption is justified when the sequence repetition time, TR, is 
much longer than T 2 or when external magnetic field gradients 
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Figure 2.1. Different pulse trains used for repetitive excitation of the 
transverse magnetization. 



(sometimes called crusher or spoiling gradients) are applied after signal 
acquisition in order to disperse the transverse magnetization in a 
sample. A general case, when the steady-state transverse magnetization 
is nonzero immediately before the beginning of an r.f. pulse, is consid- 
ered in a classic paper by Ernst and Anderson [1]. 

During an excitation the effective magnetic field in the rotating 
reference frame is given by (see Chapter 1): 



5 eff 



= [K - io L y/Y + s, f . 



2il/2 



(2.1.1) 
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where ui is the frequency of the r.f. field and lj l is the local Larmor 
frequency. Let uj l = *yB + 8uj, where B is uniform and 8u> is defined by 
magnetic field nonuniformity or chemical shift effect. In the following 
derivations we assume that the irradiating r.f field is applied at the 
resonant frequency (i.e., u = jB ). We also assume that the amplitude 
of the r.f. field is large enough such that 6w < -yB T{ . Under the latter 
condition the effective magnetic field in Eq. (2.1.1) is equal to B l{ and 
is independent of the local Larmor frequency of nuclei. 

Using the Bloch equations (see Chapter 1) and assuming that the 
duration of the excitation pulses is shorter than 7\ and T 2 , we obtain 
the following equations for the longitudinal magnetization M z : 

MM = M , 
M Z {TR) = M [(l - e~ TR ^) + e- TR ^ cos a], 
M Z (2TR) = M [(l - e' TR ^)(l + e~ TR ^ cos a) + e' 2 ™^ cos 2 a], 
M z (nTR) = M [(1 - e -TRIT^ e -kmi Tl cos * Q + e -nTR /Tl ^ 1 _ 

In the preceding equations the 2-axis is chosen in the direction of B ; 
M z (nTR) is the longitudinal magnetization immediately before the 
(n + 1) pulse; a is the angle of rotation of the magnetization (known 
as the flip angle) in the rotating reference frame (see Chapter 1); M 
is the equilibrium magnetization. Using the formula for geometric 
progression the expression for M z (nTR) can be written as 

M z (nTR) = Mo [(1 - e~ TR ^) l ~ 6 ""^' COB " Q! + e -™/* cos" a] - 
L l-e 1K/l1 cosa J 

(2.1.2) 

This equation describes the longitudinal magnetization as a function 
of T x relaxation time, repetition time, flip angle, and number of 
excitations. As the number of excitations increases, the longitudinal 
magnetization reaches its steady-state 

where E 1 = e~ TR ^ Tl . Note that in the case when a = tt/2, the steady- 
state is established after the first pulse regardless of the ratio TR/T X . 
In a general case, however, the relaxation to the steady-state in 
(2.1.3) occurs much more slowly and the relaxation rate depends 
upon a as well as TR/T X (Figure 2.2). 
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Figure 2.2. Dependence of the longitudinal magnetization M z (nTR) o 
the number of excitations at 7\ of 800 ms and 77? of 80 ms. 



Consider the equation for the steady-state transverse magnetization 
at an arbitrary time TE after an excitation. In the case when dephasing 
of spins can be neglected (i.e., when SuiTE -C 1), the steady-state trans- 
verse magnetization in the rotating reference frame can be expressed 
as M tr (TE) = M Z {TR) e~ TE/T2 sin a. Using (2.1.3) we obtain 



M tr {TE) = M E 2 sin a - 



(2.1.4) 



where E 2 = e~ TEITl . Using this equation it is easy to verify that the 
steady-state transverse magnetization decreases with decreasing TR. 
This conclusion also follows from the observation that short TR does 
not allow enough time for relaxation of the longitudinal magnetization 
in a sample after an excitation pulse. In this case the sample is 
frequently referred to as saturated. Equation (2.1.4) shows that M tr 
increases with decreasing TE because of decreased T 2 decay. These 
simple results turn out to be extremely important because they can 
be used to enhance contrast in MR images (see Chapters 4 and 8). 
The dependence of M tr on flip angle is discussed in section 2.2. 



Spin-Echo Sequence 

To avoid signal loss due to magnetic field nonuniformity, 180 degree 
pulses are frequently included in pulse sequences used in NMR imaging 
and spectroscopy. A spin-echo pulse sequence is shown in Figure 2.1(b). 
The dynamics of magnetization in this case can be studied by utilizing 
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the same approach we used to derive Eq. (2.1.3). However, to avoid cum- 
bersome derivations we will consider at this time only the equations 
describing the steady-state magnetization. 

Let M z _ ( TR + TE/2) and M z+ (TR + TE/2) denote the steady-state 
longitudinal magnetization immediately before and after a 180 degree 
pulse, respectively. From the Bloch equations it follows that 

M Z _(TR + TE/2) = \M Z {TR) cos a- M ] e - TE/2Tl + M , 

M Z+ (TR + TE/2) = -M Z _(TR + TE/2), 

M Z (TR) = [M Z+ {TR + TE/2) - M ] e ^ TR - TE ^l T ^ + M , 

where M Z (TR) is the steady-state magnetization immediately before an 
excitation pulse. Using these equations we obtain 

-TR/Ti _ n-iTR-TEM/T! 



,-TR/Ti 



(2.1.5) 



1 cos a 

From (2.1.5) it follows that the steady-state transverse magnetization at 
time TE after excitation is given by 

1 . ,-ra/T, _ 2p -(TR-TE/2)/T 1 

M tr (TR) = M e- TE/T > sina— —- ^ . (2.1.6) 

l + e -lR/l 1 CQSa 

Inversion-Recovery Sequence 

A basic inversion-recovery pulse sequence is shown in Figure 
2.1(c). During each cycle the first 180 degree pulse (referred to below 
as the inversion pulse) inverts the longitudinal magnetization, which 
then recovers during the time interval TI, known as the inversion 
time. The subsequent excitation and refocusing of the magnetization 
are essentially the same as in a spin-echo sequence (Figure 2.1(b)). It 
is clear that M z is zero immediately after a 90 degree pulse. Following 
the pulse the longitudinal magnetization recovers during time TE/2. 
Taking into account that M z is subsequently inverted again by a 180 
degree pulse applied at time TE/2 after the 90 degree pulse we obtain: 

M Z (TR) = M [l - 2e-< TR - TI - TE W T > + e ^ TR - T ^% (2.1.7) 

where M Z (TR) denotes the steady-state longitudinal magnetization 
immediately before an inversion pulse. 

The steady-state transverse magnetization at time TE after a 90 
degree pulse (Figure 2.1(c)) can be expressed as 

M tr (TE) = [M - (M Z (TR) + M ) e" 77 ^] e~ TE ^ 

= M e~ TE/T * [1 - 2e- n ^ + 2 e^ m ' TEI2) ^ - e"™/^]. (2.1.8) 
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In the case when TE < 7\ we can simplify this equation as follows: 

M tr {TE) = M e- TE ^[1 - 2e~ TI/T i + e-™^}. (2.1.9) 

According to Eq. (2.1.9), it is possible to null the signal from specific 
nuclei with known T x value by choosing the inversion time to be 

r7=Tlln [ l + e 2 ^ ]' (2 - L10) 

This feature of inversion-recovery sequence is often used in diagnostic 
MRI in order to eliminate the signal from a chosen tissue such as fat or 
cerebrospinal fluid (CSF). For example when a lesion in the brain par- 
enchyma is bright on a T 2 -weighted image (see Chapter 4) but obscured 
by the even brighter CSF in the adjacent ventricle, suppression of the 
CSF signal through the use of inversion-recovery sequence can 
improve the contrast between the brain tumor and the surrounding 
tissue. 

2.2. ERNST ANGLE 

Steady-state transverse magnetization established in the presence 
of r.f. pulses depends on a number of parameters, which can be divided 
into two different groups. The first group is composed of parameters 
such as M , Tj, and T 2 relaxation times that are defined by the proper- 
ties of the excited nuclei, their temperature and surroundings, and the 
strength of the applied magnetic field. The second group consists of 
user-controlled parameters, such as TR, TE, and flip angle, which 
can be chosen to maximize the observed signal or enhance image con- 
trast (see Chapter 4). It is easy to verify that the maximum transverse 
magnetization in (2.1.4) and (2.1.6) is achieved when TR > 7\ and 
a = 7r/2. However, in order to ensure desirable image contrast or 
shorten scan time (see Chapters 4 and 7) it is often necessary to use 
TR shorter than 7\. In the latter case Eqs. (2.1.4) or (2.1.6) can be 
used to determine the optimum flip angle, known as the Ernst angle, 
which corresponds to the maximum M tT (Figure 2.3). 

Let us first consider Eq. (2.1.4) which describes the steady-state 
transverse magnetization in the pulse sequence in Figure 2.1(a). 
Using Eq. (2.1.4) it is easy to show that the Ernst angle, a E , satisfies 
the following equation [1] 

cosa £ -exp(-Ti?/T 1 ). (2.2.1) 

From Eq. (2.2.1) it follows that q £ increases with TR. In the case 
when TR » T x we have a E ss tt/2 and M tr {TE) « M exp(-TE/T 2 ). 
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Conversely, in the case when TR <§: T x the Ernst angle is approximately 
equal to {2TR/T i ) 1/2 and M tr (TE) * M (TR/2T 1 ) 1/2 exp(-TE/T 2 ). 

The choice of flip angle becomes particularly important in the case 
of rapid gradient-echo MR imaging (see Chapter 7) which is performed 
with short TR to decrease scan time. Note that when imaging with 
TR <: 7\ and large flip angle (i.e. a « 90 degrees) the signal is propor- 
tional to TR/T l . Conversely, when imaging with the Ernst flip angle, 
the signal varies as (TR/T-^) 1 ^ 2 . As a result, to avoid a significant 
loss of signal, gradient-echo imaging with short TR is typically imple- 
mented with low flip angle that does not significantly deviate from the 
Ernst angle. 

From Eq. (2.1.6) it follows that the Ernst angle for a spin-echo 
sequence is defined by the equation 

cos a E = - exp(-TR/T{). (2.2.2) 

Note that, in contrast to the previous case, the Ernst angle for a spin- 
echo sequence is greater than n/2. 

2.3. SPATIALLY SELECTIVE EXCITATION 

Magnetic resonance imaging is typically used to image an object 
in three dimensions. To achieve spatial localization it is a common 
practice to first selectively excite transverse magnetization in a thin 
slice of material, which is subsequently imaged in the two remaining 
directions. The design of spatially selective pulses is an interesting 
and important subject that has been studied by many investigators. 
To elucidate the principles of spatially selective excitation, we will 
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consider the dynamics of the transverse magnetization in the presence 
of a sine shaped r.f. pulse. The approach used in this section is based on 
the small-flip-angle approximation that makes it possible to solve the 
Bloch equations analytically [2-4]. 

Suppose that a short r.f. pulse with duration T p < T 2 , T x is applied 
in the presence of an external magnetic field with linear gradient: 
B = k(B + Gz), where k is a unit vector in the z direction. During 
the excitation the dynamics of the nuclear magnetization in the 
rotating reference frame is described by the following equations: 

dM * n *, 

at 

— / = --yGzM x +-yB l{ (t)M z , 
at 

assuming that the r.f. field has only ^-component, B r{ (t). 

The small-flip-angle approximation is based on the assumption 
that during excitation the magnetization remains close to the z-axis 
so that {M 2 X + M]f 12 < M z « M . Let 

M xy = M x +jM y . (2.3.1) 

Neglecting the difference between M z and M we obtain: 

^ = -hGzM xy +hB r , i .(t)M . (2.3.2) 

The solution of the last equation with initial condition M xy (0) = is 
given by 



M xy =nM / B r . f .(s)exp 



j-yz / G(u)du 



(2.3.3) 



This equation defines the transverse magnetization as a function of the 
time-varying r.f. field and applied gradient. In principle, Eq. (2.3.3) can 
be used to design r.f. pulses and gradient waveforms for excitation of 
arbitrary spatial profiles of the transverse magnetization in a sample. 
Several interesting examples of spatially selective excitation in one 
and two dimensions have been described by Pauly et al. [4]. 

An important practical example is the excitation by a smc-shaped 
r.f. pulse 

£ r .^=A Sin "f;; o) , (2.3.4) 
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where A and u b are constants. The spectrum of the pulse in (2.3.4) is 
given by 

B rl M = { COnSt - eM °> tf M«** (2.3.5) 

assuming that uj b > 0. A smc-shaped pulse with the spectrum denned in 
Eq. (2.3.5) will excite only the nuclei with Larmor frequencies centered 
near the irradiating frequency w = 7#o within the bandwidth of the 
pulse, 2u b } In the presence of a static, linear gradient, Larmor frequency 
of nuclei becomes spatially dependent: u L = u + ^Gz. Because of the 
spatially varying lo l , excitation of the transverse magnetization by a 
smc-shaped pulse occurs only inside a slice of material that contains 
nuclei with Larmor frequencies \co L - u \ < uj b . It then follows that the 
excited slice thickness is given by 

Sz = 2\u b /-yG\. (2.3.6) 

We can use Eq. (2.3.3) to determine the transverse magnetization 
generated by a sine pulse. Under the condition that G remains constant 
during the pulse, this equation can be written as 

T p /2 

M xy (T p )=j'yAM exv(- J -^^ J eMhGzu) 3 ^ du, 

-T p /2. 

(2.3.7) 

where T p = 2t is the duration of the pulse. In the case when T p » l/wj,, 
the integral in Eq. (2.3.7) can be approximated by Fourier transform 
of the sine function and the resulting expression for M xy is given 
by 

(jwyAMo jGT/2 . flr , 
M xy (T p ) = \-^ b — e '> lf l^l<^ (2.3.8) 

I 0, if | 7 Gz| > u; b 

Notice that according to the above equation the transverse mag- 
netization acquires a phase ip = -"/GT p z/2. Because of the accumu- 
lated phase, the signal from the excited slice will be very small if no 
measures to cancel this phase are taken. Fortunately, refocusing of 

1 A smc-shaped r.f. pulse applied in the absence of external field gradients provides an 
example of a frequency (spectrally) selective pulse that only excites a limited range of 
frequencies. In MRI frequency selective pulses are often used for chemical shift imaging 
(see Chapter 6). 
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Figure 2.4. Schematic of spatially selective excitation. After the excita- 
tion r.f. pulse, the polarity of the slice-select gradient is reversed in 
order to refocus the transverse magnetization in the excited slice. 

the magnetization can be easily achieved by reversing the polarity of 
the slice-select gradient after the excitation (Figure 2.4). It is easy to 
verify that the area of the refocusing gradient lobe must be half of 
the area of the gradient lobe used for slice selection in order to null 
the phase accumulated during excitation. Figure 2.5 shows an example 
of a 25 degree sine pulse applied along the x-axis in the rotating frame 
of reference and the corresponding components of the transverse 
magnetization calculated numerically. The deviations from the ideal 
excitation profile are mostly due to a small number of side lobes of 
the sine function used in this example. Although the slice profile can 
be improved by using a larger number of side lobes, the penalty is 
prolonged pulse duration which is undesirable for a number of applica- 
tions such as flow imaging, rapid gradient-echo imaging etc., when 
short excitation pulses (a few milliseconds) are needed to reduce 
dephasing of spins or shorten repetition time. 

The small-flip-angle approximation makes it possible to design r.f. 
pulses with flip angles on the order of 7r/2. Spatially selective pulses 
with flip angles greater than it/2 are frequently designed by using 
numerical solutions of the Bloch equations. Alternatively, r.f. pulses 
can be designed by using an approach independently suggested by 
Shinnar and Le Roux [5]. 

An important problem in the design of spatially selective pulses is 
the deviation from the intended excitation profile caused by the hetero- 
geneity of the r.f. field. Distortion of slice profile can be particularly 
severe in the case when r.f. pulses are generated by using surface r.f. 
coils which produce an extremely inhomogeneous B r f field. The effect 
of r.f. inhomogeneity can be reduced by using so called "adiabatic 
pulses" [6,7] which are designed to compensate for the spatial variations 
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Figure 2.5. Example of spatially selective excitation: (a) amplitude of a 
25 degree sine pulse applied along the x axis in the rotating reference 
frame; (b) transverse components of the magnetization in the rotating 
reference frame after the excitation and gradient refocusing. 

in B r{ and are characterized by simultaneous modulation of the r.f. 
amplitude and frequency. 
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CHAPTER 3 



Basic Techniques for 2D and 3D MRI 



3.1. IMAGE RECONSTRUCTION FROM DISCRETE 
SAMPLES 

Because image reconstruction in MRI typically involves the use of 
Fourier transform (FT), we need to begin our discussion by describing 
the basic principles of Fourier analysis. The continuous Fourier trans- 
form of a function f(x) is denned as 

s(k x ) = J~j( x y k * x dx, (3.1.1) 

where j = y/-[. If the Fourier transform S(k x ) is known, then the func- 
tion f(x) is given by inverse Fourier transform of S(k x ): 



«■»-££ 



S(k x )e- Jk ' x dk x . (3.1.2) 



For the purpose of convenience, we will refer to f(x) as the object. The 
process of measuring S(k x ) will be referred to as sampling, and the 
function S(k x ) itself will be referred to as the signal. 

In practice it is not possible to obtain S(k x ) as a function of 
continuous k x . Instead S{k x ) is measured only at a finite number of 
sampling points. To demonstrate how discrete samples of the signal 
can be used for image reconstruction, let us assume that S(k x ) in 
Eq. (3.1.1) is measured at N locations: 

k x = 2nn/L x . (3.1.3) 
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In this equation n is an integer, -N/2 <n< N/2, and N is assumed to be 
even; L x is a positive constant, known as field-of-view, which defines the 
size of the imaged region. Using Equations (3.1.1) and (3.1.3) we obtain 

S(n) = S(k x (n)) = f°° f(x)e j2mx/L * dx. (3.1.4) 

We now define the reconstructed image intensity, I, as the inverse 
discrete FT of the signal in (3.1.4) [1]: 

ra=2V/2-l 

I(p) = {\/N) £ S(n)eM-j2nnp/N). (3-1-5) 

n=-N/2 

In this equation p is an integer, -N/2 <p< N/2. 1 After substituting 
S(k x ) from Eq. (3.1.4) into Eq. (3.1.5) and using the formula for geo- 
metric progression, we obtain the following equation describing the 
relationship between the object and its image: 

I(p) = r f( x )PSF(x - x p ) dx, (3.1.6) 

where x p = pL x /N and 

v ' Nsm(irx/L x ) v ' 

From Eq. (3.1.6) it follows that the image intensity is given by the 
convolution of the object f{x) and PSF(x), known as point-spread 
function. The term point-spread function is used because PSF(x) repre- 
sents the image of a point source described by the Dirac delta function, 
S(x) 2 

These results can be generalized to describe Fourier reconstruc- 
tion from discrete samples in two or three dimensions. For example, 
a two-dimensional image intensity is given by the convolution of 



1 Discrete FT of the signal can be computed very efficiently by using the approach 
described by Cooley and Tukey in 1965 and known as the fast Fourier transform (FFT). 
The majority of FFT algorithms require that the number of data points be an integer 
power of two. As a result, the number of signal samples in applications which entail 
the use of FFT is most likely to be a power of two. 

2 The Dirac delta function has the properties: 

if x ^ then 6(x) = 
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the object and a two-dimensional point-spread function PSF(x,y), 
which is given by the product of one-dimensional point-spread func- 
tions PSF x (x) and PSF y (y). 

Spatial resolution 

Image intensity in Eq. (3.1.6) can be expressed as a sum of weighted 
elements f(x) dx with PSF(x) as the weighting function [2]. It is there- 
fore clear that the image intensity has contributions from different 
locations in the object. The ability of an image to accurately reproduce 
spatial variations in the object is characterized by the spatial resolu- 
tion of the image. Spatial resolution is often defined as the minimum 
distance between two points in an object at which they can still be 
distinguished from one another in the image. Based on Eq. (3.1.6) 
spatial resolution, 6x, can also be defined as the effective width of a 
point-spread function [3]: 



By using the latter definition and Eq. (3.1.7) it can be shown that 

6x = LJN. 

Aliasing 

A one-dimensional point-spread from Eq. (3.1.7) is plotted in 
Figure 3.1. Notice that the PSF(x) oscillates such that its amplitude 
reaches maximum at the origin x = and rapidly decreases with 
increasing distance from the origin. Since PSF(x) is also a periodical 
function with period L x , the reconstructed image intensity I(p) is 
defined primarily by a sum of the equally weighted object intensities at 
x = —r^- + mL x , where to = 0, ±1, ±2 — This phenomenon (known as 
aliasing) results from discrete sampling of S(k) and in general prevents 
complete recovery of f(x). However, under the circumstances typical 
for MR imaging the object has a limited spatial extent such that 

(\f(x)\ >0, if -a/2<x<a/2 

\\f{x)\ = 0, otherwise. [ ' ' 

According to the sampling theorem [1], in this case aliasing can be 
avoided and f{x) in principle can be fully recovered by satisfying the 
Nyquist criterion for the sampling interval Ak: 

Ak<—. (3.1.10) 
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Figure 3.1. One-dimensional point-spread function with iV = 32. The 
solid and dashed lines represent the real and imaginary parts of the 
PSF, respectively. 

By examining this equation and Eq. (3.1.3) we conclude that in order to 
prevent aliasing the selected field-of-view L x should be equal to or 
greater than the length of the object, a. 

3.2. FREQUENCY AND PHASE ENCODING 

In practice MR imaging is most frequently performed by using 
the approach suggested by Kumar et al. [4]. This approach employs a 
sequence of pulsed magnetic field gradients which are applied during 
a free induction decay or a spin echo to ensure that the signal is 
given by Fourier transform of the transverse magnetization in the 
sample. At this stage we are ready to discuss two basic imaging 
techniques for Fourier encoding: frequency encoding and phase 
encoding. 

Frequency Encoding 

Frequency encoding is implemented by acquiring signal in the 
presence of an external magnetic field gradient. The purpose of the gradi- 
ent, known as the frequency-encoding or readout gradient, is to make the 
Larmor frequency of nuclei spatially dependent during signal acquisition. 

The signal acquired in the presence of a readout gradient is 
composed of components with frequencies from a narrow range, 
known as the signal bandwidth, around the Larmor frequency w at 
the center of the field-of-view. The frequency w is also known as the 
reference frequency. Each of the signal components is produced by 
spins from a certain location in the object. Because the typical signal 
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bandwidth is much smaller than w , 3 it is important to remove the refer- 
ence frequency from the signal in order to distinguish between the 
spectral components within the bandwidth. In order to do this a special 
technique, referred to as the phase-sensitive detection (see Appendix), is 
used to shift the signal down in frequency by i% As a result the output 
signal consists of primarily low frequency components. 4 

At this stage we consider signal acquisition in the presence of a 
readout gradient, G x = dB z /dx, given by 

f -G ,*, to < t < t + TJ2 
G x = I G 0iX , t + TJ2 <t<t + 3T s /2 (3.2.1) 

[ 0, elsewhere 

where Go,* is a constant (Figure 3.2). Suppose that the signal acquisi- 
tion starts at t + TJ2 and ends at t + 3TJ2. The phase of the trans- 
verse magnetization accumulated during the interval [t , t + T s /2] in 
the presence of the first (defocusing) lobe of the readout gradient can 
be written as 

4> (x) = 1 xJ t ° + S G x dt = - 7 G ^T s x/2. (3.2.2) 

The readout gradient reversal at t = t + T s /2 initially causes refocusing 
of spins, which in turn leads to a signal increase, known as the gradient 
echo, such that the maximum signal occurs in the center of the acquisition 
interval (Figure 3.2). It can easily be shown that during the acquisition 
the phase of the transverse magnetization is given by 

</>(x, t) = >yG 0jc (t - t )x - ~/G , x T s x. (3.2.3) 

Suppose that N samples of the signal are collected at t = t + T s /2 + 
pr w , where an integer p changes from to N - 1. The sampling interval 
T w — T s /N is known as the dwell time and T s is referred to as the 
acquisition or readout time. Using Eq. (3.2.3) we can express the 
phase of magnetization as 

<t>(x, t{n)) = -/G 0tX T w nx, (3.2.4) 

where n=p- N/2 and t(n) = t + T s + nr w . 



3 For example, the Larmor frequency of 'H nuclei at 1.5 Tesla is about 63.86 MHz, 
whereas typical signal bandwidth in MRI is less than 100 Khz. 

4 Strictly speaking, during phase-sensitive detection the signal is shifted in frequency by 
+/ - lo . Therefore, the output signal also includes high-frequency components that are 
centered at 2w . However, these components are subsequently removed from the signal 
by passing them through a low-pass filter. 
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Figure 3.2. Formation of a gradient-echo: (a) readout gradient wave- 
form; (b) accumulated phase; (c) NMR signal. 



To understand how NMR signal obtained with frequency encoding 
can be used for image reconstruction we need to consider the equation 
describing the relationship between the signal and magnetization in 
the object. We will assume for simplicity that T s <C T 2 and that the 
effect of magnetic field inhomogeneities during data acquisition can 
be neglected. We will also assume that the transmitter and receiver 
coils used for excitation and signal detection, respectively, generate 
uniform r.f. fields. From the principle of reciprocity (see Chapter 1), 
it then follows that the sampled signal in a one-dimensional case can 
be written as 

S(n) = S(t{n)) = £ f M xy {x)e^ {xAn)) dx = £ f M xy (x)e jk * {n)x dx. (3.2.5) 

In this equation £ is a constant that, for simplicity, will be given the 
value of one in the derivations below, M xy is the (complex) transverse 
magnetization in the object (see Appendix), and 



k x {n) = 7G 0iX t ii 



(3.2.6) 



Frequency and Phase Encoding 



From Eq. (3.2.5) it follows that the acquired signal is defined by Fourier 
transform of the magnetization in the object. Therefore, according to 
the results obtained in the previous section, inverse discrete FT of 
the signal in Eq. (3.2.5) will reconstruct M xy with spatial resolution 
L x /N, where 

L x = 27r/Aft x = 2ir/{-yG . x T w ). (3.2.7) 



Sampling 

To understand better the important relationship between the field- 
of-view, readout gradient, and dwell time in Eq. (3.2.7), we need to 
recall that according to the sampling theorem [1] a continuous func- 
tion f(i) can be reconstructed fully from a series of discrete samples 
provided that: a) FT{f(t)} is zero for all frequencies \v\ > f max ; and b) 
sampling interval does not exceed l/2u max . Conversely, if sampling 
interval exceeds l/2z/ m(M then some of the high-frequency components 
of f{t) cannot be determined. The insufficiently high sampling rate in 
the latter case, often referred to as under sampling, causes aliasing, 
which generally makes it impossible to reconstruct a continuous 
function faithfully from its discrete samples. 

The signal bandwidth in MR imaging is defined by the readout 
gradient and field-of-view 

2l/max= l^h. ( 3. 2 . 8) 

The signal is digitized at a rate of r^ 1 by using a special device called 
analog-to-digital converter (ADC). By examining Eq. (3.2.7) in terms of 
the sampling theorem, we realize that this equation ensures that the 
sampling rate is sufficient (i.e., t w = l/2i/ max ) to avoid aliasing when 
reconstructing a continuous NMR signal from its discrete samples. 

Phase Encoding 

Two-dimensional spatial encoding in MRI is normally achieved 
through the use of an additional gradient, known as the phase-encoding 
gradient, which is perpendicular to the frequency-encoding gradient. 
In the presence of the phase-encoding gradient, G y = dB z /dy, applied 
prior to readout, the transverse magnetization acquires a phase 



yf " h G y dt, 



4>(y)=7yJ t G y dt, (3.2.9) 

where t ph is the duration of the gradient. As a result of frequency and 
phase encodings the NMR signal in a two-dimensional case is given by 



■//' 
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2D Fourier transform of the magnetization: 

S= If M xy (x,y)e jk * in)x+j * M dxdy 

f M xy (x,y)e jk " {n)x+jk > y dxdy, (3.2.10) 

where k y = 1 f t +tph G y dt. 

Phase encoding requires repetitive excitations of the transverse 
magnetization in the object in order to collect signals at different 
values of k y . In the initially proposed phase-encoding scheme, changes 
in k y were achieved by varying the duration of the phase-encoding 
gradient while keeping its strength constant [4]. A more common 
approach is to vary the strength of the phase-encoding gradient in a 
step-like fashion; that is, 

G y = G , y m, (3.2.11) 

while its duration, t p h, is kept constant [5]. In Eq. (3.2.11) G 0<y is a 
constant, and an integer m changes from -M/2 to M/2 - 1 (M is the 
total number of phase-encoding steps). 

Using Eq. (3.2.11), the signal obtained with frequency and phase 
encodings can be written as 

S(n,m)= J J M xy (x,y)e ik ' {n)x+jk ' im)y dxdy, (3.2.12) 

where k x (n) = jG 0x T w n and k y {m) = -yG 0y t ph m. After all phase encod- 
ings are implemented, image reconstruction is performed by computing 
inverse discrete FT of S(n,m). Using results from the previous section 
we find that the reconstructed image intensity is given by the con- 
volution of the transverse magnetization in the sample and a two- 
dimensional PSF. The resultant spatial resolution in the x-y plane is 
(L x /N) x (Ly/M), where L y is the field-of-view in the direction of the 
phase-encoding gradient: 

L y = 2n/Ak y = 2n/( 7 G J ph ). (3.2.13) 

Note that the reconstructed image intensity is a complex quantity 
which can be written in polar form Ae J , where A and 8 are the 
amplitude and phase of the intensity, respectively. It is the amplitude 
of the intensity that is normally used to display MR images (often 
referred to as magnitude images), while the phase information is 
neglected. However, in several instances (e.g., imaging of flow, imaging 
of static magnetic field inhomogeneities etc.) the reconstructed phase 
of the intensity may also be used. Basic aspects of image display are 
discussed in the Appendix. 
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3.3. BASIC GRADIENT AND SPIN-ECHO PULSE 
SEQUENCES 



Gradient-Echo S 

A gradient-echo pulse sequence begins with excitation of the 
transverse magnetization in a slice of material. After the excitation 
the polarity of the slice-select gradient (shown as G z in Figure 3.3) is 
reversed in order to refocus spins in the slice (see Chapter 2). The 
excited slice is subsequently imaged in the x- and ^-directions using 
frequency and phase encodings. The pulse sequence including spatially 
selective excitation, phase encoding and readout is repeated M times 
with different values of the phase-encoding gradient as described in 
the previous section. If TR denotes the time interval between consecu- 
tive excitations (repetition time), then the total scan time, T scan , equals 
MxTR. 

Reconstructed image intensity is defined by the sum of signals 
produced by spins from the volume 

SV = 6x x Sy x L tl , (3.3.1) 

where 6x = L x /N and 6y = L y /M, L s! is the slice thickness. This result 
follows from the fact that image intensity is given by the convolution of 
the transverse magnetization in a slice and a two-dimensional point- 
spread function that has widths Sx and 6y in the x- and j-directions, 
respectively. The volume SV is known as an imaging voxel and the 
area 6x x Sy is known as a pixel. 



-J 



"Y 



Figure 3.3. Basic gradient-echo sequence. Readout, phase-encoding 
and slice-select gradients are shown as G x , G y , and G z respectively. 
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Gradient-echo imaging often suffers from signal loss caused by 
dephasing of spins in the presence of magnetic field inhomogeneities. 
An empirical formula frequently used to evaluate gradient-echo 
image intensity is 

I<xe~ TE/T \ (3.3.2) 

where TE, known as the echo time, is the time interval between the center 
of the excitation r.f. pulse and the gradient echo (see Figure 3.3). Notice 
that T\ is used in Eq. (3.3.2) instead of T 2 in order to take into account 
the effect of dephasing of spins within a voxel. Since dephasing of spins 
increases with increasing TE, gradient-echo images are often acquired 
with very short echo time (a few milliseconds) in order to reduce signal 
loss. The effect of intravoxel dephasing is discussed further in Chapter 6. 
Gradient-echo imaging with short TR was originally suggested as 
a technique for rapid in vivo MRI in order to decrease scan time and 
reduce artifacts due to the motion of the object [6-8]. Since its 
inception rapid gradient-echo imaging has found a number of impor- 
tant clinical applications (e.g., abdominal imaging, MR angiography, 
and contrast-enhanced MRI). Rapid gradient-echo imaging is dis- 
cussed further in Chapter 7. 

Spin-Echo Sequence 

The main difference between gradient-echo and spin-echo pulse 
sequences is that the latter also includes a 180 degree pulse that 
causes formation of a spin echo during signal acquisition (Figure 
3.4). Initial dephasing of spins in the excited slice can be achieved by 
applying the defocusing lobe of the readout gradient before the 180 
degree pulse. As in gradient-echo imaging, the area of the defocusing 
lobe equals one-half of the area of the refocusing lobe. The position 
of the refocusing gradient lobe is typically chosen in such a way that 
the peak of the spin echo occurs in the center of the acquisition 
interval. Note that because the 180 degree pulse inverts the phase of 
the magnetization, the defocusing and refocusing lobes of the readout 
gradient have the same polarities. 

Because of spin refocusing following a 180 degree pulse, spin-echo 
imaging is generally less sensitive to the presence of magnetic field 
inhomogeneities than gradient-echo imaging. Spin-echo image inten- 
sity can be approximated as follows: 

Ioce- TE ' T \ (3.3.3) 

where TE is the time between the excitation pulse and the formed 
spin echo (Figure 3.4). Equation (3.3.3) shows that, in contrast to 
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Figure 3.4. Basic spin-echo sequence. 



t-echo imaging, spin-echo image intensity depends on the ratio 
of TE to T 2 rather than T\. Because local T 2 values are found to be 
tissue dependent, spin-echo imaging can provide excellent contrast 
between different tissues in diagnostic MRI (see Chapter 4). 

It should be noted that the use of 180 degree pulses in spin-echo 
imaging results in longer echo time and repetition time relative to 
those in gradient echo imaging. As a result spin-echo imaging is gener- 
ally slower than gradient-echo imaging. An approach that significantly 
shortens scan time in spin-echo imaging is discussed in Chapter 7. 

2D Multi-Slice Imaging 

In the case when repetition time TR exceeds the time required for 
spatially selective excitation, phase encoding, and readout, a number 
of slices can successively be imaged during one TR interval [9]. Figure 
3.5 shows examples of different excitation schemes for multi-slice 
scanning with two-dimensional frequency and phase encodings, often 
referred to as conventional or standard MR imaging. The maximum 
number of imaged slices is given by 

TR (3.3.4) 



Nmax = 



TE + 6T' 



where 6T defines the additional time (approximately half of the read- 
out time) needed to complete signal acquisition and switch gradients. 
For example, assuming TR of 800 msec, TE of 70 msec, and ST of 10 
msec, we obtain N max = 800/(70 + 10) = 10 slices. 
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Figure 3.5. Sequential (a) and interleaved (b) excitation schemes for 
multi-slice imaging. 

To prevent artifacts and avoid a reduction in the signal-to-noise 
ratio in multi-slice imaging, simultaneous excitation of nuclei from 
different slices should be avoided. This requirement is difficult to 
satisfy in practice because of imperfections of excitation pulses. 
Experience shows that signals from neighboring slices often interfere 
with each other. One solution to the problem of slice interference is 
to use a gap between imaged slices. Another practical solution is to 
alternate slices during excitation (Figure 3.5b). The slice interference 
can also be reduced by using long r.f. pulses, which in general produce 
better slice profiles than those created by short r.f. pulses. 

3D Imaging 

The minimum slice thickness in single slice or multi-slice 2D 
imaging is limited by the maximum gradient strength and duration 
of the spectrally selective r.f. pulse used for excitation of the transverse 
magnetization. Typical slice thickness in 2D imaging is about 1-5 mm. 
Thinner slices are usually acquired by using Fourier encoding in three 
dimensions (3D FT) [10]. The principal difference between 2D imaging 
and 3D imaging is that the latter includes an additional phase- 
encoding gradient that is applied in the direction of the slice-select 
gradient. A typical 3D FT gradient-echo sequence is shown in Figure 
3.6. The spatial resolution in the z direction is L s //Q, where L si is the 
slice thickness and Q is the number of phase-encoding steps in the z- 
direction. Compared with 2D multi-slice scanning, the 3D approach 
allows higher spatial resolution and higher signal-to-noise ratio (dis- 
cussed in Chapter 5). Three-dimensional imaging is normally imple- 
mented with relatively thick slabs (i.e., 1-10 cm) while using phase 
encoding in the slab direction to achieve high spatial resolution. As a 
result, 3D imaging can be performed with very short excitation 
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Figure 3.6. Basic gradient-echo sequence for 3D imaging. 

pulses, which makes it possible to significantly shorten TE and TR. 
This property of 3D imaging is often used in MRI of flow, because 
short TE generally reduces artifacts caused by flow-induced dephasing 
of spins (see Chapter 6). The major limitation of 3D imaging is its 
relatively long scan time, T scan , which is defined by the chosen TR 
and the total number of phase encodings: T scan = M x Q x TR. 

3.4. K-SPACE 

Tweig [11] and Ljunggren [12] have suggested a very useful inter- 
pretation of the signal acquisition in Fourier MRI as sampling along 
trajectories, k = {k x (t),k y (t),k z (t)}, in a spatial frequency space, or 
k-space. The fe-space trajectories are described by the following 
equations: k x (t) = 7 j l G x dt, k y {t) = 7 f G y dt and k z {t) = 7 /d G z dt. 
For example, in 2D imaging one of the components of k produced by 
the phase-encoding gradient is constant during acquisition, whereas 
another component produced by the readout gradient changes linearly 
with time. The corresponding trajectories in £-space are represented by 
a number of parallel lines (Figure 3.7). Several frequently used sam- 
pling schemes that allow rapid acquisition of MR images using a 
small number of &-space trajectories such as echo-planar, fast spin- 
echo and spiral imaging are discussed in Chapter 7. 

In general, arbitrary trajectories in &-space can be produced by 
time-varying gradients applied during signal acquisition, assuming 
compliance with hardware limitations such as maximum available 
gradient amplitude, G max , and maximum slew rate, (dG/dt) max . How- 
ever, the time during which the signal can be acquired along a 
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Figure 3.7. K-space trajectories in conventional 2D imaging with 
frequency and phase-encoding gradients in the x and y directions, 
respectively. 

chosen trajectory is limited by the effects of T 2 decay and magnetic 
field nonuniformity. In actual practice, these two effects can be signifi- 
cant enough during MR scanning to cause image distortions that can 
vary from mild to severe depending on the chosen sampling scheme. 

The signal amplitude at low spatial frequencies in the central 
region of /z-space is typically much higher than the signal amplitude 
at high spatial frequencies on the periphery of fe-space. This is demon- 
strated in Figure 3.8, which shows an image of a uniform phantom and 
the corresponding absolute values of the acquired &-space data. In 
general, low spatial frequencies are needed to reconstruct the bulk of 
the object, and high spatial frequencies are required to accurately 
reproduce rapid spatial variations in the object (e.g., sharp edges). 
The range of spatial frequencies in &-space is inversely proportional 
to spatial resolution. For example, in 2D imaging with in-plane resolu- 
tion 6x x Sy the signal is sampled inside a (2ir/6x) x (2ir/8y) rectangle 
in fe-space. Because high spatial frequencies beyond a certain limit 
are not acquired, reconstructed images often suffer from truncation 
artifacts (also known as Gibbs ringing) characterized by a series of 
ripples originating at intensity discontinuities (Figure 3.9). 

The appearance of truncation artifacts as ripples in images can be 
explained by the oscillations of PSF (Figure 3.1). In the regions with 
slow varying magnetization, the result of the convolution of the trans- 
verse magnetization in the object and PSF will approximate M xy well. 
In contrast, when sharp spatial changes (discontinuities) in M xy are 
present, the oscillations in PSF will cause oscillations (ripples) in 
image intensity propagating away from the discontinuities [3]. As 
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Figure 3.8. (a) Plot of the magnitude signal, \S(k x ,k y )\, produced by a 
homogeneous spherical phantom; (b) reconstructed image of a trans- 
verse section of the phantom. 




Figure 3.9. Sagittal images of the head: (a) 256 (readout) x 192 (phase) 
matrix size and 0.9 mm x 1.25 mm (in-plain) spatial resolution; (b) 256 
(readout) x 64 (phase) matrix size and 0.9 mm x 3.75 spatial resolution. 
The image in (b) displays significant truncation artifacts due to low 
resolution in the phase direction. 
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spatial resolution increases ripples with significant amplitudes move 
closer and closer to discontinuities creating an impression of 
decreased Gibbs ringing. 
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CHAPTER 4 



Contrast in MR Imaging 



The observer's ability to differentiate between different structures 
in images depends on image contrast. MRI has become an important 
diagnostic tool in medical imaging because it provides the necessary 
contrast between various soft tissues required to identify pathologic 
processes. Image contrast, C(A,B), between two different structures, 
A and B, within the object may be defined as follows: 

CtA.Bl-Ky 1 '", (40.1) 

where 1(A) and 1(B) are the image intensities of A and B, respectively; 
I re f is an arbitrary reference intensity [1]. In the case of MRI, there is no 
standard I re f, therefore, the term "image contrast" is often used to 
characterize the difference between the intensities 7(A) and 1(B). 

Variations in local relaxation times 7\ and T 2 , and the amount of 
hydrogen nuclei per unit volume (known as proton density or N(H)) are 
the predominant sources of contrast in proton MRI. The use of local 
relaxation times for identification of pathology by MRI was inspired 
by the observation that malignant tissue can have longer 7\ and T 2 
than normal tissue [2]. Although the relaxation properties of J H 
nuclei in tissue are not completely understood, there is evidence that 
the pathological states of different tissues are characterized by specific 
biochemical processes that can alter relaxation times and proton 
density [3,4]. 

In this chapter we consider different mechanisms of contrast in 
1 H MRI as well as techniques used to improve contrast in MR images. 
The emphasis of the following discussion is on the relationship between 
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image contrast and NMR relevant tissue parameters, such as inherent 
relaxation times and proton density. 

4.1. MAIN CONTRAST MECHANISMS IN X H MRI 

The dependence of MRI contrast on intrinsic relaxation times 
and proton density can be studied by using the equation describing 
steady-state transverse magnetization (M tr ) in a given pulse sequence. 
For simplicity, we will focus our discussion on spin-echo imaging with 
90 degree excitation pulses. By using Eq. (2.1.6) we obtain 

M tr = M e~ TE ^(l + e~ TR ^ - 2e -( TR -™ T >). (4.1.1) 

Taking into account that M oc N(H) and the fact that spin-echo 
images are typically acquired with TE < 27\, we can express image 
intensity, /, as 

Io,N(H)e- TE ^(l-e- TR ^). (4.1.2) 

According to this equation, image intensity is directly proportional to 
proton density, N(H), and 7\ and T 2 related modulations of the 
intensity are denned by the ratio of TR to 7\ and TE to T 2 , respectively. 
The proton density variations among soft tissues are typically within a 
few percent (except for brain tissues), whereas 7\ and T 2 relaxation 
times in different tissues vary by more than 100% (see Table 4.1). It 
is this large degree of variation in relaxation times that creates the 
greatest contrast in MRI. 

Tj Contrast 

When tissues with different values of 7\ relaxation time experi- 
ence a train of excitation pulses, tissues with short 7\ will restore 
their longitudinal magnetization more rapidly following an excitation 

Table 4.1. Approximate values of relaxation times at 1.5 Tesla [4]. 
Tissue Ti (ms) T 2 (ms) 



Gray matter 


920 


White matter 


790 


Cerebrospinal fluid (CSF)* 


2650 


Kidney 


650 


Liver 


490 


Skeletal muscle 


870 



* ^ and T 2 values for CSF are from [1]. 
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pulse and will have greater steady-state transverse magnetization than 
tissues with long T lt provided that variations in T 2 and N(H) are 
sufficiently small. As a result, tissues with short 7\ seem brighter in 
MR images with dominant T x contrast, known as 7\ -weighted 
images, than tissues with long 7\. For example, the relatively low 
image intensities of the CSF and gray matter in Figure 4.1 are due to 
their longer 7\ as compared to Tj in the white matter. 

7\ -weighted images are normally acquired with TE less than T 2 
in the tissues of interest in order to reduce T 2 modulation of the 
signal. T x contrast can be further enhanced by adjusting the repetition 
time TR. For example, by examining results from Figure 4.2 we can 
conclude that 7\ contrast between tissues with close 7\ values 
diminishes at either very short TR < 7\ or very long TR > T 1; 
whereas the maximum contrast is achieved at the intermediate 
values of TR « T v 

T 2 Contrast 

Variations in T 2 relaxation time among different tissues can be used 
to produce images with predominant T 2 contrast. These images are 
characterized by the bright appearance of tissues with long T 2 and 
dark appearance of tissues with short T 2 . Because T 2 related modulation 
of image intensity is defined by the factor exp(-TE/T 2 ), it is clear that 
echo time TE is the adjustable parameter that controls T 2 contrast. In 
practice, images with predominant T 2 contrast, known as T 2 - we ightecl 
images, are acquired with TE approximately equal to or longer than 
the shortest T 2 in the tissues of interest, and with long repetition time 
(i.e., TR > Ti) needed to reduce Tj modulation of the signal. Figure 
4.3 shows T 2 -weighted spin-echo images of the brain. Notice that due 
to the relatively short T 2 of the white matter the transverse magnetiza- 
tion decays more rapidly in the white matter than in the gray matter 
and CSF. This causes the dark appearance of the white matter in the 
T 2 -weighted images of the brain. On the other hand, the very bright 
appearance of the CSF in Figure 4.3 is due to its long T 2 and high 
proton density. 

Proton Density Images 

Variations in proton density in different tissues can provide addi- 
tional contrast in MR images. Proton density "weighted" images are pro- 
duced using short echo time TE < T 2 and long repetition time TR > T x 
in order to reduce T 2 and T x modulations of the signal, respectively. 
Figure 4.4 demonstrates proton density contrast between the white 
matter and gray matter in spin-echo imaging of the brain. The apparent 
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Figure 4.1. 7\-weighted images of the brain: FOV = 24x24 cm, 
slice thickness = 2 mm, matrix size = 256 (readout) x 192 (phase encod- 
ing), TE = 20 msec. (a) TR = 150 msec, (b) TR = 450 msec, (c) TR = 
750 msec. Notice decreased contrast between brain tissues and low 
signal-to-noise ratio in (a) caused by the increased saturation of the 
magnetization in these tissues. 
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Figure 4.2. Steady-state longitudinal magnetization in spin-echo 
imaging as a function of TR. Af z>1 and M z2 correspond to an arbitrary 
chosen T x of 0.7 s and 1 s, respectively. 




Figure 4.3. T 2 -weighted brain images: FOV = 24 x 24 cm 2 , slice 
thickness = 2 mm, matrix size = 256 (readout) x 192 (phase encoding), 
TR = 3 s. (a) TE = 75 msec; (b) TE = 150 msec. 
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Figure 4.4. Proton density-weighted image of the brain: FOV = 24 x 
24 cm 2 , slice thickness = 2 mm, matrix size = 256 (readout) x 192 
(phase encoding), TR = 3 s, TE = 20 msec. 

low image intensity in white matter is primarily due to its low proton den- 
sity as compared to gray matter (see Table 4.2). Cerebrospinal fluid, 
which has higher proton density than that of white matter and gray 
matter, also has a dark appearance in the image. The explanation of 
the low image intensity of CSF is that it has a very long T\ that causes 
saturation of the CSF even at the relatively long TR of 3 seconds used 
to acquire the image in Figure 4.4. This example demonstrates that 
pure proton density contrast among all imaged tissues might be difficult 
to achieve in practice unless very long TR and short TE are used. 

A number of abnormalities such as tumor and edema are character- 
ized by increased values of T 1; T 2 , and N(H) relative to their values in 

Table 4.2. Relative proton densities in 
the brain [1]. 



Gray matter 


0.69 


White matter 


0.61 


Cerebrospinal fluid 


1.00 


edema 


0.86 



4.1. Main Contrast Mechanisms in 1 H MRI 



normal tissue [3]. Spin-echo image intensity decreases with increasing 7\ 
because of the increased saturation of spins. Conversely, increased T 2 
and proton density lead to increased image intensity. Consequently, 
the appearance of abnormalities in MR images depends on the image 
weighting. For example, brain tumors can have a dark appearance in 
T 1 -weighted images and a bright appearance in T 2 and proton density- 
weighted images as compared to surrounding tissue. In practice both 
T v , T 2 -, and proton density-weighted images are often acquired in 
diagnostic MRI in order to improve detection of pathology. 

T 2 Contrast 

T 2 -weighted images are obtained by using gradient-echo imaging 
with echo time TE comparable to T 2 in tissue (typically 30-50 ms). 
T 2 -weighted images frequently exhibit loss of signal caused by dephas- 
ing of spins in a non-uniform magnetic field, and thereby have inferior 
contrast between soft tissues relative to 7V we ighted spin-echo images. 
On the other hand, the sensitivity of T 2 -Vfeighted images to the 
presence of magnetic field inhomogeneities has proven to be extremely 
useful for functional MR imaging (FMRI) of the brain. 

One of the main physical effects exploited in functional magnetic 
resonance imaging is the dependence of the magnetic properties of 
hemoglobin in red blood cells on the blood oxygen level: deoxyhemo- 
globin is paramagnetic, and oxyhemoglobin is diamagnetic. Ogava 
and his collaborators [5,6] have demonstrated that the difference in 
susceptibility between venous blood having a high concentration of 
deoxyhemoglobin and surrounding brain tissue causes dephasing of 
nuclear spins, which in turn gives rise to the dark appearance of veins 
in gradient-echo images. The resulting image contrast is known as 
blood oxygen level dependent contrast (BOLD). FMRI based on BOLD 
contrast [7-9] employs T 2 -weighted images to record changes in the 
signal caused by neuron activation due to various stimuli (e.g., photic 
stimulation, finger movements etc.). During neuron activation, blood 
flow in the activated brain regions increases to satisfy the increased 
demand for oxygen. The hypothesis is that the increase in blood flow 
is accompanied by an increase in oxygen supply that exceeds the 
increase in oxygen consumption. Because of the temporary imbalance 
between the supply and consumption of oxygen, the concentration of 
deoxyhemoglobin in the activated regions decreases relative to that in 
the resting state. As a result, the susceptibility difference between 
blood and its surroundings also decreases, which in turn causes the 
observed increases in intensity in T 2 -weighted images of the brain. 
Since changes in image intensity during activation are small (a few 
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percent at 1.5 Tesla), it is normally necessary to acquire a number of 
images during the stimulation and resting states and to use statistical 
analysis to differentiate between the activated and non-activated 
brain regions. 

4.2. CONTRAST AGENTS 

Certain materials, known as contrast agents, can enhance MR 
image contrast by altering T lt T 2 , and T 2 relaxation times. Contrast 
agents are frequently used in diagnostic MRI in order to achieve 
better assessment of local physiologic and anatomic conditions or to 
improve detection of malignancy. Because contrast agents are usually 
administered internally, they must possess low toxicity and be easily 
excreted from the body. 

Unlike contrast agents used in nuclear medicine and clinical 
radiography, MRI contrast agents affect the signal indirectly via 
interaction with hydrogen nuclei. Most of the commonly used MRI 
contrast agents significantly alter 7\ and T 2 relaxation times in 
tissue because of dipole-dipole interaction with water protons. Based 
on the theoretical model by Bloembergen, Purcell, and Pound [10], 
changes in T x and T 2 caused by a contrast agent can be described by 
the following equations: 

i/r.-i/r,, + «,»., 

l/T 2 = l/T 2fl + R 2 n a , 

where Tj and T 2 are the observed relaxation times; T 10 and T 2 are the 
relaxation times in the absence of the agent; n a is the agent's con- 
centration; i?i and R 2 are the agent's relaxivities. 

Paramagnetic ions with a number of unpaired electrons such as 
Gd 3+ and Mn 2+ are frequently used as MRI contrast agents because 
of their high relaxivities. To reduce the toxicity of these ions for in 
vivo applications, they are usually chelated to special molecules. Gd- 
DTPA, the best known MRI contrast agent, is a chelate composed of 
gadolinium ions and diethylenetriaminepentaacetic acid (DTP A). 
Because chelation of gadolinium ions decreases the effective amount 
of water molecules that can directly interact with the paramagnetic 
ions, the i?i and R 2 relaxivities of Gd-DTPA are smaller than the 
corresponding relaxivities of free gadolinium ion. The reported i? x 
and R 2 relaxivities of Gd-DTPA at 0.5-1.5 Tesla are approximately 
4.5kg -mmol' 1 • s" 1 and 6.0 kg -minor 1 • s" 1 , respectively [11]. Gd- 
DTPA is typically administered at a dose of 0.1-0.2 mmol/kg. 



4.2. Contrast Agents 




Gd-DTPA (mmoiykg) 

Figure 4.5. Relative change in spin-echo image intensity versus Gd- 
DTPA concentration. In this example precontrast 7\ and T 2 are 
800 msec and 80 msec, respectively. TR = 400 msec and TE = 20 msec. 



Changes in MR image intensity after administration of gadolinium 
chelates occur due to the effects of T x and T 2 shortening [11-13]. These 
two effects are competing: T x shortening leads to increased image 
intensity and T 2 shortening causes a decrease in intensity. The result- 
ing image intensity is a nonlinear function of contrast agent concen- 
tration in tissue and is dependent on other parameters such as 
intrinsic tissue relaxation times, TR, TE, and flip angle. From the 
results shown in Figure 4.5, it follows that image intensity initially 
increases with increasing concentration of Gd-DTPA because of the 
dominant effect of T x shortening. As the Gd-DTPA concentration 
continues to increase, the effect of T 2 shortening eventually becomes 
dominant causing a net reduction in image intensity. 

Gadolinium-enhanced MRI has been used successfully for the 
detection of tissue abnormalities. For example, in the brain, adminis- 
tered Gd-DTPA usually remains in the vasculature in normal tissue 
but can extravasate into the interstitial space in malignant tumors 
because of the locally disrupted blood-brain barrier [14]. This allows 
detection of brain tumors based on the observed intensity changes in 
post-contrast brain images. Another example is contrast-enhanced 
MRI of the breast. A number of studies have demonstrated intense, 
rapid signal enhancement in malignant breast tumors as compared to 
the enhancement in normal tissue [15,16]. Figure 4.6 shows pre- and 
post-contrast gradient-echo images of a patient with breast carcinoma. 
Notice significant signal enhancement in the tumor after administration 
of Gd-DTPA, whereas only relatively small changes are noticeable in 
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Figure 4.6. Contrast-enhanced MRI of the breast. The patient is a 74- 
year-old woman diagnosed with infiltrating ductal carcinoma of the 
right breast: (a) pre-contrast and (b) post-contrast sagittal MR 
images of the right breast. Imaging parameters: FOV = 18 x 18 cm 2 , 
slice thickness = 2 mm, imaging matrix — 256 x 128, 77? = 20 msec, 
77? = 4 msec, flip angle = 45 degrees. Tumor lesion in (b) is indicated 
by an arrow. 

normal tissue. Other clinical applications of gadolinium-based contrast 
agents include cardiac, vascular, and perfusion MR imaging [14]. 

Susceptibility Agents 

Contrast agents that significantly alter the local magnetic field in 
the specimen are frequently referred to as susceptibility agents. Such 
agents typically include different iron oxides [17]. Macroscopic 
samples of iron are characterized by the existence of magnetic 
domains. In their natural state the magnetic moments of these domains 
are oriented randomly, but in the presence of an external magnetic 
field, magnetic moments are aligned along the field. Iron oxide crystals 
used as MR contrast agents are extremely small (<50nm); because of 
their small size, individual crystals usually contain a single magnetic 
domain. Similar to ferromagnetic materials, single-domain particles 
align their magnetic moments along an external field. As a result, 
magnetic susceptibility in a system of such particles usually far exceeds 
the susceptibility of paramagnetic materials. On the other hand, similar 
to paramagnetic materials, the net magnetization in a system of single- 
domain particles vanishes after the magnetizing field is removed. To 
characterize the unique magnetic properties of the single-domain crys- 
tals of iron oxides, they are often referred to as superparamagnetic. 



4.3. Magnetization Transfer Contrast 



Susceptibility agents are characterized by relatively small R 1 
relaxivity compared to paramagnetic chelates. The main effect pro- 
duced by susceptibility agents is markedly increased local magnetic 
field gradients, resulting in significantly reduced T\ relaxation 
time and increased attenuation of the NMR signal due to diffusion 
[18,19]. 

4.3. MAGNETIZATION TRANSFER CONTRAST 

In 1989, Wolff and Balaban discovered a novel form of MRI con- 
trast known as magnetization transfer contrast (or MTC). MTC is 
based on the effect of magnetization transfer between two distinct 
pools of protons found in biological tissues [20,21]. These pools can 
be described as: a) freely mobile water protons Ciif) characterized by 
relatively long T 2 (>10ms); and b) restricted motion macromolecular 
protons Cflj.) in proteins and lipids with very short T 2 (<100/xs) 
(Figure 4.7). The magnetization transfer effect that exists due to 
dipole-dipole interaction and chemical exchange between the mobile 
and macromolecular protons, can be described by the modified Bloch 
equations for the longitudinal magnetizations of the X H^ and *!!,. 
pools [22,23]: 

dM z , M . f -M l£ _ kfM ^ + KM ^ 



dt T u 

dM z , M^-M^ T 
~dt~ 'r lr 



k r M z , + k f M zI , 



where indices f and r denote mobile and restricted motion protons, 
respectively; M is the equilibrium magnetization; k f and k r are the 
relaxation rate constants describing the coupling between the two 
pools of protons. 



H f ,linewidth< 100 Hz 




H r ,linewidth > 10 kHz 



Figure 4.7. Schematic of 1 H^ and x H r tissue spectra. 
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The main contribution to the observed NMR signal comes from 
the mobile protons; the signal from the macromolecular protons is 
normally negligibly small due to their short T 2 . Magnetization transfer 
affects the signal indirectly via interaction between the two pools of 
protons. In general, attenuation of the signal depends on the degree 
of saturation of the 1 K r pool. However, the basic features of MTC 
can be described by assuming completely saturated macromolecular 
protons (i.e., M zr = 0). In the latter case we obtain from Eq. (4.3.1) 

dM zf M of -M zf 

The solution of this equation is given by 

M Ztf = M 0/ ^(l - e- l ' T ^) + M zJ (t = 0)e^ T —, (4.3.3) 

where T happ = T 1:f /(1 + k f T lf ). Equation (4.3.3) shows that longitudi- 
nal relaxation of M zf in the presence of magnetization transfer is 
defined by the apparent relaxation time T liapp , which is shorter than 
Tij. The exchange between the two pools also results in a decreased 
value of the asymptotic magnetization: 

M s/ = M 0/ %^ = -§/ (4.3.4) 

where M s f = M z j(t — > oo). 

Using equation (4.3.3) in the case of repetitive r.f. excitations with 
flip angle a and repetition time TR (considered previously in Chapter 
2), we can obtain the following equation for the steady-state longitudi- 
nal magnetization (see Eq. (2.1.3)): 

M zJ (TR) = M s/ l '— (4.3.5) 

(1 -e IK ' I ^pp cos a) 

From this equation it follows that the signal at time TE after excitation 



can 

S oc M sJ e-™^ sina (1 :4 /Ti ' } (4.3.6) 

(1 - e iK ' J| -wcosa) 

By examining Equations (4.3.4) and (4.3.6), we find that the effect of 
magnetization transfer on the observed signal depends upon the 
product of kf and 7\. This result is illustrated in Figure 4.8, which 
shows image intensity as a function of flip angle at different values 
of kfT x . It is interesting that the Ernst angle that corresponds to the 



4.3. Magnetization Transfer Contrast 




' 20 40 60 80 

Flip angle (degrees) 

Figure 4.8. Image intensity versus flip angle shown for T x jTR = 20. 
Solid line: k f T x = 0, dashed line: k f T x = 0.5, dotted line: k f T x = 1.0. 



maximum signal from the l YL f pool increases in the presence of 
magnetization transfer due to shortening of the apparent T x . 

Wolff and Balaban, and other investigators have experimentally 
demonstrated that magnetization transfer between mobile protons 
and macromolecular protons can significantly decrease NMR signals 
from certain tissues such as skeletal muscle, and gray and white 
brain matter [20,24]. Because the magnetization transfer effect is 
weak in some tissues such as blood, probably due to low concentration 
of macromolecules, and stronger in others such as muscle, it can 
enhance contrast in MRI. For example, the magnetization transfer 
effect can be used to improve visualization of blood vessels in MRI 
by decreasing the signal from brain parenchyma without significantly 
affecting the signal from blood [24]. Reference [25] is recommended for 
a more detailed discussion of MTC and its applications in diagnostic 
MRI. 

MTC in MR imaging is typically achieved by selective saturation 
of the 1 H r pool. In the original approach by Wolff and Balaban [20] the 
1 H r pool is saturated by continuous off-resonance r.f. irradiation with- 
out affecting the 1 Hf pool. In another technique selective saturation of 
the 1 H r pool is achieved by intense on-resonance r.f. pulses [24,26]. A 
basic 2D gradient-echo pulse sequence utilizing this technique is 
shown in Figure 4.9. The pulse sequence includes three additional 
spatially nonselective r.f. pulses that "prepare" magnetization for 
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Figure 4.9. A magnetization transfer pulse sequence. The first 
three pulses providing saturation of the macromolecular protons are 
followed by a standard 2D gradient-echo pulse sequence. 



subsequent imaging. The first 90 degree pulse excites the transverse 
magnetization in the specimen. The magnetization is subsequently 
refocused by the 180 degree pulse and eventually rotated back to 
the z-axis by the second 90 degree pulse. Because of the extremely 
short T 2 of the macromolecular protons, a short delay (~1 ms) between 
the first and second 90 degree pulses is sufficient to greatly decrease the 
magnetization of the J H r pool without significantly affecting the 
magnetization of the 1 H^ pool. 

Recently Dixon et al. demonstrated that MTC can be observed in 
conventional multislice imaging [27]. To explain the observed effect 
we need to remember that spatially-selective excitation generates the 
transverse magnetization for the 1 Hf pool only within the imaged 
slice without affecting the mobile protons in other slices. However, 
because of the extremely broad 1 H. r resonance the excitation pulse 
causes saturation of the macromolecular protons in these slices. As a 
result, the mobile water protons lose some of their longitudinal 
magnetization via interaction with the macromolecular protons, and 
subsequent excitations of these slices therefore produce less signal. 
The resultant effect is a small decrease in intensity in multislice 
imaging as compared to single slice and 3D imaging [27]. 
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4.4. DIFFUSION- WEIGHTED IMAGES 

An additional attenuation of the signal in proton MRI results from 
diffusion of water molecules (see Chapter 1). Different microscopic 
structures such as cell membranes play the role of microbarriers that 
can significantly alter random motion of water molecules in tissue. 
Because of these microbarriers diffusion of water molecules in tissue 
is usually restricted to a limited region in contrast to self-diffusion in 
pure water. This explains why the reported diffusion coefficients in 
different tissues are 2 to 3 times smaller than the diffusion coefficient 
in water (D w 2.2 x 10" 5 cm 2 /s at 25°C). Because of the restricted diffu- 
sion the mean square displacement, (r 2 ), of water molecules in tissue 
does not increase linearly with time as predicted by the Einstein equa- 
tion, (a- 2 ) = 6Dt. Since the probability of encountering microbarriers 
increases with time, the diffusion coefficient tends to decrease at 
longer times. Moreover, different tissue microstructures often exhibit 
anisotropy, which in turn results in anisotropy of diffusion (i.e., depen- 
dence of diffusion on direction) [28,29]. 

Since signal attenuation due to diffusion is tissue-dependent, it can 
provide useful image contrast in MRI. Diffusion-weighted imaging 
(DWI) is typically performed by using the Stejskal-Tanner method [30]. 
A basic DWI spin-echo pulse sequence shown in Figure 4.10 incorporates 
two strong gradients providing diffusion-related attenuation of the 
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Figure 4.10. 2D DWI spin-echo sequence with diffusion-weighting pro- 
vided by two strong gradients applied in the z-direction symmetrically 
with respect to the refocusing 180 degree pulse. In principle, diffusion- 
sensitizing gradients can be applied in any direction. 
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signal. Diffusion weighting of the signal in the Stejskal-Tanner 
sequence is given by the factor e~ bD , where b is dependent on the 
gradient strength, duration and the time interval between the gradients 
(see Chapter 1). Therefore, when image contrast is primarily defined by 
diffusion, tissues with large values of D have darker image appearance 
than tissues with smaller values of D. 

MR images acquired with different diffusion weighting can be 
used for in vivo quantitative measurements of the diffusion coefficient 
in different tissues. Since in practice it is often difficult to discriminate 
between diffusion and other motions (e.g., blood flow) also contri- 
buting to signal attenuation, the diffusion coefficient measured in 
vivo can differ from the actual diffusion coefficient. To distinguish 
between the measured and the actual diffusion coefficients, the 
former is often termed the apparent diffusion coefficient, or ADC [31]. 

Because it uses strong diffusion gradients, DWI frequently suffers 
from artifacts caused by eddy currents. The mechanism of these artifacts 
can be understood from Figure 4.11. Distortions of the gradient waveform 
caused by eddy currents lead to an additional attenuation of the signal. 
Since the resulting image distortions vary depending on the strength of 
the applied gradients, they tend to decrease accuracy in the calculation 
of the diffusion coefficient from a series of images acquired with different 
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Figure 4.11. Eddy currents produced by switching of strong diffusion- 
sensitizing gradients cause distortions of the ideal gradient waveform 
(dashed line). These gradient distortions cause image artifacts and 
make it more difficult to measure diffusion coefficient accurately. 



diffusion gradients. In practice, implementation of DWI normally 
requires special self-shielded gradient coils or the use of electronic 
compensation techniques to reduce the effect of eddy currents (see 
Chapter 9). 

In a number of studies it has been suggested that DWI can be used 
for characterization of tissue perfusion and detection of neurologic dis- 
orders [29,31]. Moseley et al. [32] have demonstrated that ischemic 
brain injury is associated with an early and significant reduction in 
the ADC of water. Because of the changes in ADC, brain regions 
with a perfusion deficit often seem hyperintense in diffusion-weighted 
images. This suggests that DWI of the brain may be useful as a diagnos- 
tic tool for the early detection of stroke. 
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CHAPTER 5 



Signal-to-Noise Ratio in MRI 



The signal used for reconstruction of MR images is always 
corrupted by noise caused by randomly fluctuating currents in the 
receiver coil and the imaged object. The quality of MR images can be 
severely degraded by noise, making it difficult to distinguish between 
different structures in the object. The three most important parameters 
defining image quality are spatial resolution, image contrast, and 
signal-to-noise ratio (SNR). The first two parameters were discussed 
in Chapters 3 and 4. The third parameter, SNR, is defined as the 
ratio of the mean image intensity in a chosen region of interest 
(ROT) to the square root of the noise variance: 

_ mean intensity 
%/noise variance 

SNR in MRI depends upon a number of factors including strength 
of the static magnetic field, type and characteristics of r.f. coils, 
imaging parameters (e.g., image resolution and matrix size), and 
chosen pulse sequence. The following discussion focuses on the rela- 
tionship between SNR and imaging parameters as well as the field 
strength [1-8]. For a discussion of SNR in different pulse sequences 
and SNR performance of r.f. coils, see references [1,9-16]. 

5.1. NOISE VARIANCE IN MR IMAGES 

2D Imaging 

The induced voltage in the NMR receiver coil is composed mainly 
of components with frequencies distributed in a narrow range around 
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the frequency u> = jB . Acquisition of the NMR signal is normally 
performed by using phase-sensitive detectors (see Appendix) that shift 
the signal frequencies by +/ - u . The output of a phase-sensitive detec- 
tor is a sum of low and high frequency components centered at zero and 
2u Q , respectively. The high frequency components of the output signal, 
V, are subsequently removed by passing the signal through a low-pass 
filter,^. This process can be described as a convolution 

S(t, m)= f $(t - t') V{t\ m)dt\ (5.1.1) 

where index m defines a phase-encoding cycle, m = —M/2, . . . , M/2 — 1. 
By using an analog-to-digital converter the resultant signal, S(t, m), is 
subsequently sampled at a rate of 1/t w . 

In the presence of noise, the output of a phase-sensitive detector 
can be written as 

V(t,m) = V s (t,m) + V(t,m), (5.1.2) 

where V s denotes the true signal and V(t, m) denotes the noise part of 
the signal. Following phase-sensitive detection, the low-pass filtering 
of the output signal removes all signal components with frequencies 
higher than a certain cutoff frequency. To choose the cutoff frequency, 
we need to determine the range of frequencies in the part of the output 
signal that consists of low frequency components. In 2D imaging the 
true signal is described by the following equation (see Chapter 3 and 
Appendix): 

V s (t, m) oc L sl ff M xy {x,y) exip(j-fG x xt + hG y {m)yt ph ) dxdy. (5.1.3) 

In this equation M xy is the transverse magnetization; G x and G y are 
the frequency and phase-encoding gradients, respectively; t ph is the 
duration of the phase-encoding gradient; L sl is the slice thickness. 

Under the circumstances typical for 2D imaging the transverse 
magnetization is zero outside of the chosen field-of-view (FOV), 
L x x L y , in the x-y plane. That is, 

M xy {x,y) = 0, if |x| > L x /2 and/or |y| > L y /2. (5.1.4) 

By using Equations (5.1.3) and (5.1.4) it can be shown that the spectral 
components of the signal V s are zero for all frequencies \v\ > v s /2, 
where v s = 1/t w — \^G x L x /2tx\ is known as the sampling bandwidth. 
Therefore, the true signal remains unchanged after passing through 
a low-pass filter that has a cutoff frequency of i/ s /2: 

*W=^. (5-1.5) 



5.1. Noise Variance in MR Images 



In 2D imaging with N readout points and M phase encodings, the 
noise contribution to a N x M image intensity can be written as 

" ^^^-^^-''^^in', m'). (5.1.6) 

In this equation n and n' change from -N/2 to N/2 - 1; to and to' change 
from -Af/2 to Af/2 - 1; S mise (ri,m) = J$(nT dw - t')V(t',m')dt' is the 
noise signal after low-pass filtering and digitization. To calculate the 
noise variance, we consider white noise that has zero mean and auto- 
covariance 

l 0, TO ^ TO 

where .E[x(£)] is the expected value of x(t), a is a constant, and 6 is the 
Dirac delta function. From Equations (5.1.5-5.1.7) it follows that the 
noise variance in an image can be expressed as 



Var{r,) = E[\r,(n,m)\ 2 ] 

-^$>p(-; 27 m(n'-»")/A0^ 



Taking into account that the interval between sampling points, t w , is 
the inverse of the sampling bandwidth, u s , we finally obtain 

v ^=jrk = JTTS (5 - L9) 

where T s = N/ v s is the total acquisition time. Therefore, the noise level 
measured by y/ Var(rj) is inversely proportional to the square root of the 
acquisition time. Notice that the noise variance in Eq. (5.1.9) is spatially 
independent. In principle this allows measurements of Var(rf) in an 
arbitrary region within the field-of-view if this region contains a 
sufficient number of pixels making statistical averaging possible. In 
actual practice it is often convenient to measure noise in the back- 
ground outside of the object's boundaries where noise is the only 
contributor to image intensity. 

3D Imaging 

In the case of 3D imaging with N readout points and a total of 
M x Q phase encodings, the output of the phase-sensitive detector 
can be written as 

V(t, to, q) = V s (t, to, q) + V(t, to, q). (5.1.10) 
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In this equation index q indicates an additional phase-encoding used 
for 3D imaging: q = -Q/2, . . . , Q/2 - 1 (see Chapter 3). After low-pass 
filtering and digitization, the noise signal is given by 

Snoisein', m'q') = f §{n'T dw - t') V(t\ m',q')dt'. (5.1.11) 

The noise part of a N x M x Q image is given by a three-dimensional 
Fourier transform of S noise : 

^ n' m' q' 

(5.1.12) 
In the case of 3D imaging, the noise autocovariance can be expressed as 

L 0, m ^ m and/or q ^ q 

(5.1.13) 

By using Equations (5.1.11) through (5.1.13) it can be shown that the 
noise variance in 3D imaging is given by 

Var{r)) = 



~NMQ M-Q-T.' v ' ; 

5.2. SNR IN 2D AND 3D MRI 

Based on the results of the previous section, we consider SNR 
dependence on the basic imaging parameters such as voxel volume, 
matrix size, and acquisition time in 2D and 3D imaging. The obtained 
results demonstrate that 3D imaging in general provides higher SNR 
as compared to 2D imaging. 

2D Imaging 

By using the results from Chapter 3 and assuming for simplicity 
that M xy is uniform within the excited slice we obtain 

(5-2.1) 

where / is the image intensity in the absence of noise and 
L sl -L x -L y /(N-M) is the voxel volume, 6V. From Eq. (5.2.1) and 
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Figure 5.1. Images of an apple in (a) and (b) were obtained with the 
same parameters except for FOV, which was reduced by a factor of 4 
in (b) as compared to (a). The reduction in FOV decreased the SNR 
and lowered the overall image quality in (b). 



Eq. (5.1.9) we obtain that the SNR in 2D imaging is given by 1 



SNR oc MsySVy/ 



(5.2.2) 



This equation shows that SNR is proportional to voxel volume. 
Although decreasing voxel volume generally improves visualization 
of small structures in images, it causes a proportional decrease in 
SNR. It should be noted that the decreased SNR can result in unaccep- 
tably low image quality, defeating the initial purpose of the decreased 
voxel size (Figure 5.1). 

The SNR dependence on voxel volume, number of phase-encoding 
steps, and acquisition time should be kept in mind when choosing 
imaging parameters. For example, if one decides to decrease L x and 
L y by a factor of 2 (without altering matrix size and T s ) in order to 
increase in-plane spatial resolution, then the SNR will be decreased 
by a factor of four. Conversely, if the same increase in spatial resolu- 
tion is achieved by correspondingly increasing matrix size (e.g., 
using N' = 2N and M = 2M) while maintaining FOV and T s , then 



1 Note that the SNR in Eq. (5.2.2) is a complex quantity. In practice SNR is frequently 
calculated in magnitude MR images. The characteristics of noise in magnitude images 
are discussed in [17]. 
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the SNR in the image will be reduced only by a factor of 2\/2 at the 
expense of a twofold increase in scan time. 

3D Imaging 

Image intensity in 3D MRI is proportional to the product of M xy 
and the voxel volume, 8V = L x ■ L y ■ L s! /(N ■ M ■ Q). By using Eq. 
(5.1.14), the expression for the SNR in 3D imaging can be written as 



SNR oc M xy SVy/MQT s . (5.2.3) 

After comparing Equations (5.2.3) and (5.2.2) one can conclude 
that the SNR in 3D imaging is ^/Q times greater than the correspond- 
ing SNR in 2D imaging. This result is due to more statistical averaging 
of the noise in 3D imaging as compared to 2D imaging. The actual 
increase in SNR in 3D imaging can be very significant. For example, 
when imaging with the same voxel volume and acquisition time, the 
SNR in 3D imaging with Q = 64 will be eight times greater than the 
SNR in 2D imaging. 

The Effect of Signal Averaging on SNR 

Suppose that NSA identical (except for noise) frequency and phase- 
encoded signals are acquired and subsequently averaged together. 
Because averaging reduces the noise variance by a factor of NSA 
without altering the true signal, it is clear that averaging increases 
SNR by a factor of (NSA) l/2 . 

The number of times the signal can be averaged in 2D imaging is 
given by 

NSA ^^Ffn' (524) 

where T scan is the total scan time, TR is the repetition time, M is the 
number of phase-encoding steps, N si is the number of slices, and N' s i 
is the number of slices imaged during one TR interval. By multiplying 
the expression for the SNR in Eq. (5.2.2) by (NSA) 1/2 we obtain 



T< ■ T arnn 

s scan (525) 



y/N tl ■ TR 

We would like to compare the SNR in 3D imaging with that in 2D 
imaging under the condition that the scan time is the same. The 
number of times the signal can be averaged in 3D imaging is 
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By combining Equations (5.2.6) and (5.2.3), the SNR in 3D imaging can 
be expressed as 



— . — -. (5.2.7) 

y/TR 

In summary, the SNR in 2D and 3D imaging is proportional to 
the voxel volume and the square root of the acquisition time. Three- 
dimensional imaging in general has higher SNR than 2D imaging 
given the same M xy , voxel size and acquisition time. The SNR in 2D 
and 3D imaging increases as the square root of NSA at the expense 
of prolonged scan time (which is proportional to NSA). After examin- 
ing Eqs. (5.2.5) and (5.2.7), we conclude that in the presence of signal 
averaging, the SNR in 2D imaging approaches that in 3D imaging as 
the number of slices imaged during a single TR interval, N' sh 
approaches the total number of slices, N sl . 

5.3. SNR FIELD DEPENDENCE 

Experience shows that SNR in MR images varies depending on the 
strength of the static magnetic field, B . Because of its importance for 
MRI, the SNR dependence on the field strength has been studied by 
many investigators [1-4]. The following approach, used to derive the 
SNR field dependence, mainly employs the arguments put forward by 
Hoult et al. [1,2]. 

To describe the SNR dependence on B , we need to consider how 
the NMR signal and noise variance depend upon B . As shown in 
Chapter 1, the signal is proportional to the square of B . By examining, 
Equations (5.1.9) and (5.1.14) the reader can realize that the depen- 
dence of the noise variance on B is "hidden" in the constant a. 

Because the source of noise in MRI is randomly fluctuating currents 
in the sample and the receiver coil, a is given by the Johnson formula: 

a 2 = 4kT(R sample + R coil ). (5.3.1) 

In this equation, k is Boltzmann's constant, T is the temperature, and 
Rsampie an d R co u are tne sample and coil resistances, respectively. Due 
to the "skin effect," the coil resistance is proportional to the square 
root of B [1]. To determine the field dependence of R samp i e we need 
to consider thermal losses in the sample. 

The average power of energy losses in the sample is given by 

P= J Xs (E 2 )dV, (5.3.2) 
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where \ s is the conductivity of the sample, E is the induced electric field, 
and (E 2 ) denotes the average value of E 2 . On the other hand, the expres- 
sion for P can also be written as 

P = (I 2 )R samp ie, (5.3.3) 

where I is the oscillating current used for excitation of the transverse 

magnetization in the sample, and (I 2 ) denotes the average value of/ 2 . 

Assuming that I — I cos ujt, the induced magnetic field in the sample is 

B(xl cos cut. (5.3.4) 

From the equation VxE = — it follows that the induced electric 

a u • cdt 

field is 

E <xcol sinut. (5.3.5) 

Therefore, (E 2 ) is proportional to the square of the r.f. frequency u. At 
resonance u is equal to the Larmor frequency of spins 7B0. According 
to Equations (5.3.2) and (5.3.5), in this case the average losses in the 
sample are proportional to the square of B ; therefore R samp i e oc B\. 

The resultant dependence of the noise variance on B can be 
described as 

Var{rj) oc a s B 2 + a c B l J 2 , (5.3.6) 

where q s and a c are two constants that depend on the geometry of the 
sample and the coil, and their conductivities. Because the true signal is 
proportional to Bq, the SNR dependence on the field strength is given 
by[l] 

SNR oc B ° = . (5.3.7) 

<Ja s Bl + a c Bl /2 

This equation shows that thermal noise at low magnetic field strength 
is primarily defined by the resistance of the coil and the corresponding 
SNR oc B 7 ' 4 . Conversely, as B increases, resistance of the sample 
eventually exceeds that of the coil 2 . In the latter case the SNR becomes 
proportional to Bo- 
lt should be noted that there exist several additional factors which 
can influence the SNR dependence on B . First, because 7\ relaxation 
time decreases with decreasing field strength [18], SNR at low field 



2 The majority of whole body MR scanners operate at field strengths between 0.3 and 
2 Tesla. The reported measurements of the SNR field dependence [3] clearly indicate 
that for these field strengths R sam pie > Rcoil- 



strength benefits from the smaller saturation of imaged material rela- 
tive to that at high field strength. Second, the dependence of the 
sample conductivity on r.f. frequency becomes important at high field 
strengths when SNR oc B Q /y/x^. That is, according to [19] changes in 
the tissue conductivity alter the linear relationship between SNR 
and field strength at B greater than 1.5 Tesla. Another factor, which 
becomes important when imaging at high field strength, is that 
penetration of r.f. field inside the object decreases with increasing 
B . The effect of a decreased r.f. penetration on image intensity and 
SNR is discussed in [20]. 
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CHAPTER 6 



Image Artifacts 



Experience shows that MR images are degraded frequently by 
different artifacts that occur for a variety of reasons. Historically, the 
quality of MR images was often low because of the artifacts caused by 
heterogeneity of the main field, eddy currents, nonlinearity of imaging 
gradients and other scanner-related problems. More recently, the fre- 
quency of occurrence and severity of these artifacts in MR images 
have been reduced significantly due to technical improvements in 
MRI hardware. A detailed description of hardware-related artifacts is 
beyond the scope of this book. Instead, the following discussion focuses 
on several major artifacts that are caused by factors predominantly 
associated with the imaged object: magnetic field nonuniformity, fat/ 
water chemical shift, flow and motion. 

6.1. ARTIFACTS DUE TO MAGNETIC FIELD 
NONUNIFORMITY 

During MR imaging the magnetic field in the object can be non- 
uniform for various reasons including nonuniformity of the external 
magnetic field, presence of metal implants and spatial variations in 
magnetic susceptibility of the object. Magnetic field nonuniformity 
causes two major types of artifacts in MR images: geometric distortion 
and signal loss due to intravoxel dephasing [1-7]. 

Geometric Distortion 

To describe geometric distortion caused by magnetic field non- 
uniformity, it is convenient to express the magnetic field in the 



Chapter 6 Image Artifacts 



specimen as a sum of two terms: 

B = B + B', (6.1.1) 

where B is uniform and B' is a function of coordinates. Because MR 
imaging is typically performed in a very homogeneous magnetic field 
such that |B'| <C |B |, the components of B' perpendicular to B can 
be neglected because they cause only a minor perturbation of 
Larmor frequencies of spins. Thus, we can assume that B' = kB', 
where k is the unit vector in the direction of B . 

For simplicity we will confine our discussion to the case of 2D 
gradient-echo imaging considered in Chapter 3. The magnetic field B' 
gives rise to the spatially-varying phase of the transverse magnetiza- 
tion <p = w't, where J = jB 1 . A general equation describing the signal 
acquired in the presence of T 2 relaxation and dephasing of spins 
caused by the field inhomogeneity B' can be written as 

S(k x , k y ) = ff M xy ( X , y)e J^ + Jk y y + jAt + TE)-(t + TE)/T 2 dx dy ^ (aL2) 

where k x = •yGj and k y = "fG y t ph ; G x and G y are the readout and phase- 
encoding gradients, respectively; t = nr w and n changes from -N/2 to 
N/2 - 1; TE is the echo time (see Chapter 3). To simplify our analysis 
we neglect T 2 relaxation and assume that B' can be approximated as 
a sum of constant and linear terms: B' = a + G' x x + G' y y + G' z z. Using 
these assumptions we obtain from (6.1.2) 

M ( x-6B'/G x \ 

S{k ^ = II " i + ajo' eik ' x+ikyy+lJTE dxd ^ ( 6 - L3 ) 

where 8B' = a + G' y y + G' z z. 

Geometric distortion is produced by intrinsic field gradients G' x , G' y 
and G' z that alter the dependence of Larmor frequencies of spins on 
their actual locations in the specimen. By examining (6.1.3) we can con- 
clude that the image distortion can be described by the following: 
spatially varying image displacement in the readout direction by 
8B'/G X , magnification in the readout direction and weighting of the 
image intensity by the factor (1 + G^/G.,.)- 1 . The strength of the read- 
out gradient is typically high enough to allow only minor geometrical 
distortion in conventional MR imaging. However in some cases, non- 
uniformity of the magnetic field can be large enough (e.g., due to the 
presence of metal objects [7]) to cause significant degradation of 



6.1. Artifacts Due to Magnetic Field Nonuniformity 



Another artifact caused by magnetic field nonuniformity is distor- 
tion of the slice profile. In Chapter 2 we have shown that spatially 
selective excitation in a uniform magnetic field creates a nonzero 
transverse magnetization in a slice of material with thickness defined 
by the bandwidth of the excitation pulse and the applied gradient. 
Assuming that 2u b is the excitation bandwidth and G z is the slice- 
select gradient, we obtain that the boundaries of the slice are: 
Zi = -Vb/dGz) and z 2 = uJb/{iG z ) (see Chapter 2). In the presence of 
a magnetic field inhomogeneity, B'(x,y,z), the slice boundaries are 
defined by the following equations: 

-^- 7 B , (x,v,2: 1 ) w b --fB'(x,y,z 2 ) 
Zl = — and z 2 = — . (6.1.4) 

From Eq. (6.1.4) it follows that magnetic field nonuniformity gives rise 
to nonplanar slices of varying thickness. The degree of slice distortion 
depends on magnetic field nonuniformity and slice-select gradient: a 
large B' or small G z can result in significant distortions of the slice 
profile. 

It should be noted that no image distortion due to magnetic field 
inhomogeneities exists in the phase-encoding direction. Consequently, 
distortion free images can in principle be acquired by using phase 
encoding in three orthogonal directions. However, this acquisition 
technique is rarely used because it is prohibitively slow. 

Intravoxel Dephasing 

Magnetic field inhomogeneities cause dephasing of nuclear spins 
during data acquisition. Spin dephasing in turn leads to a loss of 
NMR signal. The resulting effect is a noticeable reduction in image 
intensity (Figure 6.1). Mathematically, the signal loss is defined 
by the factor e JL0 TE averaged over the voxel volume. It is instructive 
to consider the effect of intravoxel dephasing in the presence of an 
isotropic gradient, G' x = G' y — G' z , in a cubic voxel, A 3 x. In this case 
the reduction in image intensity is given by the factor 

I i f^ )TE A]^m (615) 

where Aw = jG' x Ax/2. For example, an isotropic gradient of 1 ppm/cm 
causes a loss of more than 50% of the signal in gradient-echo imaging 
at 1.5Tesla with TE of 30 msec and 2 mm 3 voxel size. The good news is 
that decreasing voxel size and shortening echo time reduces loss of 
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Figure 6.1. Images of a uniform phantom in (a) and (b) were acquired 
using a gradient-echo pulse sequence with the same imaging param- 
eters. Visible reduction in image intensity in (b) is due to a metal 
paper clip placed in the vicinity of the phantom. 

signal. The bad news is that by decreasing voxel size we proportionally 
reduce SNR in images. 

Equation (6.1.5) is interesting because it clearly indicates that 
magnetic field inhomogeneities generally cause a nonexponential 
decay of signal as a function of echo time in gradient-echo imaging. 
In many instances this decay can be significant enough to cause a 
drastic reduction in image intensity (Figure 6.1). It is therefore impor- 
tant that the effect of magnetic field nonuniformity can be minimized 
by acquiring spin-echo signals. Because of the reduced dephasing of 
spins, the resulting spin-echo images generally have better contrast 
than gradient-echo images. 



6.2. CHEMICAL SHIFT ARTIFACTS 

Because of the chemical shift effect, hydrogen nuclei in water and 
fat have different Larmor frequencies even if the magnetic field in 
media is highly uniform. Fat/water chemical shift is about 3.5 ppm. 
For example, at 1.5 Tesla the difference between the resonant frequen- 
cies of water, v w , and fat, i/f, is approximately 220 Hz. Chemical shift 
between fat and water causes two different types of artifacts: (a) mis- 
registration of fat or water in MR images; (b) visible reduction in 
image intensity as a result of destructive interference between the 
fat and water signals. 
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Figure 6.2. Spin-echo images of a phantom composed of a vial with fat- 
containing vegetable oil placed in a bigger tube filled with water. 
Images in (a) and (b) were obtained at 1.5 Tesla with readout gradients 
of 0.6 gauss/cm and 0.15 gauss/cm, respectively. The increased spatial 
displacement of fat in (b) caused a signal void in the regions that 
artificially "lost" fat. 

As mentioned earlier, by using phase-sensitive detection the 
frequency of the observed NMR signals is lowered by the reference 
frequency u> . Using Eq. (6.1.3) it can be shown that if the Larmor 
frequency of water is chosen as the reference frequency, then fat in 
an MR image will be shifted in the readout direction by 
2tt6is 



6x = 



(6.2.1) 



relative to its actual location in the object. In this equation 
6v = i/ w - Vf and G x is the readout gradient. At a field strength of 
1.5 Tesla and a readout gradient of 1 gauss/cm, the misplacement of 
fat is only about 0.5 mm. However, at higher magnetic field strengths 
or lower gradient strengths the misplacement of fat can become signif- 
icant enough to cause visible image distortion (see Figure 6.2). 

In spin-echo imaging, the acquired phase difference between the 
water and fat magnetizations is canceled out following a 180° refocus- 
ing pulse. As a result water and fat are "in-phase" at echo time TE. In 
gradient-echo imaging, water and fat can be either in-phase or out-of- 
phase depending on the product of bv and TE. For example, the maxi- 
mum destructive interference between water and fat signals occurs 
when TEbv = (n + 1/2), where n = 0, 1, 2 . . . (Figure 6.3). A common 
artifact caused by destructive interference between water and fat 
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Figure 6.3. As fat and water magnetizations precess with different 
frequencies, the resulting signal exhibits oscillatory behavior. The 
signal intensity has its minima at t = r, 3t, 5t . . ., where r = 
\/2{v w - u f ) « 2.2 msec at 1.5Tesla. 

signals is an apparent signal void in voxels that contain about the same 
amounts of fat and water. This artifact often gives rise to artificial 
black boundaries between fat and muscle regions in gradient-echo 
images (Figure 6.4). 

Dixon [8] has suggested a modification of standard spin-echo 
acquisition in order to obtain images depicting only water or only fat. 




Figure 6.4. Gradient-echo image of the abdomen obtained at 1.5Tesla 
with TE of 2.2 msec. Boundaries between fat and muscle regions 
appear black due to interference between the fat and water signals. 
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Figure 6.5. Basic pulse sequence in the Dixon method. In one of two 
images, acquired with r = 0, water and fat are in-phase. Another 
image is acquired with r = l/2Sv so that the water and fat magnetiza- 
tions have opposite directions. 

In the Dixon method two images are acquired. One of these images is 
produced from signals S w+ f, obtained with a spin-echo sequence that 
ensures that signals from water, S w , and fat, S f , are in-phase; therefore, 
S w+ f = S w + Sf. Another image is produced from signals, S w _f, 
acquired during a gradient echo shifted by ±l/26is relative to the spin 
echo. In the latter case water and fat are out-of-phase; therefore, 
S w _f = S w — Sf (Figure 6.5). Subsequently the water-only image is 
reconstructed from the sum S w+ f + S w _f, and the fat-only image is 
reconstructed from the difference S w+ f — S w _f [8,9]. 

The major problem in the Dixon method is the effect of magnetic 
field inhomogeneities on the signals from fat and water. In the Dixon 
method the net signal from water and fat being out-of-phase is given 
by S w _f = (S w - Sf)e j6<t> , where 6<p is the phase shift (accumulated 
during time \/28v) due to magnetic field inhomogeneities. If 6(f> is not 
negligibly small then the Dixon approach can cause artifacts because 
of incomplete fat/water decomposition. 1 

Better differentiation between water and fat can be achieved in the so-called three-point 
Dixon method [10], which utilizes an additional measurement to separate phase shifts 
caused by magnetic field nonuniformity from those caused by the fat/water chemical shift. 
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A technique known as chemical shift selective (CHESS) imaging 
discriminates between fat and water by first selectively exciting and 
subsequently destroying the magnetization of the unwanted spectral 
component (e.g., fat) [11]. For example, the fat magnetization can be 
excited (without affecting the water magnetization) by a frequency- 
selective 90 degree pulse with narrow bandwidth applied in the 
absence of external field gradients. After the excited magnetization 
in fat is destroyed by applying spoiling gradients, the water magnetiza- 
tion is imaged. 

Another technique for selective excitation of fat or water employs 
a family of pulses, known as binomial pulses. A binomial pulse can be 
defined as a combination of equally spaced pulses such as 1-1 or 1-1, 
1-2-1 or 1-2-1, where the numbers indicate the relative pulse flip 
angles and the bar over a number indicates a 180° phase-shifted pulse 
[12,13]. For simplicity we consider selective excitation of fat by a 
pulse sequence 1-1, known as the jump-and-return pulse. Suppose 
that the first pulse has a flip angle of 45 degrees and is applied in the 
x-direction in the rotating reference frame. Fat and water magnetiza- 
tions immediately after the pulse are shown in Figure 6.6(a). If the 
interval between pulses is l/28v, the second 45 degree pulse applied 
along the negative x-axis (Figure 6.6(b)) returns the water magnetiza- 
tion to the z-axis while placing the magnetization of fat into the trans- 
verse plane (Figure 6.6(c)). 

Recently, several techniques for simultaneous spectral and spatial 
excitation have been introduced [14,15]. As in the Dixon method, mag- 
netic field nonuniformity presents a problem for selective excitation of 
fat or water. As a result these techniques, as well as CHESS and excita- 
tion by binomial pulses, require sufficient magnetic field homogeneity 
in order to achieve good discrimination between water and fat. 

Another MR approach for fat/water decomposition employs the 
difference in 7\ relaxation times of fat and water 2 in order to null 
the signal from fat while still retaining significant signal from water. 
In practice this approach is implemented by using inversion-recovery 
imaging. If T x of fat is indicated as T X f at , then the inversion time 
needed to null the signal from fat is given by (see Chapter 2) 

w - r * ta [lT=SPnyjd- (622 » 



2 The reported T x values for fat vary from 150 msec to 400 msec [16], while the typical T r of 
water-based tissues at 1.5 Tesla is greater than 500 msec. 
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Figure 6.6. Excitation of fat by a 1-1 pulse sequence: (a) fat and water 
magnetizations immediately after the first 45 degree pulse applied in 
the x-direction; (b) fat and water magnetizations immediately before 
the second 45 degree pulse applied along the negative x-axis; (c) 
immediately after the second pulse the water magnetization is along 
the z-axis while the fat magnetization is in the transverse plane. 



Since T lfat can vary depending on patient and anatomical site, the 
inversion time TI is often adjusted manually during the MR scan 
until the TI that causes the greatest decrease in the fat signal is 
found. Note that this approach is insensitive to the presence of 
magnetic field inhomogeneities because it does not rely on the fat/ 
water chemical shift. 



6.3. ARTIFACTS DUE TO MOTION AND FLOW 

Because in conventional MRI in vivo the time required for spatially 
selective excitation, phase, and frequency encoding is normally about 
100 msec or less, patient's motion during a single phase-encoding cycle 
can be considered minor. However, the time (a few minutes) needed to 
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Figure 6.7. Spin-echo image of the head. The visible image blurring was 
caused by the patient motion during imaging. 



acquire a complete set of 128 or 256 phase-encoded signals is long com- 
pared with the time scale for various physiologic motions such as 
respiration, cardiac motion, etc. As a result, motion-related changes in 
the specimen during different phase-encoding cycles create a potential 
source of artifacts in MR images. 

Motion can perturb both amplitude and phase of the nuclear mag- 
netization. Perturbations of the magnetization amplitude are caused 
by spins entering and exiting a voxel during MR imaging. A common 
result of such perturbations is image blurring, which can vary from 
very mild to very severe depending on the range and velocity of 
patient's motion (see Figure 6.7). 

To examine motion-induced changes in the phase of the magnetiza- 
tion, let us consider for simplicity one-dimensional motion of spins in 
the direction of an applied gradient (e.g., G x ). The trajectory of a 
spin can be expanded in a series 

x (t) = x(0) + v x (0)t + ^~ L + ---, (6.3.1) 
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where x(0), v x (0), and a x (0) are the position, velocity, and acceleration 
of the spin at t = 0, respectively. The motion-induced phase can be 
written as 

4>{t) = 7 *(0) f G x dt! + 7^(0) / Gj dt' + ^°i f G x t' 2 dt! ... 

+ -yx n (0) r Gxt <n dt , + (632) 

o 

The first term on the right side of this equation describes the phase 
shift that occurs in the absence of motion. The second term is caused 
by motion with constant velocity. The third and higher order terms 
result from nonstationary motion. We will discuss two important arti- 
facts resulting from the motion-dependent phase: image "ghosts" that 
occur in the presence of periodical motion of spins [17-19] and signal 
loss due to distribution of spin velocities in a voxel [20, 21]. 

Image Ghosts 

One of the most common motion-induced artifacts is the appear- 
ance of image "ghosts," which normally occurs in the presence of 
periodical motion such as cardiac motion, respiratory motion, or pulsa- 
tile blood flow (see Figure 6.8). Regardless of the specific nature of the 
periodical motion, the resultant image artifacts appear as a series of 
ghost images of the moving structures. To understand the mechanism 
of these artifacts, we will consider a specific case when image ghosts 
are generated by pulsatile flow. 

An analytical description of the ghosting artifact can be obtained 
by using a model [19] in which velocity of spins v x during phase- 
encoding cycles is given by the sum of a constant term v 0x and a 
periodical modulation. That is, we assume 

v x = v o,x + i>i,* sm(2wmf p TR + 6), (6.3.3) 

where v lx and f p are the amplitude and frequency of modulation, 
respectively; TR is the repetition time; m is an integer that defines 
a phase-encoding cycle (see Chapter 3); 6 is an arbitrary phase. Flow 
pulsations in the presence of the readout gradient G x alter the phase 
of the magnetization according to the Eq. (6.3.2). Assuming that the 
excitation pulse is applied at time t = 0, the flow-induced phase shift 
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Figure 6.8. Gradient-echo images of a phantom composed of tubing 
with flowing liquid (water with added Gd-DTPA) immersed in a 
water bath. The phantom was connected to an electric pump. Image 
in (a) depicts stationary flow in the phantom. The artificial ghosts 
seen in (b) are caused by flow pulsations with frequency of 1 Hz. 

of the gradient echo during the m lh phase-encoding cycle is given by 

fam) = fa+fa{m), (6.3.4) 

where fa = jM hx v^ x , fa{m) = ^M lx Vi x sm(2irmf p TR + 9); M liX = 
Jq G x t' dt' is the first order moment of the gradient waveform. Note 
that the second term on the right side of Eq. (6.3.4) describes the 
periodically varying phase induced by the motion of spins. We shall 
show that this periodicity causes image ghosts in the direction of the 
phase encoding. 

Taking into account the motion-induced phase of spins, the signal 
in standard gradient-echo imaging can be written as 

S(k x , k y {m)) = If M xy (x,y)e jk ' +jk ^ m)y+ ^ +jMm) dxdy. (6.3.5) 

The factor exp(jfa(m)) can be expanded in a series 

e JMm) = g J,(>yM u v ltX )e i *"' mf > TR+iW , (6.3.6) 

Z = -oo 

where J/ is a Bessel function of order I. Using the last two equations we 
obtain: 

S(k x ,k y (m)) = JJ ]T J^M^v^M^y-lSy) 

x e*x*+yM»0y+M+;w dx dy, (6.3.7) 
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Figure 6.9. An example of pulsatile ghosts in MRI. The visible image 
ghosts in a gradient-echo image of the chest were produced by blood 
flow pulsations in the aorta. 

where by = f p TRL y , L y = 2nm/k y (m) is the field-of-view in the phase- 
encoding direction. The resulting image of the structure with pulsatile 
flow contains ghosts described by the terms with nonzero I in the series 



An example of image ghosts induced by blood flow pulsations is shown 
in Figure 6.9. The distance between the ghosts, 6y, is defined by the 
frequency of velocity modulation, repetition time, and field of view. 
The relative intensity of the I th ghost, J t ('yMi tX v lx )e jW , is a function of 
the amplitude of modulation and the first moment of the gradient wave- 
form. Ghosts that are located outside of the field of view are folded back 
into the image because of aliasing; therefore, image ghosts cannot be 
removed by increasing the separation between the ghosts. If ghosts 
and other structures in images overlap, then the resultant image inten- 
sity becomes a sum of the complex intensities of the ghosts and the over- 
lapping structures. Under the condition that ghost and background 
magnetizations are out of phase, the ghost can appear darker than its 
surroundings, whereas in other cases the in-phase condition ensures 
the appearance of bright ghosts in images [19]. 
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Flow-Induced Signal Loss 

As shown earlier, motion of spins in the presence of applied gradi- 
ents modulates the phase of magnetization. Since normally there is a 
spread of spin velocities within a voxel, the resultant effect is a loss 
of signal from the voxel due to dephasing of spins. To model this 
effect we assume that moving spins have one-dimensional velocity dis- 
tribution that is uniform between —v max and +v max , where v max is the 
maximum velocity of spins. In the presence of a gradient G x applied 
in the direction of flow the signal from the moving spins is given by 

-D„ 

where M lx = $ G x t' dt' and £ is a constant that does not depend on velo- 
city. Equation (6.3.8) shows that the flow induced signal loss is small 
when \ , yMi x v max \ <C 1. For \-yMi x v max \ > 1 the signal loss becomes signif- 
icant. It is important to realize that signal loss depends on the velocity 
distribution. For example, one can show that a Gaussian velocity distri- 
bution leads to signal loss that is exponentially dependent on the square 
of the product of the first moment of the gradient waveform and the 
width of the velocity distribution. 

6.4. MRI TECHNIQUES FOR MOTION ARTIFACT 
REDUCTION 

Motion-induced artifacts such as blurring and image ghosts 
represent a major problem in diagnostic MRI because they impede 
visualization of small structures in images. In this section we discuss 
several techniques used to reduce motion artifacts in MRI. Because 
each of these techniques alone has its limitations, they can be used 
in combination to achieve better results. 



Physical Restraint 

Motion-induced artifacts can be reduced by decreasing the range 
of motion that causes the artifacts. In diagnostic MRI a reduction in 
patient motion is often achieved by using different immobilization 
devices. Suspended respiration is another popular technique used to 
eliminate respiratory motion of the patient within a breath-holding 
interval during imaging [22]. Due to its simplicity some form of 
physical restraint is frequently used during clinical MRI exams, 
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although success of this approach ultimately depends on the patient's 
cooperation. 

Signal Averaging 

Signal averaging is the approach in which signal acquisition is 
repeated a number of times for each phase-encoding. The repeatedly 
acquired signals are subsequently averaged together in order to 
reduce the contribution from different physiological motions of the 
patient. The number of times the signals are averaged together will 
be referred to as NSA. Signal averaging can be performed by acquiring 
a series of data sets, each composed of M phase encodings. This 
approach is known as serial averaging. An alternative strategy, 
known as parallel averaging, is to repeat acquisitions for each phase 
encoding before changing the phase-encoding gradient. In the case of 
periodical motion, parallel averaging is ineffective when NSA x TR 
is much smaller than the period of motion because the averaged signals 
are acquired at the same stage of the motion cycle. Serial averaging is 
generally more effective for reducing image ghosts because it allows a 
longer interval between repetitive measurements thereby increasing 
the likelihood that the averaged signals are acquired at different 
stages of the motion cycle. 

Averaging is a simple method for motion artifact reduction in 
MRI. Averaging does not require any special equipment and normally 
can be implemented with every pulse sequence. The major disadvan- 
tage of averaging is prolonged scan time, which is directly propor- 
tional to NSA. Another disadvantage is that averaging generally 
does not reduce motion-induced blurring in images. 

Respiratory Gating 

This approach is based on constant monitoring of patient position 
during respiratory motion by using a special transducer attached to 
the patient [23]. The obtained information is used either to trigger ima- 
ging pulse sequence when the patient is in a proper position or to reject 
a number of signals acquired during patient motion. Respiratory 
gating reduces ghosting and blurring artifacts because it minimizes 
motion-induced variations in the signal. The major disadvantage of 
respiratory gating is that it can double or triple scan time in diagnostic 
MRI. Scan time with respiratory gating can be reduced at the expense 
of increased blurring and ghosting artifacts if respiratory gating is 
implemented only during a fraction of data acquisition [24]. 
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Cardiac Gating 

Cardiac motion and related blood flow pulsations present a ser- 
ious problem for in vivo MR imaging with scan times longer than 
the duration of the cardiac cycle. The basic idea of cardiac gating 
is to synchronize MR imaging with the cardiac cycle so that the 
same magnetization is imaged during successive phase encodings. 
In practice, cardiac gating is implemented by monitoring the 
electrocardiogram (ECG) signal produced during the cardiac cycle 
including contraction (systole) and relaxation (diastole). The ECG 
signal (a few mV) can be measured by applying electrodes to the 
body. Alternatively, cardiac motion can be monitored by recording 
an optical signal from a peripheral pulse of arterial blood (e.g., in a 
fingertip). 

Detected ECG signal can be used to synchronize successive data 
acquisitions by repeatedly triggering a pulse sequence at a certain 
stage of the cardiac cycle. Consequently, in this approach the shortest 
repetition time, 77?, is equal to the cardiac period. One disadvantage of 
ECG triggering is that it results in a relatively long TR of about 1 sec, 
which in turn leads to long scan times in diagnostic MRI. Another dis- 
advantage is the appearance of artifacts caused by variations in 77? 
due to irregular heart beats. 

In an alternative approach, known as retrospective gating, MR 
imaging is performed continuously with fixed TR, which is much 
shorter than the period of the cardiac cycle [25]. The information 
about cardiac motion during imaging is constantly recorded. After 
completion of data acquisition, this information is used for interpola- 
tion and sorting of the data according to the corresponding phases of 
the cardiac cycle. Subsequently a number of images depicting different 
stages of cardiac motion are reconstructed. The major advantage of 
MR imaging with retrospective gating is its ability to depict the 
dynamics of the entire cardiac cycle. 



Ordered Phase-Encoding Measurements 

Ghosts that occur in MR images acquired in the presence of peri- 
odic motion are caused by periodic modulation of the magnetization 
during image acquisition. To eliminate image ghosts Perman et al. 
[19] suggested using a random order of phase encoding which, in the 
presence of periodic motion, leads to a random modulation pattern. 
The resulting effect of the randomly ordered phase encoding is a 
noise-like appearance of the motion artifact in contrast to the more 
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conspicuous appearance of image ghosts in the case of a sequential 
(standard) order of phase encoding [19]. 

In another approach called ROPE (respiratory ordering of phase 
encoding), the phase-encoding order is denned by the patient position 
[26]. As in the case of respiratory gating, patient displacements in 
ROPE are monitored by a special transducer attached to the patient. 
The collected information is used to acquire spatial frequencies in the 
phase direction, k y , at certain stages of the patient motion. In the case 
of periodical motion (e.g., in the phase-encoding direction), the displace- 
ment of the patient, y, becomes a periodical function of k y if the standard 
order of phase encoding is used. This periodicity gives rise to image 
ghosts in the phase-encoding direction. The objective of ROPE is to 
make y{k y ) a monotonic rather than a periodic function of k y by choos- 
ing the sequence of phase encoding according to the measured patient 
motion. Although it can be difficult to implement in practice, ROPE sig- 
nificantly reduces image ghosts without increasing scan time (unlike 
averaging and respiratory gating). 

Selective Saturation 

In many instances motion-induced artifacts can be significantly 
reduced by selectively saturating moving structures which can be 
omitted from the images. A slab of tissue (oblique, parallel or perpen- 
dicular to the imaging plane) that contains such structures, can be 
saturated by applying a series of frequency selective r.f. pulses in the 
presence of the slice-select gradient. Because the magnetization of 
spins in the saturated slab is small, the moving structures in the slab 
do not generate significant artifacts in images. It should be noted 
that selective saturation is frequently used to improve differentiation 
between arteries and veins in MR images by saturating a slab (normally 
adjacent to the imaged volume) that contains these blood vessels. Since 
the blood entering the image volume through the slab has a low magne- 
tization, it appears dark in images and causes minimal image distortion 
due to blood flow pulsations (Figure 6.10). 

Gradient Moment Nulling 

Motion induced phase shifts are linearly dependent on the 
moments of the gradient waveform: 



= I t m Gdt', 
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Figure 6.10. Schematic of selective saturation of flow. The arrows indi- 
cate the direction of flow. The flow entering the imaged volume from 
above has a low magnetization; therefore, it appears darker than the 
unsaturated flow entering the image volume from below. 

n = 1, 2, . . . The gradient moment nulling technique suppresses motion 
artifacts by using specially designed gradient waveforms with M„ = 
[27-29]. For example, if the gradient pulses have zero M x and M 2 , then 
the accumulated phase no longer depends on velocity or acceleration. 
Although gradient waveforms can be designed to null any chosen 
number of moments, nulling of higher-order moments requires a 
number of additional gradient lobes that significantly increase the 
minimum TE of a pulse sequence. As a result, gradient moment nulling 
is most often implemented by using gradient waveforms that only null 
the first-order moment, M x (Figure 6.11). 



Fractional Echo 

Since motion-dependent phase shifts occur only in the presence of 
magnetic field gradients, they can be reduced by shortening the 
duration of the gradient waveforms. To implement a pulse sequence 
with a compact readout gradient, Nishimura et al. [30] suggested using a 
bipolar gradient with a very short dephasing lobe applied prior to data 
acquisition. As a result, during the signal acquisition, only a fraction of 
the gradient echo is collected (Figure 6.12). Although fractional echo 
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Figure 6.11. Examples of gradient waveforms with first order moment 
nulling used for slice selection (a and b) and readout (c and d). 



acquisition tends to decrease SNR due to decreased readout time, the 
reduced dephasing of spins can actually increase the resultant SNR. 

In the fractional-echo approach, the acquired spatial frequencies 
are used to synthesize the missing part of the data. The reconstruction 
algorithms can be understood by considering an ideal case when image 
intensity is described by a real function. In this case, the &-space data, 
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Figure 6.12. Standard (a) and fractional echo (b) acquisitions : 
gradient-echo imaging. 
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S, are Hermitian: 

S(-k x , -k y ) = S*(k x , ky), (6.4.1) 

where S*(k x , k y ) is the complex conjugate of S(k x , k y ). From Eq. (6.4.1), 
it follows that acquisition of only half of spatial frequencies in the 
readout direction makes it possible to synthesize the unsampled half 
of the data. In practice magnetic field inhomogeneities, motion-induced 
phase shifts and other factors give rise to complex image intensity 
characterized by a spatially varying phase. Image reconstruction from 
partially acquired data in the latter case is discussed in Chapter 7. 

Note that compared with conventional gradient-echo acquisition, 
fractional-echo acquisition can be implemented with shorter echo 
time. This makes fractional-echo acquisition a useful technique for 
rapid gradient-echo imaging performed with very short TE and TR 
(see Chapter 7). 

Navigator Echoes 

The navigator approach is based on the assumption of a global 
linear motion of the imaged object. In this approach, object motion is 
recorded by acquiring special "navigator" echoes [31,32]. Motion- 
induced phase shifts of these echoes are subsequently used for phase 
correction of the acquired signals. 

To describe this approach we consider, for simplicity, the case of 
one-dimensional motion in the readout direction. In this case the 
acquired signal in 2D imaging can be expressed as 

S(k x ,k y (m)) = ff M xy (x - Ax(m),y)e jk * x+jk > {m)y dxdy. (6.4.2) 

In this equation the readout and phase-encoding gradients are in the x- 
and ^-directions, respectively; Ax(m) is the displacement at the time of 
a m th phase encoding. Equation (6.4.2) can be written as 

S(k x ,k y (m)) = e ik ** x{m) ff M X y(x,y)e jk * x+jk > {m)y dxdy 

= e jk '^ {m) S (k x ,k y (m)), (6.4.3) 

where S (k x ,k y (m)) is the signal from the same object at rest. From 
this equation it follows that the signal from a linearly moving object 
equals the signal from the object at rest multiplied by a phase factor 
e Jk x Ax(m)^ T k e ph ase f ac tor can be determined from an additional "navi- 
gator" echo acquired before applying the phase-encoding gradient 
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Figure 6.13. Gradient-echo pulse sequence with acquisition of a naviga- 
tor echo in the readout direction. 

(Figure 6.13). The navigator signal, S nav ^ m (k x ), acquired during the m th 
phase encoding is given by 

Snav,m( k x) = e jk ^ x(m) S navfi {k x ), (6.4.4) 

where S naU fi(k x ) is the reference navigator signal (e.g., acquired before 
the first phase encoding). By using navigator echoes, the signal from 
the object at rest can be restored as follows: 

S (k x ,k y (m)) = ^ avfi{ ll x) S(k x ,k y (m)). (6.4.5) 

This equation demonstrates that navigator signals can be used to 
eliminate the phase factor caused by motion in the direction of the 
navigator's gradient. In general, three orthogonal navigator gradients 
are required for phase correction in the case of three-dimensional 
linear motion. 



REFERENCES 

[1] J.M.S. Hutchison, R.J. Sutherland, J.R. Mallard. "NMR imaging: 
image recovery under magnetic field with large nonuniformities," 
J. Phys. E: Sci. Instrum., 11, 217 (1978). 



Chapter 6 Image Artifacts 



[2] K.M. Ludeke, P. Roschmann, R. Tischler. "Susceptibility artifacts 

in NMR imaging," Magn. Reson. Imaging, 3, 329 (1985). 
[3] K. Sekihara, S. Matsui, H. Kohno. "NMR imaging for magnets 
with large nonuniformities," IEEE Trans. Med. Imag., MI-4, 193 
(1985). 
[4] F. Farzaneh, S.J. Riederer, N.J. Pelc. "Analysis of T2 limitations 
and off-resonance effects on spatial resolution and artifacts in 
echo-planar imaging," Magn. Reson. Med., 14, 123 (1990). 
[5] J. Frahm, K.D. Merboldt, W. Hanicke. "The effects of intravoxel 
dephasing and incomplete slice refocusing on susceptibility con- 
trast in gradient-echo MRI," J. Magn. Reson. B., 109, 234 (1995). 
[6] J.R. Reichenbax, R. Venkatesan, D.A. Yablonskiy, M.R. Thomp- 
son, S.Lai, E.M. Haacke. "Theory and applcation of static field 
inhomogeneity effects in gradient-echo imaging," JMRI, 7, 266 
(1997). 
[7] B. Romner, M. Olsson, B. Ljunggren, S. Holtas, H. Saveland, 
L. Brandt, B. Persson. "Magnetic resonance imaging and aneur- 
ysm clips," J. Neurosurg., 70, 426 (1989). 
[8] W.T. Dixon. "Simple proton spectroscopic imaging," Radiology, 

153, 189 (1984). 
[9] J. Szumowski, D.B. Plewes. "Separation of lipid and water 
MR imaging signals by chopper averaging in the time domain," 
Radiology, 165, 247 (1987). 

[10] G.H. Glover, E. Schneider. "Three-point Dixon technique for true 
water/fat decomposition with B inhomogeneity correction," 
Magn. Reson. Med., 18, 371 (1991). 

[11] A. Haase, J.Frahm, W. Hanicke, D. Matthaei. ' a H NMR chemical 
shift selective (CHESS) imaging," Phys. Med. Biol., 30, 341 (1985). 

[12] V. Sklenar, Z. Starcuk. "1-2-1 pulse train: a new effective method 
of selective excitation for proton NMR in water," J. Magn. Reson., 
50, 495 (1982). 

[13] P.J. Hore. "Solvent suppression in Fourier transform nuclear 
magnetic resonance," J. Magn. Reson., 55, 283 (1983). 

[14] C.H. Meyer, J.M. Pauly, A. Macovski, D.G. Nishimura. "Simulta- 
neous spatial and spectral selective excitation," Magn. Reson. 
Med., 15, 287 (1990). 

[15] F. Schick. "Simultaneous highly selective MR water and fat 
imaging using a simple new type of spectral-spatial excitation," 
Magn. Reson. Med., 40, 194 (1998). 

[16] J. Szumowski, J.H. Simon. "Proton chemical shift imaging," 
Magnetic Resonance Imaging. D.D. Start and W.G. Bradly, eds., 
Mosby-Year Book, 1992. 



References 



[17] C.L. Schultz, R.J. Alfidi, A.D. Nelson, S.Y. Kopiwoda, M.E. 
Clampitt. "The effect of motion on two-dimensional Fourier trans- 
formation magnetic resonance images," Radiology, 152, 117 (1984). 

[18] M.L. Wood, R.M. Henkelman. "MR image artifacts from periodic 
motion," Med. Phys., 12, 143 (1985). 

[19] W.H. Perman, P.R. Moran, R.A. Moran, M.A. Bernstein. "Arti- 
facts from pulsatile flow in MR imaging," J. Comput. Assist. 
Tomogr., 10, 473 (1986). 

[20] G.L. Naylor, D.N. Firmin, D.B. Longmore. "Blood flow imaging by 
cine magnetic resonance," J. Comp. Assist. Tomogr., 10, 715 (1986). 

[21] C.B. Ahn, S.Y. Lee, 0. Nalcioglu, Z.H. Cho. "The effects of random 
directional distributed flow in nuclear magnetic resonance 
imaging," Med. Phys., 14, 43 (1987). 

[22] R.R. Edelman, P.F. Hahn, R. Buxton, J. Wittenberg, J.T. Ferrucci, 
S. Saini, T.J. Brady. "Rapid MR imaging with suspended res- 
piration: clinical application in the liver," Radiology, 161, 125 
(1986). 

[23] V.M. Runge, J. A. Clanton, C.L. Partain, A.E. James. "Respiratory 
gating in magnetic resonance imaging at 0.5Tesla," Radiology, 
151, 521 (1984). 

[24] M.L. Wood, R.M. Henkelman. "Suppression of respiratory motion 
artifacts in magnetic resonance imaging," Med. Phys., 13, 794 
(1986). 

[25] G.H. Glover, N.J. Pelc. "A rapid-gated cine MRI technique," 
Magn. Reson. Annual, 299 (1988). 

[26] D.R. Bailes, D.J. Gilderable, G.M. Bydder, A.G. Collins, D.N. 
Firmin. "Respiratory ordering of phase encoding (ROPE): a 
method for reducing motion artifacts in MR imaging," J. Comp. 
Assist. Tomogr., 9, 835 (1985). 

[27] P.M. Pattany, J.J. Phillips, L.C. Chiu, J.D. Lipcamon, J.L. Duerk, 
J.M. McNally, S.N. Mohapatra. "Motion artifact suppression 
technique (MAST) for MR imaging," J. Comput. Assist. Tomogr., 
11, 369 (1987). 

[28] E.M. Haacke, GW. Lenz. "Improving image quality in the pre- 
sence of motion by using rephasing gradients," AJR, 148, 1251 
(1987). 

[29] J.G. Pipe, T.L. Chenevert. "A progressive gradient moment nul- 
ling design technique," Magn. Reson. Med., 19, 175 (1991). 

[30] D.G. Nishimura, A. Macovski, J.I. Jackson, R.S. Hu, C.A. Stevick, 
L. Axel. "Magnetic resonance angiography by selective inversion 
recovery imaging using compact gradient-echo sequence," Magn. 
Reson. Med., 8, 96 (1988). 



Chapter 6 Image Artifacts 



[31] R.L. Ehman, J.P. Felmlee. "Adaptive technique for high-definition 
MR imaging of moving structures," Radiology, 173, 255 (1989). 

[32] H.W. Korin, F. Farzaneh, R.C. Wright, S.J. Riederer. "Compensa- 
tion for effects of linear motion in MR imaging," Magn. Reson. 
Med., 12, 99 (1989). 



CHAPTER 7 



Rapid MR Imaging 



The need to reduce motion artifacts and overall scan time in diag- 
nostic MRI provided the initial impetus for the development of rapid 
MR imaging. Rapid MRI with scan times of 0.1-1 sec per image has 
found a number of important clinical applications such as functional, 
perfusion, and cardiac imaging. Because of its vast potential and 
already proven importance for biomedical imaging, rapid MR scanning 
has been an active area of MRI research. 

A reduction in MRI scan time is commonly achieved by using two 
different approaches individually or in combination. The underlying 
strategy in the first approach is to shorten the time needed to collect 
NMR signals required for image reconstruction. This can be achieved 
by using a number of developed techniques such as gradient-echo 
imaging with short repetition times [1-3], echo-planar imaging [4,5], 
fast spin-echo imaging [6], and spiral imaging [7,8]. These imaging tech- 
niques will be discussed in more detail in the following sections of this 
chapter. In the second approach, known as partial k-space acquisition, 
a reduction in scan time is achieved by collecting only a fraction of full 
k-space data required to implement the conventional reconstruction 
algorithm. Image reconstruction in this approach is performed by 
using additional assumptions about the properties of the imaged 
object. Two of the partial k-space acquisition techniques will be dis- 
cussed in the last section of this chapter. 

7.1. RAPID GRADIENT-ECHO IMAGING 

Scan time in conventional MR imaging is given by the product of 
repetition time, TR, total number of phase encodings, and number of 
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times the signal is averaged, NSA. Assuming that NSA = 1, only two 
options exist for reducing scan time in the conventional approach: 
decreasing the number of phase encodings and shortening TR. Because 
in practice the minimum number of phase encodings is limited by the 
required spatial resolution, short scan times in rapid MR imaging 
with standard frequency and phase encodings are achieved by using 
very short 77? (normally on the order of 10 msec). 

Techniques for rapid gradient-echo imaging can be divided into 
two different groups. The first group is composed of pulse sequences 
that employ r.f. and gradient spoiling to destroy the transverse magne- 
tization after each signal acquisition while maintaining the steady- 
state longitudinal magnetization. The second group consists of pulse 
sequences with gradient refocusing used to sustain both transverse 
and longitudinal components of the magnetization in equilibrium. 

"Spoiled" Gradient-Echo Imaging 

Successful implementation of conventional MR imaging requires 
the existence of a steady-state in which the same magnetization in the 
object is created after each r.f. excitation. It has been shown [9-11] 
that a steady-state can be established by a train of repetitive r.f. 
pulses even when the local magnetic field is nonuniform, or when 
identical external gradients are applied between successive excita- 
tions. On the other hand, in conventional MR imaging the phase- 
encoding gradient is varied from one phase-encoding cycle to the 
next. It is therefore clear that the phase of the transverse magnetization 
also changes from cycle to cycle. The resultant dynamics of the magneti- 
zation during imaging critically depends on the interval between 
successive excitations. If TR is much longer than T 2 , then the trans- 
verse magnetization is negligibly small prior to each excitation. In 
this case a steady-state exists despite the fact that the phase of the 
transverse magnetization varies from excitation to excitation. On the 
other hand, when TR is smaller than T 2 the variations in the phase of 
magnetization do not allow the existence of a steady-state. 

To reduce nonequilibrium effects, the transverse magnetization in 
the sample can be destroyed (spoiled) by additional external gradients 
applied after every signal acquisition. These gradients (also known as 
crusher gradients) cause dephasing of spins, which in turn reduces 
the average transverse magnetization in a voxel. The most effective 
spoiling of the transverse magnetization is achieved when in addition 
to gradient spoiling the phase of r.f. pulses is varied from excitation 
to excitation [11,12]. A basic 2D gradient-echo sequence, known 
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Figure 7.1. 2D FLASH sequence with a spoiling gradient applied in the 
z-direction. 

as FLASH (fast low-angle shot) or spoiled GRASS (gradient 
recalled acquisition in the steady-state), with gradient "spoiling" is 
shown in Figure 7.1. An example of a FLASH image is shown in 
Figure 7.2. 

Gradient-echo image intensity in FLASH can be approximated as 

ZociY(ff)exp(-r.E/r2)sintt 1 ~ El , (7.1.1) 

1 — A 1 cosa 

where N(H) is the proton density, E 1 = exp(-Ti?/T 1 ), and a is the flip 
angle. When imaging with TR significantly longer than T 1 the longitu- 
dinal magnetization almost completely recovers during the interval 
between successive excitations. In this case 90 degree excitation 
pulses provide the greatest signal. Conversely, when repetition time 
is much shorter than 7\ the magnetization does not have enough 
time for recovery. In the latter case, flip angles smaller than 90 degrees 
should be used to maximize the signal (see Chapter 2). 

It should be noted that decreasing 77? results in SNR penalty. The 
maximum signal, achieve d when imaging with the Ernst flip angle (i.e., 
cosa = ^), varies as y/TR/l\ when 77? <C T v Therefore, we can 
conclude that at short TR the SNR is proportional to VTR. However, 
this reasoning does not take into account the fact that imaging with 
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Figure 7.2. 256 (readout) x 160 (phase encoding) FLASH image of a 
transverse section of the chest. Imaging parameters: field-of-view = 
34 cm x 34 cm, slice thickness = 7 mm, TR = 40 msec, TE = 2.2 msec, 
flip angle = 40 degrees. 

ultrashort TR (a few milliseconds) requires the use of very short read- 
out time, which further reduces the SNR. 

Equation (7.1.1) shows that in the case when TR < T|, image con- 
trast is denned by 7\ and proton density. The lack of T 2 contrast and 
low SNR are the main limitations of short TR spoiled gradient-echo 
imaging. To address the issue of T 2 -dependent contrast, Haase [3] has 
suggested using a modified FLASH pulse sequence in which nuclear 
magnetization in a sample is specially prepared prior to imaging. A 
magnetization-prepared pulse sequence used to acquire T 2 -weighted 
images contains an additional sequence of 90, 180, and 90 degree 
pulses (Figure 7.3). After the last 90 degree pulse flips magnetization 
back to the z-axis, the longitudinal magnetization is given by 

M 2 = M exp(-r/T 2 ), (7.1.2) 

where r is the interval between the first and the last preparation 
pulses. Crusher gradients are then used to destroy any residual trans- 
verse magnetization that exists due to incomplete refocusing of spins 
following the 180 degree pulse. After the preparation period, "snap- 
shot" FLASH imaging is performed with low flip angle (a < 5°) and 
short repetition times (a few milliseconds) in order to decrease effects 
caused by 7\ relaxation. Other modifications of FLASH allow acquisi- 
tion of inversion-recovery and diffusion-weighted images [3, 13]. FLASH 
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Figure 7.3. Magnetization-prepared snapshot FLASH sequence for 
acquiring T 2 -weighted images. 



imaging of the prepared magnetization is frequently referred to as MP- 
RAGE (magnetization prepared rapid acquisition gradient-echo). 



Gradient-Echo Imaging 

The idea of this approach, frequently referred to as FISP (fast 
imaging with steady-state precession) or SSFP (steady state free 
precession), is to create conditions under which both longitudinal 
and transverse magnetization are in equilibrium [14, 15]. The problem 
of accumulated phase changing from excitation to excitation due to the 
varying phase-encoding gradient is solved by applying a gradient 
(sometimes called rewinding gradient) with reverse polarity in the 
end of each repetition period. In principle, only refocusing along the 
phase-encoding direction is necessary because the slice-select and 
frequency-encoding gradients do not vary from cycle to cycle. How- 
ever, in order to null the accumulated phase rewinding gradients can 
be applied along all axes. 

Hinshow [16] has shown that in the case when the phase of r.f. pulses 
alternates from cycle to cycle, the steady-state magnetization can be rela- 
tively large even in the case of extremely short TR. Neglecting magnetic 
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field inhomogeneities, FISP image intensity obtained with alternate- 
phase excitations is given by 

ZocMoe-" ^ , F (1 " gl)SlnQ FF , (7-1.3) 

1-{E 1 - E 2 ) cos a - E X E 2 

where E 2 — exp(-TR/T 2 ) [16]. From this equation it follows that con- 
trast in FISP images acquired with long repetition time (i.e., TR > T 2 ) 
is similar to contrast in spoiled gradient-echo imaging. At very short 
repetition time (i.e., TR <C T 2 , Tj) we obtain from Eq. (7.1.3) 

T M sina ,„ .. 

I + T l /T 2 -(T l /T 2 - I) cos a K ' 

One important conclusion drawn from this equation is that the signal 
does not depend on repetition time. Therefore, FISP imaging with very 
short TR results in much smaller SNR penalty as compared to FLASH. 
Another important result is that image contrast in FISP depends on the 
ratio of T 2 /T 1 . For example, when imaging with a — 90° image inten- 
sity varies as (1 + T l /T 2 )~ l . This explains why tissues with large 
T 2 /T l ratios appear very bright in FISP images acquired with large 
flip angles. 

7.2. ECHO-PLANAR IMAGING 

Echo-planar imaging (EPI), originally suggested by Mansfield 
[4, 5], is one of the fastest MRI techniques currently used. A basic 
EPI pulse sequence is shown in Figure 7.4(a). Similar to conventional 
MR imaging, the signal acquisition in EPI is performed along parallel 
lines in k-space (Figure 7.4(b)). However, unlike conventional imaging 
EPI utilizes a very efficient sampling scheme by acquiring a number of 
lines in k-space after each excitation pulse. As a result EPI images can 
be produced with only a few excitations of the magnetization in the 
object. The fastest EPI scanning, known as single-shot EPI, is achieved 
when all required signals are collected after a single excitation pulse. 
Next we consider several EPI pulse sequences used in MR imaging. 

Gradient-Echo EPI (GRE-EPI) 

In GRE-EPI, a number (e.g., 128 or 256) of equidistant gradient 
echoes (echo train) with opposite polarities are produced by reversals 
of the frequency-encoding gradient (Figure 7.4(a)). The formation of 
an echo train in EPI can be understood by considering the phase of 
the transverse magnetization accumulated in the presence of the 
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Figure 7.4. (a) Gradient-echo EPI pulse sequence; (b) EPI k-space 
sampling (arrows indicate the direction of sampling). 

time-varying readout gradient. The peak of each gradient echo occurs 
in the center of the corresponding acquisition interval when the 
accumulated phase approaches zero (Figure 7.5). The echo-train time 
(ETT) is given by the product of the time interval between consecutive 
gradient echoes, referred to as echo spacing (ESP), and the number of 
gradient echoes generated after an excitation pulse. 

Each line in the GRE-EPI k-space corresponds to a phase-encoded 
gradient echo. Phase encoding in EPI is normally implemented by 
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Figure 7.5. Diagram of the readout gradient reversals and the 
corresponding phase reversals in EPI. 



applying "blips" of the phase-encoding gradient between reversals of 
the frequency-encoding gradient (Figure 7.4(a)). Each phase-encoding 
blip defines a new line in k-space. In contrast to the standard acqui- 
sition scheme, successive lines in the EPI k-space are sampled in 
opposite directions (Figure 7.4(b)). 

The major drawback of the highly efficient k-space sampling in 
GRE-EPI is frequent image artifacts caused by T 2 * decay or fat/water 
chemical shift. Image artifacts can be particularly severe in the case 
of single-shot EPI. Consider, for example, EPI with the following 
parameters: 256(readout) x 128(phase encoding) imaging matrix, 
256 kHz (±128 kHz) sampling bandwidth. Given these parameters 
each line in k-space is acquired during 256 x 10" 3 /256sec = 1msec. 
Assuming that the gradient ramp time (i.e., the time required to 
change the gradient from zero to its maximum amplitude) is much 
smaller than the acquisition time, we obtain ESP « 1 msec and 
ETT fa 128 msec. Because the calculated ETT for single-shot EPI is sig- 
nificantly greater than characteristic T| values in tissues (30-60 msec), 
gradient echoes at the tail of the echo train would be heavily T 2 *- 
weighted in contrast to gradient echoes in the beginning of the train. 
Suppose, for example, that phase encoding in EPI is implemented in 
such a way that low spatial frequency components of the signal are 
acquired at the beginning of the echo-train when the effect of T\ 
decay is relatively small, and high spatial frequencies are acquired 
at the tail of the echo-train when the signal significantly decreases 
due to T\ decay. The reduction in the amplitudes of high frequency 
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components amounts to low-pass filtering of the signal resulting in 
decreased spatial resolution and increased blurring in EPI images. 

Because of long readout time required to generate a number of 
gradient echoes after each excitation pulse, EPI images can display sig- 
nificant geometric distortion caused by magnetic field inhomogeneities 
and fat/water chemical shift. The potential severity of image distortion 
can be assessed by considering the effect of fat/water chemical shift 
in single-shot EPI. The phase difference between the water and fat 
magnetizations accumulated between consecutive gradient echoes is 
27rESP(i/ w - Vf), where v w and ly are the resonance frequencies of 
water and fat, respectively. The accumulated phase causes spatial 
shift (misplacement) of fat in the direction of the phase encoding 
rather than in the readout direction; the converse is true in the case 
of conventional imaging. The magnitude of the fat misplacement in 
single-shot EPI images is L 3 ,ESP(^ U) - u f ), where L y is the field of 
view in the phase-encoding direction. At a field strength of 1.5Tesla, 
the misplacement of fat in images acquired with 128 phase encodings 
and ESP of 1 msec is approximately 28 pixels. Because of such signifi- 
cant misplacement of fat, EPI is typically implemented with suppres- 
sion of the fat signal. 

Shortening of readout time is desirable because it decreases EPI 
image distortions. Readout time can be reduced by decreasing echo 
spacing, decreasing the number of phase encodings at the expense of 
spatial resolution and increasing the number of excitations (referred 
to as shots) at the expense of increased scan time. Shortening echo 
spacing can be achieved by either increasing sampling bandwidth, 
which requires strong gradients and causes a decrease in SNR, or by 
reducing the number of readout points at the expense of decreased 
spatial resolution in the readout direction. In practice, single-shot 
EPI with high spatial resolution is performed by using specially 
designed gradient systems providing both high gradient strength and 
short ramp time (about 100 microseconds). It should be noted that 
compared with rapid gradient-echo imaging considered earlier, 
single-shot GRE-EPI generally allows shorter scan time (i.e., 
~100msec per image) due to its higher sampling efficiency and higher 
SNR. These advantages of gradient-echo EPI have proven extremely 
useful for functional imaging of neuron activation in the brain (see 
Chapter 4). 

Because it uses a number of excitations to acquire image data, 
multi-shot EPI has shorter readout time compared with single-shot 
EPI. As a result, multi-shot EPI images display smaller distortion 
and loss of signal than single-shot EPI images. On the other hand, 
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increased scan time makes multi-shot EPI more susceptible to motion 
artifacts than single-shot EPI. 

Spin-Echo EPI 

Spin-echo EPI (SE-EPI) is implemented by using a refocusing 180° 
pulse after spatially selective excitation. A typical SE-EPI pulse- 
sequence is shown in Figure 7.6(a). By choosing the appropriate ESP 
and phase-encoding order, the center of k-space is normally sampled 
at the peak of the spin echo. Like in conventional spin-echo imaging, 
contrast in SE-EPI images depends on the echo time, TE, which defines 
T 2 -weighting of the center of k-space. Experience shows that due to 
refocusing of spins SE-EPI images are less sensitive to the presence 
of magnetic field inhomogeneities than GRE-EPI images. Compared 
with conventional spin-echo imaging, SE-EPI generally requires a sig- 
nificantly longer readout time. Consequently, SE-EPI is more suscepti- 
ble to artifacts caused by different T 2 weighting of the acquired signals 
and to artifacts caused by magnetic field inhomogeneities. 

Diffusion-Weighted EPI 

Diffusion-weighted EPI (DW-EPI) is performed with an additional 
pair of strong pulsed gradients applied on either side of a 180 degree 
pulse (Figure 7.6(b)). The attenuation of image intensity due to diffu- 
sion is given by the factor exp(-6D), where b depends on the gradient 
strength and duration, and the time interval between the diffusion 
gradients, and D is the diffusion coefficient [17]. Since EPI employs a 
significant number of imaging gradients, the equation for b should 
also include the contribution from these gradients. The diffusion coef- 
ficient can be calculated from a series of images acquired with different 
values of b (but the same 7\ and T 2 -weightings). Single-shot DW-EPI 
has found a number of applications in neuroimaging because DW- 
EPI is fast (~100msec per image), and therefore, almost free from 
motion-induced artifacts that often degrade conventional diffusion- 
weighted images. 

7.3. FAST SPIN-ECHO IMAGING 

Fast spin-echo imaging proposed by Hanning et al. [6] employs a 
number of 180 degree pulses to produce a train of spin echoes after 
each excitation of the transverse magnetization. To understand the 
process of formation of a spin-echo train, consider a simple model 
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Figure 7.6. (a) SE-EPI sequence; (b) DWI-EPI sequence. 
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Figure 7.7. Formation of a spin-echo train in the presence of 180 degree 
pulses. 



system composed of two noninteracting spins. Suppose that these spins 
have different Larmor frequencies u^ and w 2 . After a 90 degree excita- 
tion pulse, the phase of each spin will increase linearly with time: 
fa = ujit, fa = u> 2 t (Figure 7.7). A first 180° pulse, applied at time r 
after the excitation, inverts the phases of spins. Following the pulse 
a spin-echo is produced at time 2r when the phases fa and <fr 2 approach 
zero. Successive 180 degree pulses applied at t = r, 3r, 5r, It,..., 
(2N - 1)t cause a series of spin echoes at t = It, At, 6t, 8t, . . . , 2Nr 
(Figure 7.7). Spatial information can be encoded in the spin echoes 
by applying frequency and phase-encoding gradients between succes- 
sive 180 degree pulses. A basic fast spin-echo (FSE) pulse sequence, 
also referred to as RARE (rapid acquisition with relaxation enhance- 
ment) or turbo spin-echo, is shown in Figure 7.8. 

Like EPI, fast spin-echo imaging utilizes a very efficient sampling 
scheme by acquiring several lines in k-space after a single excitation. 
Because each line in the FSE k-space corresponds to a spin echo, 
FSE acquisition reduces artifacts due to magnetic field inhomo- 
geneities and fat/water chemical shift as compared to EPI. Because 
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Figure 7.8. Diagram of a basic FSE sequence. 

gradient echoes can be produced more rapidly than spin echoes, the 
echo spacing in EPI is usually several times shorter than the interval 
between successive spin-echoes (about 5-15 msec) during FSE acqui- 
sition. The relatively long echo spacing in FSE acquisition gives rise 
to image artifacts caused by different T 2 -weighting of the acquired 
signals [18, 19]. T 2 decay also limits the number of echoes that can be 
acquired during single-shot FSE imaging. For example, echo-train 
length in a single-shot FSE acquisition with 128 phase encodings and 
echo spacing of 10 msec would be 1280 msec. Since typical T 2 relaxation 
time in tissue is about 100 msec, the spin echoes at the tail of the echo 
train would produce very little signal. This example can explain 
why high resolution FSE images are usually obtained by using 
multishot acquisition of 4-16 successive spin echoes produced after 
each excitation of the magnetization in a sample. 

Contrast in fast spin-echo imaging is defined by the T 2 -weighting 
of central lines in k-space corresponding to small values of the 
phase-encoding gradient. A parameter defining T 2 -weighting of FSE 
images is the time interval, known as the effective echo time, between 
an excitation pulse and the center of a spin-echo with zero phase encod- 
ing. For example, heavily T 2 -weighted FSE images are acquired with 
effective echo time equal to or greater than T 2 in tissue. In contrast, 
reduced T 2 -weighting is achieved with effective echo time significantly 
shorter than T 2 (Figure 7.9). 



Chapter 7 Rapid MR Imaging 




Figure 7.9. Dependence of T 2 contrast in FSE imaging on effective echo 
time: (a) TE = 15 msec, (b) TE = 75 msec, (c) TE = 115 msec. FSE images 
of the brain were acquired with eight spin echoes per excitation using 
the following imaging parameters: matrix size = 256 (readout) x 
192 (phase encoding), field-of-view — 24 cm x 24 cm, slice thickness = 
2 mm, TR = 3 sec. 



FSE imaging has become widely used clinically because it pro- 
vides the required contrast between different tissues while being 
significantly faster than conventional spin-echo imaging. Contrast in 
FSE images can be enhanced further by using appropriately prepared 
magnetization. For example, in inversion-recovery (IR) FSE imaging 
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the magnetization in the object is inverted first by a 180° pulse. Follow- 
ing the inversion, the magnetization experiences Tj relaxation during 
a chosen delay after which FSE acquisition begins. Similar to con- 
ventional inversion-recovery imaging, IR FSE is frequently used to 
improve image contrast by nulling signals from certain tissues such 
as CSF. 

7.4. SPIRAL IMAGING 

Spiral imaging first demonstrated by Ahn et al. [7] is a very 
effective approach for rapid MRI. In contrast to EPI, spiral imaging 
utilizes smooth k-space trajectories that can be implemented on MR 
scanner without special hardware. Spiral acquisition reduces 
artifacts due to flow and motion as compared to conventional imaging 
with Cartesian k-space trajectories mainly because the center of k- 
space is acquired immediately after excitation with very little time 
for phase perturbations. Because of these attractive features, spiral 
imaging has been used for MR angiography, cardiac, and functional 
imaging. 

In the following discussion, we restrict ourselves to the case of 2D 
spiral trajectories that can be described by using polar coordinates 
(k,<p). The relationship between Cartesian (k x ,k y ) and polar coordi- 
nates is given by 

k x = k cos cf>, k y = k sin (f>. (7.4.1) 

Next we consider Archimedean spirals, which provide uniform radial 
sampling of k-space. The Archimedean spirals are described by the 
following equation [20]: 

k(t) = A<p(t), (7.4.2) 

where A is a constant. The radius of the acquired circular region in k- 
space, k max , is defined by the required spatial resolution, Ar. k max = 
ir/Ar. If a spiral trajectory reaches k max after N rev revolutions, then 
the maximum polar angle is (f> max = 2nN rev and from Eq. (7.4.2) it follows 
that 

A = 2N~Ar- < 7 - 4 " 3 > 

To avoid aliasing when imaging with Lx L field-of-view, the radial 
distance between acquired k-space points cannot exceed 2n/L. From 
this condition and Equations (7.4.2) and (7.4.3), it follows that spiral 
imaging with N x N matrix and M interleaves requires at least 
N/(2M) revolutions for each of the spiral trajectories. 
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The Cartesian components of the gradient, G, and components of 
the slew rate, dG/dt, are denned by the time derivatives of the vector 
k. That is, 

G x = 7 _1 /Jv = 7 _1 A(i(cos 6 - 6&a\d>), 

* ' (7.4.4) 

G y = ~f k y = 7 A(f){sm.(j) + 4>cos,(j)). 

The equation describing the linear "velocity" along a spiral trajectory 
in k-space can be written as 



t; = -f^Gl + Gy = A4>^fl + 4> 2 . (7.4.5) 

From Eq. (7.4.5) it follows that if v = const then at large values of t 

4>{t) « pVi, (7.4.6) 

where /3 is defined by the number of revolutions and acquisition (read- 
out) time, T s , required to reach k max : (3 = 2nN rev / \/T~ s . Spirals with 
constant linear velocity (Figure 7.10) are used frequently because 
they provide high SNR [8]. Notice that Eq. (7.4.6) cannot be used 
near the origin, where it requires infinite gradient strength and slew 
rate. In practice, gradient waveforms for spiral imaging are designed 
to provide desirable sampling near the origin of k-space while satisfy- 
ing the requirements imposed by the maximum slew rate and gradient 
strength available. 

To estimate readout time in spiral imaging with L x L field-of-view 
and bandwidth v s , we consider constant linear velocity spirals: 
v = 7G0 = const. In this case we have jG = 2itv s /L. If 27rN rev 3> 1, 
then we obtain the following estimate for the readout time, T s : 

ttAT 2 

For example, in the case of spiral imaging with 128 x 128 imaging 
matrix, 4 interleaves, and 128 kHz bandwidth, we obtain T s « 
25 msec. In comparison, the readout time in conventional imaging 
with the same matrix size and bandwidth will be only 1msec. As a 
result of their long readout time, spiral images display more blurring 
due to magnetic field nonuniformity than conventional images. 
Blurring in spiral images can be reduced by using image processing 
techniques [21, 22]. Likewise, it is reduced by increasing the number 
of interleaves thereby decreasing the readout time. The major limita- 
tion of the latter approach is that it prolongs scan time which is 
linearly dependent on the number of interleaves. 
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Figure 7.10. A constant linear velocity spiral: (a) trajectory in k-space; 
(b) time dependence of k x . 



Because k-space data collected on spiral trajectories do not lie 
on a Cartesian grid, spiral image reconstruction is more complicated 
than the reconstruction procedure in conventional imaging with 
Cartesian k-space trajectories. The acquired spiral data need to be 
interpolated first onto a Cartesian grid in order to subsequently 
implement computationally efficient Fast Fourier Transform. A com- 
monly used algorithm for reconstructing images from k-space data 
sampled along spiral trajectories is described in [23]. 
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7.5. PARTIAL K-SPACE ACQUISITION 

MR imaging with partial k-space acquisition allows significant 
reduction in scan time as compared to standard MR imaging by acquir- 
ing only a fraction of full k-space data. Two of the frequently used 
partial k-space acquisition schemes, Partial Fourier (PF) [24-26] and 
keyhole imaging [27], are discussed in the following sections. 

Partial Fourier Imaging 

From the properties of Fourier transform it follows that in the case 
when the image intensity is real the corresponding signal, S(k x , k y ), is 
Hermitian. That is, 

S{-k x ,-k y ) = S*{k x ,k y ), (7.5.1) 

where * denotes complex conjugation. In PF methods only a fraction 
of the full data set is acquired and the remainder of the data is syn- 
thesized through the use of Hermitian symmetry. PF imaging with 
acquisition of one half of k-space was employed by Feinberg et al. [24] 
to reduce the number of phase encodings without decreasing spatial 
resolution. The resulting penalty was a reduction in SNR by a factor 
of V2 because only half of the signals were actually acquired. If such 
decrease in SNR can be tolerated, this approach becomes an attractive 
option for rapid MRI because it reduces scan time by a factor of 2. 
Another application of PF methods is fractional echo acquisition, 
often used to shorten TE and readout time in gradient-echo imaging. 

In practice various factors (e.g., field inhomegeneities in the 
sample, eddy currents, etc.) cause the acquired data to be non- 
Hermitian. Several methods have been suggested in order to recon- 
struct MR images from partially acquired non-Hermitian data. The 
following discussion focuses on one of the earliest and most effective 
methods proposed by Margosian [25,26] (a review of this and related 
methods can be found in [28]). 

To describe the Margosian method we consider a one-dimensional 
complex image 

/(*) - \I{x)\e^ x \ (7.5.2) 

where \I(x)\ and </>(x) are the magnitude and phase of I(x), respectively. 
Suppose that Fourier transform of I(x) 

S(k) = f I(x) exp(;foc) dx (7.5.3) 
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is known for k > 0. Inverse Fourier transform of S(k) using only half 
the data defines a new function, I ha if(x). It can be shown that I half (x) 
can be expressed as follows: 

Ihaif(x)=lJs(k)eM-Jkx) d k = I -f + i- n j 1&L*>. (7.5.4) 

-oo 

If the image phase is known, then by multiplying Ih a if( x ) by e -7 ^ we 
obtain: 

Mx) \I(x)\ j 7 \I(*>)\eM^-W , 
haif(x) e J{ >= —£- + — J -r— x dx . (7.5.5) 



The Margosian method is based on the assumption that <j>{x) is a 
slowly varying function such that the difference \<j>{x) — 4>{x')\ is 
small and can be neglected. In this case the magnitude image can 
be obtained by displaying the real part of Ihaif( x ) e~ J ^ x ' in (7.5.5). To 
calculate slowly varying phase in 2D imaging with conventional 
rectilinear k-space trajectories we can use a low resolution image 
reconstructed from a small number of central lines in k-space. In 
summary, the basic reconstruction algorithm in the Margosian 
method includes: 

(a) Acquisition of about 60-75% of phase-encoded signals and zero fill- 
ing of all unsampled locations in k-space 

(b) Fourier transform of k-space data resulting in a complex image, 
i(x,y) 

(c) Reconstruction of a low resolution phase map 

(d) Calculation of the magnitude image, \I(x,y)\, by multiplying the 
complex image intensity, I(x,y), in each pixel by the factor e"^ 
and taking the real part of the product 

The Margosian method is frequently used for rapid MRI because it 
can easily be implemented with good results (Figure 7.11). However, 
this approach can sometimes cause image distortions for two main 
reasons. First, image phase can have significant contribution from 
high spatial frequency components that cannot be determined from a 
low resolution image. Second, the presence of large spatial variations 
in phase can violate the assumption that the difference \4>{x) - <f>{x')\ in 
(7.5.5) is small. 
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Figure 7.11. (a) Axial 256 x 256 spin-echo image of the head; (b) image 
of the same slice reconstructed from 65% of the phase encodings by 
using the Margosian method; (c) difference between (a) and (b). 

Dynamic Keyhole Imaging 

To describe this approach it is convenient to express image inten- 
sity, I, at some arbitrary time t as a sum: 



7(r,0 = /(r,« ) + A7(r ) *). 



(7.5.6) 
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Figure 7.12. Schematic of keyhole imaging: (a) reference data; (b) dyna- 
mically acquired data; (c) resultant data set produced by combining the 
outlined portions of the reference data set with the dynamically 
acquired data. 



In this equation r indicates spatial location in the sample, J(r, t ) 
denotes the reference image intensity, and AI(r,i) denotes the 
temporal change in intensity. The keyhole method is based on the 
assumption that A/(r, t) is primarily composed of slow varying spatial 
components. Clearly, under this assumption AI(r,i) can be recon- 
structed accurately with low spatial resolution. As a result dynamic 
keyhole imaging frequently acquires only low spatial frequency 
components of the signal. All unsampled locations in k-space even- 
tually are filled by using the data from the reference images. In 
practice, keyhole images are reconstructed from dynamically collected 
central k-space lines corresponding to small values of the phase- 
encoding gradient and peripheral lines from the reference images 
(Figure 7.12). 

Due to a reduced number of phase encodings, dynamic keyhole 
imaging allows significantly increased temporal resolution at the 
expense of reduced spatial resolution. In general, the keyhole 
approach is adequate if changes in intensity, AI(r,t), vary slowly in 
the field-of-view. Conversely, any sharp spatial variations in intensity 
appearing as a result of perturbations (e.g., due to motion or adminis- 
tered contrast agents) of the specimen will be blurred and will cause 
truncation artifacts due to insufficient resolution in the phase-encod- 
ing direction (Figure 7.13). 
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Figure 7.13. (a) Axial 256 x 256 spin-echo image of the head; (b) image 
acquired with the same parameters as in (a) after an intervention 
was simulated by inserting fluid-filled tubing (indicated by an arrow) 
into the field of view; (c) keyhole image reconstructed from the 
reference image in (a) and 64 phase encodings acquired after the 
intervention. Note that the tubing in the keyhole image is blurred 
because of insufficient spatial resolution in the phase-encoding 
direction. 
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CHAPTER 8 



MR Imaging of Flow 



The effect of flow on the observed NMR signal was first described 
by Suryan [1] soon after the discovery of NMR. Suryan observed an 
increased NMR signal from a flowing liquid as compared to the 
signal from the same liquid at rest. Suryan correctly explained the 
observed effect by the decreased saturation of the magnetization in 
the flowing liquid. A decade later the same effect was reported by 
Singer [2, 3] who, along with Bowman and Kudracev [4], developed 
several techniques for quantitative NMR flow measurements. The 
studies of these and other investigators laid the groundwork for MR 
flow imaging. Reference [5] contains an interesting review of early 
NMR studies of flow. 

After flow-related effects in MR images were first discussed by 
Young et al. [6] and Crooks et al. [7] in the early 1980s, various techni- 
ques have been developed for MRI of flow. In general, these techniques 
can be divided into three major categories based on the effects used for 
flow imaging: (a) "time-of-flight" effects [8-21]; (b) velocity-induced 
phase [22, 23]; and (c) signal enhancement caused by MR contrast 
agents [24-30]. 



8.1. TIME-OF-FLIGHT TECHNIQUES 

Gradient-Echo Imaging 

Spatially selective r.f. pulses are used in 2D and 3D MR imaging to 
repeatedly excite magnetization in a slice of material. When flowing 
spins enter an excited slice during imaging, they normally have a 
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higher longitudinal magnetization than that of the saturated static 
material in the slice. Therefore, in-flow of unsaturated spins makes 
the flowing material in images appear brighter than the static 
material. Time-of-flight (TOF) MR imaging utilizes the difference in 
magnetization of stationary and flowing spins in order to improve 
visualization of flow in MR images. Gradient-echo TOF MRI based 
on the in-flow effect is performed with short repetition time 
TR <C Ti, needed for saturation of the static material in the excited 
slice. Flow-induced dephasing of spins is reduced by imaging with 
short TE (a few milliseconds) combined with fractional-echo acqui- 
sition or by employing gradient moment nulling techniques (see 
Chapter 6). 

In 2D TOF imaging the signals are acquired from consecutively 
excited thin slices (1.5-3 mm thickness), which are typically chosen 
perpendicular to the direction of flow in order to minimize saturation 
of the flowing material (Figure 8.1(a)). Because it is often difficult to 
assess flow conditions from conventional two-dimensional MR 
images, additional images of flow can be produced by using a post-pro- 
cessing technique, known as the maximum intensity projection (MIP). 
In this technique an imaged volume is projected onto a plane. The 
determined maximum voxel intensities along rays perpendicular to 
the projection plane are displayed subsequently as two-dimensional 
MIP images. The MIP algorithm is very effective for selectively 
displaying high-intensity regions in the imaged volume. MIP images 
of the carotid and vertebral arteries in the neck reconstructed from 
images of thin slices are shown in Figure 8.1(b) and (c). 

Since 3D imaging in general provides higher SNR and allows 
better resolution than 2D imaging, it is often chosen as the technique 
for diagnostic MR imaging of arteries, known as magnetic resonance 
angiography (MR A). Because in 3D imaging signals are acquired fol- 
lowing excitation of a relatively thick slice (slab), the flowing material 
(blood) is continuously saturated by successive excitation pulses as it 
propagates through the slab. This saturation effect is particularly sig- 
nificant under the conditions of slow flow or flow travelling parallel to 
the imaged slab. The saturation of blood is spatially dependent, vary- 
ing from its minimum near the entrance to the slab to its maximum 
near the exit (Figure 8.2). Depending on the flow velocity, TR, and 
flip angle, this effect can lead to noticeable variations in the intensity 
of blood in 3D TOF images. Several imaging techniques have recently 
been suggested in order to ameliorate contrast in 3D TOF MRA. One of 
these techniques employs r.f. excitation with a spatially varying flip 
angle that increases from its minimum value at the slab entrance to 
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Figure 8.1. (a) A gradient-echo 
image of a representative axial 
slice through the neck of a volun- 
teer. Imaging parameters: in-plane 
resolution = 0.7 x 0.8 mm 2 , slice 
thickness = 1.5 mm, TR = 23 msec, 
7^ = 5 msec, flip angle = 60°. (b) 
and (c) show MIP images of the 
left and right carotid/vertebral 
arteries, respectively, recon- 
structed from 80 images of thin 
slices (like the one shown in (a)). 




its maximum value at the slab exit [13]. Example of MIP images of the 
brain reconstructed from 3D TOF data acquired with varying flip angle 
are shown in Figure 8.3. In another approach saturation of flowing 
material in 3D TOF MRA is decreased by sequential imaging of over- 
lapping thin slabs [14]. 

In many instances the close proximity of veins to arteries might 
make it difficult to differentiate between the two in MR images. Visua- 
lization of arterial structures in TOF MRA can be improved by selec- 
tively saturating the magnetization of venous blood (see Chapter 6). 
The use of selective saturation in TOF MRA is facilitated by the fact 
that the direction of arterial blood flow is typically opposite to 
that of venous blood flow. In practice the magnetization of venous 
blood is reduced routinely by selectively saturating a slab of tissue 
through which the blood flows prior to entering the imaged volume 
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Figure 8.2. Simulated image demonstrating spatial variations in inten- 
sity of flow in 3D TOF imaging caused by saturation of the flowing 
material in the imaged slab. The saturation of the flow continuously 
increases as it propagates through the slab. 



a n 

Figure 8.3. (a) MIP image of cerebral vessels in the brain reconstructed 
from a data set obtained with 3D acquisition. Imaging parameters: 
spatial resolution = 0.6 x 0.8 x 0.8mm 3 , TR = 40msec, TE = 3msec, 
average flip angle = 30°. (b) MIP image reconstructed from a data set 
acquired with the same parameters as in (a). In addition, the pulse- 
sequence used in (b) included an off-resonance preparation pulse pro- 
viding saturation of the macromolecular protons (see Chapter 4). The 
resulting MTC enhances visualization of small vessels (like the one 
indicated by an arrow in (b)). 



Time-of-flight techniques 



(Figure 6.10). The resulting effect is the dark appearance of veins rela- 
tive to the bright appearance of arteries in TOF images. 

Spin-Eecho Imaging 

Initially excited blood remains in the imaging volume during the 
time t = d/v, where d is the slice thickness and v is the component of 
the flow velocity perpendicular to the slice. If the time interval TEJ2 
between excitation and the spatially selective 180 degree pulse in 
spin-echo imaging exceeds r, then only the magnetization of the 
static material is refocused. Therefore, the washout of flowing spins 
from the imaged slice, known as the out-flow effect, makes blood vessels 
appear darker than the surrounding tissue in the image. For example, 
spin-echo images acquired with slice thickness of 2 mm, TE of 20 ms 
would display a greatly reduced signal from flow (flow void) if the 
flow velocity exceeds 20 cm/sec. Imaging techniques based on the out- 
flow effect [15] are frequently referred to as "black blood" MRA 
because of the darker appearance of blood relative to that of stationary 
material in the MRA images. Conversely, techniques used to acquire 
images with bright appearance of blood are known as "bright blood" 
MRA. Unlike bright blood MRA, black blood MRA techniques gener- 
ate flow contrast by minimizing the signal from moving spins. As a 
result, these techniques also minimize artifacts due to blood flow 
pulsations, saturating of the blood magnetization and flow-induced 
dephasing of spins which are often present in bright blood MRA 
images. 

Inversion-Recovery Imaging 

To improve contrast between flowing and stationary spins, Nishi- 
mura et al. [16] suggested a subtraction technique that can completely 
eliminate signal from the static material. In this technique two images 
are acquired as follows: the first image is obtained following a 180 
degree pulse that selectively inverts magnetization of the spins in 
the slab and above; the second image is acquired after the magnetiza- 
tion in the entire volume is inverted (Figure 8.4). In both cases the 
same time delay, TI, is used between the inversion and excitation 
pulses. Image subtraction leads to cancellation of the signals from 
both moving and stationary spins that have experienced inversion 
pulses. The resulting image is characterized by a bright appearance 
of vessels with upstream flow entering the slab from below. The 
major advantage of this approach as compared to the gradient-echo 
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Figure 8.4. A diagram describing selective-inversion MRA with image 
subtraction: (a) first image is acquired following inversion of the mag- 
netization in the slab and above; (b) second image is acquired after 
magnetization in the entire volume is inverted; (c) image subtraction 
results in cancellation of the signal from the static material in the 
imaged slab. 

techniques described earlier is potentially better visualization of ves- 
sels with slow flow, provided that the inversion time TI is long 
enough so that complete replenishment of the flowing material in the 
slab takes place. The major limitations are long scan time, because 
two images are acquired, and increased sensitivity to patient motion 
during data acquisition, because image subtraction is used. 

In another approach selective inversion of the longitudinal mag- 
netization is used in order to null the signal from blood [17]. In this 
black blood MRA approach a nonselective 180 degree pulse first 
inverts the longitudinal magnetization in the entire volume. Immedi- 
ately after the first inversion a spatially-selective 180 degree pulse 
returns the magnetization to its equilibrium state in the imaged slice 
only (Figure 8.5). During the inversion time, TI, the blood that initially 
experienced selective 180 degree pulse leaves the slice. The longitudi- 
nal magnetization of the blood entering the slice at time TI is propor- 
tional to 1 - 2e~ T1 l TlM °° d , where T lb i ood is the spin-lattice relaxation 
time of blood. The flow signal is nulled by choosing the appropriate 
inversion time: TI = T 1Mood ln2. 

Flow Imaging with Spatial Modulation of Magnetization 

Visualization of flow-related displacements of spins in MR images 
can be achieved by using specially prepared magnetization [18, 19]. In 
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Figure 8.5. Pulse-sequence for selective-inversion black blood MR 
imaging. 

this approach, known as spin tagging, the longitudinal magnetization, 
M z , in the sample is modulated initially by an arbitrary periodical 
function. Subsequent MR imaging depicts changes in the modulation 
pattern caused by flow or motion. 

Immediately after spin tagging the longitudinal magnetization can 
be expressed as 

M z (t, t = 0)= M T ag {r), (8.1.1) 

where M is the steady-state longitudinal magnetization and T ag (r) is 
an arbitrary periodical function varying between —1 and 1. Displace- 
ments of spins due to flow with velocity v shift the initial tagging 
pattern by vt. The resultant effect of flow and 7\ relaxation on the long- 
itudinal magnetization is described by the following equation: 

M z (r, t) = M e- (/T ' [T ag (r - vt) - 1] + M . (8.1.2) 

Periodical modulation of M z can be produced by using several 

spin-tagging techniques [18-20]. A basic spin-tagging pulse sequence 

includes a pair of nonselective 90 degree pulses (Figure 8.6). After 

the first 90 degree pulse the phase of the transverse magnetization, </>, 
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Figure 8.6. Basic pulse sequence for spin-tagging imaging of flow. 

changes under the influence of the applied gradient, G x : 4> = 7X / G x dt. 
Under the assumption that the relaxation and flow effects during spin- 
tagging can be neglected, the longitudinal magnetization immediately 
after the second 90 degree pulse is given by 

M Z = M cos fee, (8.1.3) 

where k = 7 J ep G x dt. Any residual transverse magnetization is sub- 
sequently destroyed by crusher gradients applied prior to imaging. If 
two-dimensional spatial modulation of the magnetization is needed, 
it can be produced by applying a pair of orthogonal gradients (e.g., 
G x and G y ) during the time interval between the 90 degree pulses. 

Following spin tagging the modulated nuclear magnetization is 
imaged. By using equations (8.1.2) and (8.1.3) we obtain that the trans- 
verse magnetization in the sample during imaging can be written as 

M tr oc M e' T/Tl cos k(x - v x t) + M„(l - e~ T/T '), (8.1.4) 

where v x is the x-component of the flow velocity and r is the interval 
between the spin-tagging and imaging parts of the pulse sequence. 
The spatially varying transverse magnetization appears as a series of 
bright bands in images (Figure 8.7(a)). In principle, local displacements 
of spins during the interval r can be assessed by examining distortions 
(bending) of the initially straight bands (Figure 8.7(b)). 
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Figure 8.7. (a) Axial 256 x 256 spin-echo image of a glass container 
filled with poppy seeds. The image was acquired after a delay of 
120 msec following modulation of the magnetization in the container. 
The distance between bright bands in the image is approximately 
3 mm. (b) Image of the container acquired with the same parameters 
as in (a) but after the flow of seeds in the container was induced by 
periodical shaking. Reprinted with permission from [21]. 



8.2. PHASE-CONTRAST MRA 

Phase-contrast imaging of flow utilizes the phase of magnetization 
in order to differentiate between moving and static materials [22, 23]. 
In this approach phase encoding of flow is implemented by using bipolar 
magnetic field gradients. A bipolar gradient, G x , can be defined as 

-G v , t <t<t + T 

G x = { G v , t + T<t<t + T + r (8.2.1) 

0, elsewhere 

where G v , r and T are the amplitude, duration, and separation of the 
two gradient lobes, respectively (Figure 8.8). 

Assuming that the flow velocity, v x , in the direction of the applied 
bipolar gradient is constant, the accumulated phase shift can be 
written as 



h + T + 



-yG x (t)x{t)dt = -yv x G v TT. 



(8.2.2) 
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Figure 8.8. A bipolar gradient used for phase-contrast imaging of flow. 

This equation shows that the flow-induced phase is directly propor- 
tional to velocity. Consequently, the overall phase shift for stationary 
material (v x = 0) is zero. 

In phase-contrast MRA, bipolar gradients are incorporated into 
MRI pulse sequences. The velocity-dependent phase can be used for 
flow imaging if the pulse sequence is repeated with the inverted bipolar 
gradient (Figure 8.9) and the signal from static material is cancelled 
out by calculating the difference between the reconstructed images: 



A7 = I„e''* - I v e 



= 2jl v sin <f>. 



(8.2.3) 



In this equation I v is a proportionality factor, which is linearly depen- 
dent on the number of moving spins in a voxel. MR images of flow can 



Figure 8.9. Gradient-echo pulse sequence for phase-contrast MRA. 



Phase-contrast MRA 



be produced by displaying the magnitude of AI. Because of the non- 
linear relationship between AI and <f>, bipolar gradients need to be 
experimentally adjusted to insure that <\> < 1. Under this condition 
sin (f> « and AI is approximately proportional to the component of 
flow in the direction of the applied bipolar gradient. 

In many instances it might be desirable to obtain images in which 
local image intensity is proportional to the magnitude of flow rather 
than components of flow. A measure of the flow magnitude can 
be obtained if the reconstructed images of orthogonal x, y and z- 
components of flow (i.e., AI X , AI y and AI Z , respectively) are combined 
as follows: AI = ^\AI X \ 2 + \AI y \ 2 + \AI Z \ 2 . The limitation of this 
approach is a significant increase in scan time because six images 
must be acquired before an image of the flow magnitude can be 
produced. 

The phase-contrast technique for flow imaging described above 
is well suited to imaging uniform flow. However, it cannot be used to 
measure the velocity distribution function of spins. This limitation 
can be overcome by using a general method for Fourier velocity 
encoding described by Moran [22]. In this method the strength of the 
applied bipolar gradient is incremented N times in a step-like fashion: 
G v = nG , -N/2 <n< N/2. Consequently, the velocity-dependent 
phase of the magnetization becomes dependent on n. For example, 
when imaging the x-component of flow the velocity-dependent phase 
is given by 

4>(n)=>yG v TTv x ='^, (8.2.4) 

where V g = tt/^GqtT. The NMR signal generated by spins with velo- 
city distribution function i(>(v x ) can then be written as 



S(n)« 



J i>(v x )e im ' x/V8 dv x . (8.2.5) 



According to this equation the signal is proportional to the Fourier 
transform of ip(v x ). Under the assumption that V g exceeds the 
maximum velocity of spins, the velocity distribution function can be 
reconstructed with resolution of 2 V g /N by computing inverse Fourier 
transform of S(n). It should be noted that this velocity encoding 
procedure results in a substantial increase in scan time (i.e., by a 
factor of N) when used in combination with techniques for spatial 
encoding. This practical limitation often makes it difficult to utilize 
the vast potential of the Fourier velocity encoding for MR imaging 
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8.3. CONTRAST-ENHANCED MRI OF FLOW 

Recently, a number of studies have demonstrated that contrast 
between static tissue and blood in MR angiography can be improved 
by using gadolinium-containing contrast agents that significantly 
reduce the spin-lattice relaxation time of blood [24-27]. The shortening 
of T x results in decreased saturation of the blood magnetization, lead- 
ing to improved visualization of blood vessels in MR images. Because 
this technique, also known as contrast-enhanced MRA, relies on the 
effect of Ti shortening it is less flow dependent than TOF MRA. As a 
result, contrast-enhanced MRA mostly avoids artifacts due to satura- 
tion of blood magnetization, which often degrade bright blood MRA 



Currently, contrast-enhanced MRI of flow is used primarily for 
imaging of the arterial system (Figure 8.10). Practical implementation 
of contrast-enhanced MRA is based on the use of fast 3D gradient-echo 
imaging with the shortest possible TR. Short scan time in contrast- 
enhanced MRA is needed to: (a) acquire signals during the so-called 
arterial phase (~10sec) when the contrast media exists primarily in 
the imaged arteries but not in the adjacent veins; (b) minimize the 
effect of contrast media leakage from the blood into the surrounding 
tissue; (c) enable breathholding which tends to minimize image distor- 
tion due to respiratory motion. Because the arterial phase is usually 




Figure 8.10. MIP image of the aorta pro- 
duced from a 3D data set acquired 
during 34 sec after administration of 
Gd-DTPA. Imaging parameters: spatial 
resolution = 1.5 x 1.75 x 1.3 mm 3 , TR = 
7 msec, TE = 2 msec, flip angle = 45°. 



shorter than the typical scan time in 3D MRA, it is a common practice 
to acquire only signal components with low spatial frequencies (i.e., 
central portion of k-space) during the arterial phase while components 
with high spatial frequencies are acquired after contrast media has 
started to accumulate in the veins and soft tissues. Acquiring the cen- 
tral portion of k-space too early (i.e., before the arrival of the bolus of 
contrast media) or too late (i.e., after the subsequent enhancement of 
veins and soft tissue) has the adverse effect of decreased contrast 
between the arteries and their surroundings. To achieve the optimum 
image contrast, acquisition of components with low spatial frequencies 
should start immediately after the bolus of contrast media has arrived 
in the imaged arterial structures. Two approaches have recently been 
suggested in order to determine the arrival time of administered con- 
trast material [28-30]. In the first approach, the time course of the 
vascular signal is monitored after administration of a small test 
bolus of contrast media. The determined arrival time is then used as 
a time delay between the subsequent administration of the full bolus 
and the beginning of the signal acquisition [28]. An alternative 
approach is based on monitoring the signal from a tracking volume 
that can include or be adjacent to the imaged site [29, 30]. The MRA 
signal acquisition can be triggered either manually when contrast arri- 
val is detected in the images of the tracking volume [29] or automati- 
cally when the signal from the tracking volume exceeds a chosen 
threshold [30]. 
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CHAPTER 9 



MRI Instrumentation: 
Magnets, Gradient Coils, and 
Radiofrequency Coils 



MR imaging requires three different types of magnetic fields. 
First, a static and highly uniform magnetic field, known as the 
main field, is required to create the initial longitudinal magneti- 
zation in the object and to maintain Larmor precession of nuclear 
spins at constant angular frequency. Second, an r.f. field is required 
for excitation of the transverse magnetization in the object. Finally, 
highly linear magnetic field gradients that can be turned on and 
off during a relatively short time (normally less than 1msec) are 
essential for spatial localization in MR imaging. Because of the dif- 
ferent functions and required characteristics of the main field, r.f. 
field, and gradient field, they are produced by different devices 
known as magnets, r.f. coils, and gradient coils, respectively. The 
purpose of this chapter is to describe main types and characteristics 
of these devices. 



9.1. MAGNETS 

The purpose of magnets used for MRI is to produce a static, 
strong, and homogeneous magnetic (main) field, B . Before discuss- 
ing the different types of magnets later in this section, it is instruc- 
tive to first consider the required characteristics of the main field 
used for NMR imaging and spectroscopy. 
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Main Field Stability and Strength 

During MR imaging or spectroscopic studies, signals are acquired 
repeatedly over time following successive excitations of the transverse 
magnetization in a sample. As a result, temporal stability of the main 
magnetic field is extremely important. The cause of temporal varia- 
tions in B can vary depending upon the magnet type and design. 
For example, main field in a resistive magnet can fluctuate because 
of the fluctuations in the supplied power. Consequently, power supplies 
for resistive magnets must be extremely stable. Main field in a perma- 
nent magnet is extremely sensitive to temperature fluctuations, which 
makes it necessary to monitor and control the temperature of the 
magnet. Compared with permanent and resistive magnets, supercon- 
ducting magnets provide the greatest temporal stability in operation 
with temporal variations in B less than B Q x 10" 7 per hour. 

The majority of clinical whole-body MRI scanners operate at field 
strengths between 0.3Tesla and 2Tesla. More recently, the need to 
increase the sensitivity of functional MRI of the brain has stimulated 
interest in proton imaging at field strengths in the 3-8Tesla range. 
High field strengths offer the advantage of increased signal-to-noise 
ratio in both NMR imaging and spectroscopy. In addition for spectro- 
scopic applications, higher field strengths make it easier to resolve 
individual lines in NMR spectra because the difference between 
Larmor frequencies of the same nuclei in different chemical environ- 
ments (chemical shift) increases linearly with field strength. 

Main Field Homogeneity 

Magnetic field inhomogeneity is detrimental for the quality of MR 
images because of the associated geometric distortion and signal loss. 
Field inhomogeneity is particularly troublesome for NMR spectros- 
copy because it might not allow differentiation between closely located 
individual lines in the NMR spectra. The degree of magnetic field 
homogeneity is commonly defined in terms of the normalized deviation 
of the field in a given spherical volume: 

In this equation B min and B max are the minimum and maximum field 
strengths in the spherical volume, respectively, and B c is the field 
strength at the center of the sphere. Magnetic field homogeneity is 
normally expressed in parts per million (ppm), 6B/B C x 10 6 . In practice, 
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main field homogeneity for different types of magnets is typically on the 
order of 10 ppm for a 50 cm diameter spherical volume (dsv). Compared to 
permanent and resistive magnets, superconducting magnets offer the 
highest field homogeneity (about 5 ppm in a 50 cm dsv for commercially 
available systems). The required homogeneity of the main field can vary 
depending on the size of the specimen as well as the type and objectives 
of NMR studies. For example, in vivo NMR spectroscopy may require 
homogeneity better than 0.1 ppm for a 10 cm dsv [1]. 

Permanent Magnets 

The unique property of permanent magnet materials is the exis- 
tence of microscopic regions, known as magnetic domains, with non- 
zero magnetic moments. In their natural state the magnetic moments 
of different domains are randomly oriented; however, in the presence 
of an external magnetic field, the domain moments become aligned 
along the field. It is remarkable that the alignment of magnetic 
moments remains after the magnetizing field is removed, making it 
possible to create permanent magnets from these materials. 

Modern materials for permanent magnets include neodymium- 
iron-boron alloys (NdFeB), ceramic ferrite, alnico (aluminum-nickel- 
iron-cobalt) alloys and samarium-cobalt alloys. The performance of a 
given magnetic material depends on two main parameters: the strength 
of an external magnetic field required to demagnetize the material, 
known as the coercivity, and the maximum achievable product of 
magnetic field and magnetic induction, known as the energy product. 
NdFeB and SmCo materials offer the highest energy product and 
coercivity among different permanent magnet materials currently used. 

Permanent magnets for whole-body MR scanners have several 
important limitations. First, the maximum field strength currently fea- 
sible for permanent magnets is relatively low (0.2-0.3 Tesla). Second, 
permanent magnets are very sensitive to temperature fluctuations. 
For example, the temperature coefficient of magnetic field for NdFeB 
is -0.12%/°C [2]; therefore, a change in temperature by only 0.1°C 
can cause a 120 ppm change in magnetic field. Because of such high 
temperature sensitivity, practical use of permanent magnets requires 
a special system for temperature monitoring and regulation. Other lim- 
itations of permanent magnets include their heavy weight (which 
increases with field strength) 1 and high cost of permanent magnet 
materials. 

1 For example, a 0.3 Tesla permanent magnet described in [3] has a mass of 90,000 kg. 
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Figure 9.1. Schematic of an H-frame magnet. 

A popular design for whole-body permanent magnets is known as 
the H-frame. An H-frame magnet contains an iron frame that provides 
a return path for magnetic flux and supports a pair of permanent mag- 
nets (poles) generating a magnetic field in the gap where the patient is 
positioned (Figure 9.1). The main advantages of this design are easy 
patient access (compared to commonly used solenoidal electro- 
magnets) and small fringe field due to intrinsic shielding provided by 
the frame. 



Resistive Magnets 

Magnetic field in resistive magnets is produced by electric current 
in specially designed coils. The coils are made of aluminum or copper, 
materials with high electric conductivity. Two basic types of resistive 
magnets used for whole-body imaging are air-core magnets and iron- 
core magnets. The latter type employs ferromagnetic materials to 
increase the strength of B - 

Magnetic field strength in an air-core resistive magnet is propor- 
tional to the operating current. Typical operating current and electri- 
cal resistance for a 0.15Tesla air-core magnet are about 200 A and 
lohm, respectively [3]. When in use, the magnet would require a 
40 kilowatt power supply. As a result of the high power dissipating in 
a resistive magnet, it is necessary to use a special cooling system 
that absorbs the produced heat and ensures the required stability of 
the magnet temperature. The latter is extremely important because 



temperature fluctuations can alter the conductivity of the magnet coils 
and thereby cause significant changes in B . Because the power con- 
sumed by an air-core magnet is proportional to B , the requirement 
of a highly stable power supply and effective cooling system makes it 
very costly to utilize air-core magnets at field strengths above 
0.2Tesla. Resistive magnets operating at higher fields in the 0.3-0.5 
Tesla range are designed by using the iron-core configuration com- 
posed of resistive wires encompassing two iron pole bars and an iron 
frame that provides a return path for the magnetic flux. 

Superconducting Magnets 

In 1911, H. Kamerlingh-Onnes discovered that as temperature 
decreased to approximetely 4°K the resistance of mercury suddenly 
became zero. This phenomenon of zero resistance at low temperatures 
is known as superconductivity, and materials demonstrating this 
phenomenon are referred to as superconductors. The use of superconduc- 
tivity was inhibited for several decades by the fact that the utilized 
materials, known as Type I superconductors, lost their superconduc- 
tivity as applied magnetic field exceeded a relatively low threshold 
(typically smaller than 0.1 Tesla). The use of superconducting materials 
has increased dramatically since the discovery in the late 1950s of new 
compounds, known as Type II superconductors. The discovered materi- 
als can exist in the superconducting state at very high magnetic field 
strength (on the order of 10 Tesla). Because of this property, Type II 
superconductors have become the material of choice for high field mag- 
nets operating at B > 1 Tesla. All known superconductors exhibit the 
superconducting state only at temperatures below a critical tempera- 
ture, T c , which is characteristic of the material used. Since conventional 
Type I and Type II superconductors have very low T c (normally less than 
10-20°K), their use usually requires cryogenic systems that employ 
liquid helium as a coolant. New compounds, discovered in the 1980s 
and known as High Temperature Superconductors (HTS), have much 
higher T c , which makes it possible to use relatively inexpensive liquid 
nitrogen as a coolant instead of the more expensive liquid helium. 
However, HTS materials with both high mechanical stability and the 
ability to maintain the superconducting state at the current densities 
required for high field magnets have yet to be developed. 

Modern superconducting materials for MRI magnets include 
niobium-titanium (NbTi) and niobium-tin (Nb 3 Sn), with NbTi 
(T c = 9.5°K) being used most often. A schematic of a typical super- 
conducting magnet is shown in Figure 9.2. To maintain the 
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Figure 9.2. (a) A diagram of a superconducting magnet; (b) the main 
field B in a whole-body superconducting magnet is parallel to the 
long axis of the patient. 



temperature of the NbTi coils below T c they are placed in a tank 
with liquid helium. The heat flow to the liquid helium from its surround- 
ings is minimized through the use of two radiation shields mounted 
outside of the helium tank and kept at low temperatures by an auxiliary 
cryostat. In some magnets additional thermal isolation of the liquid 
helium is achieved by employing an outer tank filled with liquid 
nitrogen. 
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When activating a superconducting magnet, the current is created 
initially by an external power supply connected to the superconducting 
coils through a special switch. While the magnet is being energized, the 
switch is kept in the normal state by heating. After the operating 
magnetic field is achieved, the magnet is disconnected from the power 
supply and the heating of the switch stops. The switch then returns to 
its superconducting state required to establish a zero resistance path 
for the operating current in the magnet. 

The superconducting state can suddenly disappear as a result of 
internal heat generated in the magnet coils for a number of reasons. 
One source of heat may be friction that occurs during small displace- 
ments of the superconducting coils caused by Lorentz forces. The 
transformation to the normal state, known as quench, can be quite 
dramatic because of a rapid (a few seconds) dissipation of the magnetic 
field energy. The rapid dissipation of energy in a superconducting 
magnet causes boiling and evaporation of helium and may damage 
the magnet coils. Nowadays, as a result of the improved design, 
magnet quenching has become a rare event that occurs without 
causing permanent damage to the magnet. 

Design of Coil Arrangements for Electromagnets 

When designing electromagnets (resistive or superconducting), 
the objective is to devise coil arrangements that provide a specified 
magnetic field uniformity. To understand the basic approach for 
minimizing spatial variation in B , we consider arrangements of co- 
axial coils (Figure 9.3). Since in this case the main field B is oriented 
primarily along the coil axis (shown as the z-axis in Figure 9.3), we can 
neglect the components of B in the plane perpendicular to the axis. 
The equation describing the longitudinal component, B z , in the 
region free of current is 

AB Z = 0, (9.1.2) 

where 



dx 2 dy 2 dz 2 ' 

Using spherical coordinates r, <j>, and 6, we can express a general solu- 
tion of equation (9.1.2) in the following form: 

B z = const + Y2J2 rnp n( c °s 6 )i. A n cos m<p + B™ sin nup). (9.1.3) 
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Figure 9.3. A set of coaxial coils forming an electromagnet. 



In this equation P% are the associated Legendre polynomials, and A™ 
and B™ are arbitrary coefficients. Using the explicit dependence of 
the lower-order Legendre polynomials on (f> and 6 [4], we obtain that 
the linear field gradients in the x-, y-, and z-direction are described 
by the terms with A\, B\, and A\, respectively. The higher-order 
terms in equation (9.1.3) for which n > 1 describe any spatial 
variations in B z that are nonlinear as a function of distance from the 
origin. 

To minimize spatial variations in B z near the origin, we need to 
devise a configuration of coils for which all coefficients A™ and B™ 
with n = 1, . . . ,N - 1 are zero, where N indicates a designated order 
of field uniformity. A simple set of two identical coils placed symmetri- 
cally above and below the mid-plane of the magnet and carrying iden- 
tical parallel currents (Figure 9.3) can be used as a starting point for 
magnet design. Because of the axial symmetry and the symmetry 
with respect to the mid-plane of the two-coil magnet, all coefficients 
A™ and B™ with m > or odd n are zero and the field uniformity in 
the vicinity of the origin is defined by A\. The field uniformity can be 
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improved by employing the Helmholtz pair of coils for which the axial 
distance between coils equals the coil radius. In this case A\ is zero and 
the field uniformity is defined by A4. By adding two additional coils it is 
possible to null A®, leaving the term with A@ as the main contaminant 
for B z . The sixth-order contaminant can be eliminated in a system of 
six coils [1], which is used in practice because it provides sufficiently 
high field uniformity. In principle, by increasing the number of coils 
it is possible to null all spatially dependent contaminants of B z up to 
any order. 

Shimming 

Although magnets for MRI are manufactured with great preci- 
sion to satisfy the design requirements, including coil geometry, 
symmetric positioning of the coils, etc., small deviations from the 
specified parameters are practically unavoidable. As a result, the 
main magnetic field on a MRI scanner typically needs further adjust- 
ment, which is done at the imaging site by employing a process known 
as shimming. One approach for shimming, referred to as passive shim- 
ming, is based on the use of a number of pieces of ferromagnetic 
material (e.g., iron). Passive shimming normally requires computer 
analysis to determine the optimal locations and the number of 
ferromagnetic pieces needed to provide the required magnetic field 
uniformity. The second approach, known as active shimming, utilizes 
special shim coils (resistive or superconducting) to produce an auxili- 
ary magnetic field B shim , which cancels spatial variations in the main 
field B„. 

Another objective of shimming is to improve field homogeneity in 
the imaged specimen, which itself creates a nonuniform magnetic field. 
In the presence of a specimen, B s hi m is usually adjusted interactively 
based on the acquired signal. Interactive shimming can be performed 
manually or automatically (i.e., computer-based shimming) by varying 
the shim currents. After each variation, the signal from a chosen 
volume is acquired and examined. Initially, the signal rapidly decays 
due to significant magnetic field nonuniformity in the specimen. As 
the field uniformity improves, the signal decay becomes slower and 
the corresponding spectral linewidth decreases. The shimming stops 
after achieving the minimum linewidth. Alternatively, shimming can 
be implemented by first calculating the map of the field distribution 
in the sample and subsequently using a numerical algorithm to deter- 
mine the optimal shim currents [5-7]. Magnetic field mapping can be 
performed, for example, by collecting signals during consecutive 
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Figure 9.4. Gradient-echo pulse sequence for field mapping. 



gradient echoes (Figure 9.4). The phase difference 8(f> between the 
image intensities acquired with echo times TE l and TE 2 , respectively, 
is given by 



= 1 8B{TE X - TE 2 ) 



(9.1.4) 



where SB denotes local field inhomogeneity. After 8B is determined 
from the phase difference, the optimal currents for different shim 
coils can be calculated by minimizing the quantity 



' (6B + B shim Ydxdydz, 



(9.1.5) 



assuming that the magnetic field created by each shim current is 
known. 



Shielding 

Magnetic field outside a MR scanner (fringe field) can extend over 
a significant area and thereby interfere with other devices such as 
computers, video screens, data storage devices, etc. Consequently, 
the need to protect the sensitive equipment makes it necessary to mini- 
mize the fringe field by utilizing special magnetic screening (shielding) 
techniques. Passive shielding is one common approach used to limit 
propagation of the main field outside the magnet. It employs specially 
designed ferrous sheets placed between the magnet and the field- 
sensitive area. The second approach, referred to as active shielding, 
is based on the use of auxiliary coils. The purpose of the specially 
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designed auxiliary coils is to generate a magnetic field that counter- 
acts the main field produced by the magnet, thus reducing the resultant 
fringe field. 

9.2. GRADIENT COILS 

As discussed earlier, spatial localization in MRI is achieved 
through the use of three orthogonal magnetic field gradients. If the 
z direction is chosen along the main field B , then these gradients 
can be expressed as G x = dB z /dx, G y = dB z /dy and G z = dB z /dz. 
MR imaging generally requires a high degree of gradient linearity 
throughout the sample volume, rapid switching of the gradients and 
high gradient strength. An estimate of the minimum gradient strength 
required for two-dimensional Fourier imaging can be obtained by 
requiring that the readout time T s be smaller than T% relaxation 
time in each voxel. This condition is needed to avoid significant 
image distortion. The relationship between the readout gradient, G x , 
readout time and pixel width in the readout direction, Ax, can be 
written as 

27r - ^ (9.2.1) 

At a typical T\ of 50 msec and Ax of 1mm the minimum gradient 
required is approximately 0.05 gauss/cm. Much higher gradients are 
required for rapid MR imaging with very short T g . For example, the 
readout gradient in EPI performed with T s of 1 msec and the same spa- 
tial resolution as in the example above, is approximately 2.3 gauss/cm. 
Most commercially available whole-body MR scanners have a maxi- 
mum gradient strength between 10 and 50 mT/m (1-5 gauss/cm) and a 
switching time in the 0.1-1 msec range. 

In practice, gradient systems are composed of different coil 
arrangements used to create the x, y, and z gradients. A general math- 
ematical approach used for gradient design employs solutions of the 
Laplace equation, AB Z = 0, with boundary conditions determined by 
the applied currents, geometry and arrangement of the gradient 
coils. For the purpose of our discussion, however, we will employ a 
rather conceptual description of the subject by focusing on basic 
arrangements of gradient coils. 

A commonly used coil configuration for a z-gradient (longitudinal 
gradient) is the so-called Maxwell pair, which consists of two identical 
coils carrying currents of equal magnitude, /, in opposite directions. 
These coils are wound on a cylinder coaxial with the magnet coils 
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producing the main magnetic field. The longitudinal component of 
magnetic field created by a Maxwell pair on the axis of the magnet 
at distance z from the origin can be expressed as 



2 



- + z +a* 



(9.2.2) 



where a is the coil radius and b is the distance between the coils. By 
using the Tailor expansion for B z one can show that if the distance 
between the coils equals ay/S then the main contaminant of the 
linear field gradient near the origin is proportional to z 5 . Use of 
additional coils makes it possible to improve the gradient linearity 
by nulling the fifth-order term. However, adding more coils increases 
the gradient system inductance. An increase in inductance in turn 
prolongs the so-called rise time needed to increase the gradient 
strength from zero to its maximum value. As a result, design of gradient 
coils involves a compromise between the conflicting requirements of 
short rise time and a high degree of gradient linearity. 

Magnetic field gradients in the x- and y-directions (transverse 
gradients) are commonly produced by using two identical coil arrange- 
ments that are rotated by 90 degrees around the z-axis relative to each 
other. A frequently used design for transverse gradients, known as the 
Golay coils, is a double-saddle configuration of eight 120° circular arcs 
(Figure 9.5). The arcs are connected by eight straight segments parallel 
to the z-axis. Note that a transverse gradient (i.e., dB z /dx or dB z /dy) is 
created by the currents in the arcs without any contribution from the 
connecting segments. In practice additional arcs are often added to 




Figure 9.5. Double-saddle gradient coil. 
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improve the gradient linearity [2]. References [1,2,8] contain detailed 
considerations of the Golay coils as well as other coil arrangements 
for the transverse field gradients. 

The size of gradient coils is a very important parameter in gradient 
system design. Consider, for example, the z-gradient produced by a 
Maxwell pair. By using equation (9.2.2) it is easy to see that near the 
origin G z <x I /a 2 . This dependence implies that at fixed gradient 
strength the gradient current I is proportional to a 2 . Because the 
coil resistance is proportional to the coil radius, the voltage across 
the coil is proportional to a 3 and the power dissipation in the gradient 
system is proportional to a 5 . This indicates that a significant improve- 
ment in gradient performance can be achieved by reducing the size of 
gradient coils. 

Eddy Currents 

Rapid switching of magnetic field gradients during MR imaging cre- 
ates large changes in magnetic flux through conducting structures in 
the magnet (e.g., shim coils, main magnet coils, r.f. coils, etc.). These 
changes in magnetic flux give rise to eddy currents, which in turn pro- 
duce both temporally and spatially varying magnetic field. Through 
their magnetic field, eddy currents cause image distortion and loss of 
signal. The magnitude, duration, and spatial distribution of eddy cur- 
rents depend upon various factors including the following: gradient 
strength and rise time, geometry of the gradient coils, distance between 
the gradient coils and other conducting structures in the magnet. 

One approach for minimizing the effect of eddy currents entails the 
use of specially designed gradient waveforms that compensate for 
the magnetic field created by eddy currents [9,10]. To understand the 
basic idea of the electronic compensation, consider the time depen- 
dences of the ideal gradient waveform and the contaminant magnetic 
field produced by the eddy currents. Since the contaminant field 
tends to negate the gradient field at the beginning of the pulse, the 
gradient strength needs to be appropriately increased to provide the 
necessary overshoot at this stage of the pulse (Figure 9.6). Conversely, 
the gradient strength at the tail of the pulse needs to be appropriately 
decreased to compensate for the positive contribution of the eddy 
currents. In practice the induced eddy currents may require compen- 
sation for many milliseconds after the gradient field has changed. 
Consequently, electronic compensation of eddy currents can prevent 
the use of compact gradient waveforms making it difficult to achieve 
short echo and readout times in rapid MR imaging. 
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Figure 9.6. Schematic illustration of eddy-current compensation: (a) 
dashed and solid lines indicate the gradient waveforms with and with- 
out eddy-current compensation, respectively; (b) contaminant gradient 
produced by the eddy-currents; (c) dashed and solid lines indicate the 
resultant gradient with and without eddy-current compensation, 
respectively. 

A better solution to the eddy current problem is to shield the 
conducting parts of the magnet from the time-varying gradient field. 
In practice this approach entails the use of the so-called actively 
shielded gradient coils with minimum fringe field [11]. An actively 
shielded gradient coil actually consists of two coils: an inner (primary) 
coil, which produces the required gradient field, and an outer 
(shield) coil, which is designed to cancel the fringe field created by 
the inner coil. The mathematical framework for the design of actively 
shielded coils is discussed in [8,11]. 

9.3. RADIOFREQUENCY COILS 

The main functions of r.f. coils are to excite magnetization in a 
sample and to receive the signal produced by the excited magnetization. 
These two functions may be performed by using different transmitting 
and receiving r.f. coils. Alternatively, the same coil may be used both 
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as a transmitter and as a receiver. The requirements for r.f. coils vary 
depending on the objective of a particular MRI study. For example, 
when imaging the brain, it may be essential that a utilized coil produce 
a relatively uniform r.f. field in the imaging volume. On the other hand, 
if r.f. field uniformity is not critical, it may be desirable to use the so- 
called surface coils that are built in various sizes and shapes to 
accommodate the geometry of the sample. 

Coil Sensitivity 

From the principle of reciprocity [12], it follows that the induced 
emf in the receiver coil in general increases with an increasing ratio 
of the magnetic field B x produced by the coil to the coil current 7. 
When the ratio Bi/I, known as the coil sensitivity, is high, a given 
coil can also be effective as a transmitter. From the Biot and Savart 
law it follows that 

£-£/^ 

In this equation Idl is a current element, r is the vector between dl and 
the location where B x is measured. 2 Equation (9.3.1) shows that the 
ratio B x /7 depends only upon the coil geometry. If the z axis is taken 
in the direction of the main field B , we can express the sensitivity 
B x /7 as the sum of its transverse component, B ltr /7, and longitudinal 
component, k7? l2 /7, where k is the unit vector in the z direction. 
According to the principle of reciprocity, the induced emf in the coil 
can then be written as 

«*-- / ¥ ^-- / ¥%**> <•«> 

sample sample 

where M tr and kM z are the transverse and longitudinal magnetiza- 
tions, respectively. 

Following an excitation pulse the transverse magnetization 
exhibits oscillatory motion with frequency "fB . Since the inequality 
jB Tx » 1 is normally satisfied in NMR spectroscopy and imaging, 
the second term on the right side of equation (9.3.2) can be neglected. 
It is then clear that the coil sensitivity in the transverse plane (i.e., 
plane perpendicular to B ) is the quantity that should be maximized 
when designing r.f. coils, while the longitudinal component of Bj is 

2 Strictly speaking, equation (9.3.1) is valid at low frequencies, when the wavelength of 
Bj in air is large compared to the dimensions of the r.f. coil. 
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inconsequential. In particular, we can conclude that solenoidal coils 
are not practical for diagnostic MR imaging with superconducting 
magnets for which B is parallel to the long axis of the patient, because 
these coils produce mostly a longitudinal r.f. field when encompassing 
the patient. On the other hand, solenoidal coils are useful when 
imaging with permanent or resistive magnets that have a main field 
perpendicular to the long axis of the patient. 



Tuning and Matching 

The efficiency of a given NMR coil as a transmitter and receiver 
of electric power depends upon its inductance L, resistance R and 
capacitance C. In many instances NMR r.f. coils can be described by 
using equivalent LRC circuits (Figure 9.7a). Suppose that V e jut is 
the driving emf during an excitation pulse. The equation describing 
the current, I, through the circuit can be written as 

V e iut = L^- + RI + l, (9.3.3) 



dt 



C" 



where q is the charge stored by the capacitor. It ii 
the induced current is given by 

Ve H>*-*) 



y to v 
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Figure 9.7. (a) Circuit diagram for a typical r.f. coil; (b) impedance 
matching scheme. 
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tancf) = 

To effectively generate the required Bj field, it is generally desirable to 
have the maximum current amplitude at the frequency u> — ^B . In 
practice this can be achieved by adjusting (tuning) the coil capacitance 
in such a way that 

lC = <4- (9-3-5) 

This condition implies that the r.f. coil resonates at the Larmor fre- 
quency of the nuclei. The width of the resonance at half-height, 6w, 
can be expressed as 8uj « \/3cj /Q, where Q, known as the coil quality 
factor, is given by 

Q-^f- (9-3.6) 

As discussed previously (see Chapter 5) the signal-to-noise ratio is 
inversely proportional to the square root of R, and since a high Q indi- 
cates a low resistance R, it is desirable to design r.f. coils with high Q 
values. Note that because the effective coil resistance R is given by the 
sum of the sample and the coil resistances [13], the quality factor for a 
"loaded" coil (i.e., sample is present) will be significantly smaller than 
the quality factor measured without a sample if the coil resistance is 
quite small compared to the sample resistance. 

When operating, the NMR r.f. coil is connected with the trans- 
mitter and receiver electronics via a coaxial cable (transmission 
line). To transmit electric power to and from the coil most efficiently, 
the impedance of the r.f. coil must be equal to the impedance of the 
transmission line. The line impedance is typically 50 fi resistive. 
Because the r.f. coil resistance is usually less than 50 O, it is necessary 
to use impedance matching networks to ensure that the r.f. coil is 
equivalent to a 50 fl resistance when operating at the resonance fre- 
quency uj . The transformation of the coil impedance to 50 ft can be 
achieved through the use of various matching schemes. One such 
scheme employing external capacitors C x and C 2 is shown in Figure 
9.7b. It can be shown that for this scheme the required matching of 
the coil impedance occurs when 

**L 1 7r =[R(R -R)} 1/2 and -L- = R (R /R - 1)~ 1/2 . (9.3.7) 
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.8. Circular surface coil. 



Coil Types 

Among the various types of r.f. coils used for imaging, surface coils 
are the most basic. Surface coils offer the advantage of high sensitivity 
when used for imaging superficial structures in a specimen. An 
example of a surface coil is a circular loop of conducting wire 
(Figure 9.8). The B x field along the coil axis (taken as the x axis) is 
given by 



B,=- 



Mo-fa 



2(<r 



2^3/2 ' 



(9.3.8) 



where a is the coil radius. By examining this equation we can make sev- 
eral important conclusions. First, an r.f. field produced by a single sur- 
face coil is very inhomogeneous, making it difficult to obtain a good 
slice profile during MR imaging. As a result, surface coils are used 
primarily as receiving coils. Second, the sensitivity of a given surface 
coil rapidly decreases with increasing distance from the center of the 
coil. This explains why surface coils are used for imaging only when 
they can be placed in close proximity to the regions to be imaged. 
Third, according to equation (9.3.8) the sensitivity at the center of a 
circular surface coil is inversely proportional to the coil radius. There- 
fore, small surface coils placed on the surface of a sample provide a 
higher sensitivity near the surface compared to that of large coils. 
On the other hand, the sensitivity of small surface coils decreases 
more rapidly with increasing depth from the surface than the sensitiv- 
ity of large coils (Figure 9.9). 

When used for imaging, a single surface coil normally provides 
acceptable B 1 homogeneity only in the regions that have dimensions 
smaller than or comparable to the diameter of the coil. Consequently, 
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Figure 9.9. Sensitivity vs depth for circular surface coils with different 
radii: (a) a = 2 cm; (b) a = 4 cm; (c) a = 6 cm. 



there exist two options for imaging a large region with surface coils: (a) 
use of a single large coil; (b) use of an array of relatively small coils, 
known as the phased array coil [14]. A phased array coil consists of a 
number of noninteracting surface coils (e.g., 4 or 8) that simultaneously 
receive NMR signals. Although each surface coil effectively images only 
a small region within the field-of-view, individual images produced by an 
array of coils can be combined together to obtain an image of a large area 
with uniform sensitivity and with SNR higher than that of a single large 
surface coil. 

For many applications of MRI such as body or head imaging, a 
homogeneous transverse J3j field is essential. A transverse B x can be 
produced by a current flowing on the surface of an infinitely long 
cylinder parallel to its axis. In this case, the condition required for a 
uniform Sj is that the current density is proportional to sin^ (or 
coscp), where (p is the azimuthal angle [15]. The sinusoidal distribution 
of current density can be approximated by using a number of discreet 
wires in the design known as the saddle coil. A typical saddle coil 
(Figure 9.10(a)) consists of four linear segments and four 120° circular 
arcs on a cylindrical surface. This design approximates the sin <j> distri- 
bution by having two linear segments at 60° and 120° carrying equal 
currents in one direction and two other segments at 240° and 300° 
carrying equal currents in the opposite direction. The sinusoidal 
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Figure 9.10. (a) Saddle coil; (b) birdcage coil. 

current density can be more accurately approximated by using the so- 
called birdcage coil, which includes a large number of closely spaced 
linear segments (e.g., 16) uniformly distributed on a cylindrical surface 
[16]. Bi homogeneity in the birdcage coil (Figure 9.10(b)) is much 
higher than in a saddle coil of comparable size. It is interesting that 
the r.f. field in a birdcage coil can be made circularly polarized by 
simultaneously exciting two orthogonal modes corresponding to sin <f> 
and cos^ distributions of current. When used for signal reception, 
these modes can be combined together after one of them is phase 
shifted by 90 degrees. This doubles the NMR signal while the noise 
increases only by a factor of \/2 because the noise signals associated 
with the two modes are not correlated. Consequently, the SNR 
increases by a factor of \/2 [17]. Application of circularly polarized 
field for signal excitation is discussed in the following section. 



Quadrature Excitation 

As discussed previously (see Chapter 1), excitation of nuclear 
magnetization in a static magnetic field B = kB can be performed by 
a linearly polarized r.f. field Bj = \B X cos ut, where i and k are unit vec- 
tors in the x- and z-directions, respectively. A linearly polarized field can 
be viewed as the sum of two components rotating around the z-axis at 
the same frequency but in the opposite directions. Since the transverse 
magnetization in a sample is predominantly created by the component of 
Bj, which rotates in the direction of Larmor precession of spins, only 
half of the r.f. power is used effectively for excitation. 



Alternatively, excitation of nuclear magnetization can be 
achieved by using a circularly polarized r.f. field rotating in the direc- 
tion of Larmor procession of spins. Since in this case, known as the 
quadrature excitation, all r.f. power is effectively utilized, the use of a 
circularly polarized magnetic field lowers the required r.f. power as 
well as the power deposited in tissue by a factor of 2. 
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A. PHASE-SENSITIVE DETECTION 

As mentioned earlier, the observed NMR signal typically consists 
of different components in the narrow range of frequencies u) ± 6w, 
where 8u> <C u> . The direct sampling of NMR signal is difficult because 
it requires a very high sampling rate during analog-to-digital conver- 
sion. To lower the required sampling rate the signal can be shifted 
down in frequency by ui prior to analog-to-digital conversion while 
the phases and relative amplitudes of different signal components 
within the bandwidth 2&j are preserved. In practice such reduction 
in the signal frequency is commonly achieved by using two identical 
circuits known as phase-sensitive detectors (shown as PSD! and PSD 2 
in Figure A). 

To explain the operation of these detectors, we consider a simple 
case when the observed signal, S, is given by 

S=S cos(u;t + 4>), (A.l) 

where S , uj and <j> are the amplitude, frequency and phase of the signal, 
respectively. Initially the signal is amplified by a low-noise amplifier. 
The amplified signal is then divided equally between the two phase- 
sensitive detectors. In PSD X the incoming signal is multiplied by 
cosu t. The output signal, S u can therefore be expressed as 

Si oc S cos[(u; - uo)t + <f>] + S cos[(w + Lo )t + <j>\. (A.2) 

Because the frequency of the NMR signal is close to lo (i.e., 

Iw-WqI <w ), the term S cos[(u> + u )t + <p\ describes the rapidly 
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Figure A. A diagram of the phase-sensitive signal detection. 



varying signal component at the frequency of u + tu ss 2u> . In contrast, 
the term S cos[(u> - co )t + <f>] describes the slowly varying component 
of the output signal. The rapidly varying component is eliminated by 
low-pass filtering of S v 

In PSD 2 the incoming signal is multiplied by sinw £. The output 
signal, S 2 , is therefore given by 

S 2 oc -S sin[(w - iv )t + 0] + S B in[(w + w )i + 4>]. (A.3) 

Again the rapidly varying signal, proportional to sin[(a> + u) )t + <f>], is 
removed by low-pass filtering of S 2 . The resulting output signals Si 
and S 2 can be combined into one complex signal 

Si -jS 2 oc S Q e* u '- a <> )t+i *. (A.4) 

Let us now consider phase-sensitive detection of the NMR signal 
in two-dimensional MR imaging. Neglecting the effect of magnetic 
field inhomogeneities, the NMR signal from a sample can be written as 



Soc 



// 



M e~ 



-TE/T 2 



cos(u) t + k x x + k y y + 4>) dx dy, 



(A.5) 



where M is the equilibrium nuclear magnetization, TE is the echo 
time, k x and k y are dependent upon the readout and phase-encoding 
gradients, respectively. Taking into account equations (A.2)-(A.4) it 
is easy to show that after phase-sensitive detection and low-pass 
filtering, the complex output signal S 1 -jS 2 is given by 



Si -jS 2 = £ If M xy e ik * x+ik > y dxdy, 
where M xy = M e~ TE ^ T2+J<1 ' and £ is an arbitrary constant. 



(A.6) 
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B. IMAGE DISPLAY 

Images are typically displayed on the screen of an image device by 
using a number of screen elements (e.g., 512 x 512), referred to as dis- 
play pixels. The screen location of a pixel representing a given element 
from a two-dimensional data array is denned by the address of the 
element within the array. The pixel color is defined by the value of 
the element. There are several standard systems for reproducing 
color. One such system, known as the RGB (red, green and blue), is 
often implemented using 2 8 = 256 shades of each main color (i.e., red, 
green and blue). Consequently, the total color palette including 
shades of gray consists of 2 24 combinations of the main colors. 

Diagnostic MR images of various anatomical structures are nor- 
mally displayed by using a number (e.g., 256) of different shades of 
gray chosen according to the utilized relationship between intensity 
and brightness. Among countless possibilities for mapping intensity 
values to the displayed brightnesses, a linear relationship between 
intensity and brightness is used most often. In various instances 
image contrast can be enhanced if a limited range of intensity values 
is converted to the full range of brightnesses. In this approach inten- 
sity values smaller than an arbitrary chosen lower limit are typically 
represented by the black color, intensity values greater than an arbi- 
trary chosen upper limit are represented by the white color and the 
intermediate intensity values are displayed by using different shades 
of gray. Usually the upper and lower limits are adjusted manually or 
automatically to enhance contrast between various image structures. 
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